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lenglng  the  validity  of  l 


■ nonlinearity,  and  sorption-desorption  nonsingularity. 


itrlx)  are  probably  t! 


transfer  limitations  associated  with  aggregated  media  and/or  sorbent 
organic  matter  appear  to  predominate,  whereas  aquifer  heterogeneity 
appears  to  be  dominant  at  the  field-scale. 

Experimental  methods  for  investigating  the  sorption  kinetics  of 
organic  chemicals  were  developed  and  evaluated.  These  methods  were 

sorption.  A log-log  Inverse  relationship  was  discovered  between  the 
kinetic  sorption  constant  and  the  equilibrium  sorption  constant;  this 
relationship  was  quantified.  The  influence  of  sorbate  molecular 
structure  on  sorption  kinetics  was  recognized.  The  influence  of  organic 


investigated;  experiments  revealed  that  the  kinetic  sorption  constant  is 
log-linear ly  correlated  to  the  fraction  of  cosolvent.  A model 
explaining  the  influence  of  organic  cosolvent  on  sorption  kinetics  was 
developed.  A model  for  simulating  solute  transport  under  conditions  of 
multiple  sources  of  nonequilibrium  was  developed  and  evaluated.  A model 
simulating  solute  transport  at  the  field  scale,  where  nonideallcy  may 


INTRODUCTION 


1C  transport 
groundwater  remediation  technology  has 


groundwater  contamination,  research  on 

rapidly.  Knowledge  of  the 

essential  for  such  research.  Of  the  various 
of  the  most  important:  it  can  have  profound 

Investigation  of  sorption  nonideality  is  a relatively  new  and 
multi-disciplinary  endeavor.  * large  body  of  literature  is  developing, 

analytical  chemistry,  environmental  engineering,  and  civil  engineering, 
in  addition  to  soil  science  and  groundwater  geology/hydrology. 

The  conventional  equation  used  to  describe  advective-dispersive, 
one-dimensional,  steady-state  transport  of  a sorbing  solute  is  (Lapldus 


" 8x2 


Is  dispersion  coefficient 


(1-1) 
■s  sorbed-phase 


average  pore-water  velocity  (I/T),  , is  bulk  density  (H/L1) , I is  soil- 
water  content  (L3/La),  x is  distance  (L),  and  t is  time 


(T) . Describing 


Che  3S/3C  Cera  is  ofcen  Che  crlcics 
typical  approach  has  been  Co  assume 


seep  in  model  formulation . 


is  Kp  8C/3C.  Upon  subscicucion,  Che  cransporc  equation  becomes 

R“  - - vfx  (M) 

where  R - (1  + (p/9)  Kp)  is  Che  recardacion  factor  (Kp  - equilibrium 

Cquacion  (1-2)  may  be  written  in  Che  following  non-dimensional  form 
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(relative  Co  Che  mean  velocity  of  water).  The  solution  Co  equation  (1- 

as  P Increases,  Che  breakthrough  curve  (BTC)  approaches  a symmetrical, 
sigmoidal  shape.  However,  for  P less  Chan  approximately  10,  Che  BTC  is 


The  assumptions  employed  Co  develop  eq.  1-1  simplify  transport 
analysis  and  may  be  adequate  under  certain  conditions.  Daca  exhibiting 
behavior  that  deviates  from  that  predicted  by  Che  simple  model  have  been 
reported,  however,  challenging  the  validity  of  this  approach.  This 
nonideal  sorptive  behavior  has  been  attributed  to  several  different 
factors,  including  kinetic  sorption  reactions,  diffusive  mass  transfer 
resistances,  isotherm  nonlinearity,  and  sorptlon-desorpcion 

1-3  are  not  met,  the  A-D  model  will  not  fully  describe  solute-transport 
daca.  In  these  situations,  the  data  may  be  characterlted  as  exhibiting 
nonideality.  As  a majority  of  the  past  and  present  groundwater 

m modeling  efforts  have  employed  Che  A-D  model  and  its 
t assumptions,  it  Is  important  to  analyse  the  validity  of  this 
approach.  For  cases  where  the  A-D  model  is  not  valid,  the  effects  of 
nonideality  on  groundwater  contaminant  transport,  modeling,  and 
remediation  need  to  be  assessed. 


The  literature  is  replete  with  reports  on  the  Influence  of  soil 
structure  on  solute  transport.  Such  reports  dace  back  to  over  100  years 
ago  (Schumacher,  1864;  lawes  at  al. , 1882) , as  noted  by  White  (1985) . 

the  observed  nonideal  behavior  (e.g.,  "preferential 


The  explanation 

flow,"  "macropore  flow,"  "channeling, 

■partial  displacement,"  "subsurface  storm  flow")  is  the  result  of  ni 
unimodal  pore-size  distributions  (l.e.,  non-uniform  velocity  field) 


Rao.  1969)  is  often  employed.  Uieh  this  approach  Che 
solute  transport  occurs  by  advectlon  and  dispersion,  am 


years,  when  it  was  used  by  Gardner  and  Brooks  (1956) . Early 
mathematical  models  employing  this  conceptualization  were  presented  by 
Skopp  and  Warrick  (1976)  and  van  Genuchten  ec  al.  (1974). 

Rapid  transport  in  the  mobile  domain  is  accompanied  by  diffusive 
mass  transfer  of  solutes  between  the  mobile  and  immobile  domains,  which 
results  in  the  latter  behaving  as  sink/source  components.  Solute 
transport  in  such  syscoms,  as  described  by  breakthrough  curves  (BTCs), 
is  characterized  by  early  initial  breakthrough  and  by  ■tailing"  or 
delayed  approach  to  relative  concentration  values  of  0 or  1.  Because 
access  to  some  portion  of  the  porous  medium  is  constrained  by  diffusive 
mass  transfer,  solute  in  the  system  may  be  considered  to  be  in  a state 

transport  nonequilibrium  (THE).  One  of  Che  first  to  suggest  that  "dead 


water"  could  cause  nonideal  BTCs  was  Danckwerts  (1953). 

Behavior  attributable  to  THE  has  been  observed  in  aggregated, 
heterogeneous  (with  respect  to  hydraulic  conductivity),  and  fractured 
porous  media  os  well  as  in  macroporous  media.  Some  of  the  earliest  work 
on  solute  transport  in  aggregated  soils,  where  the  nonldealicy  observed 
in  experimentally  derived  BTCs  was  attributed  to  THE,  was  performed  by 
Biggar  and  Hielsen  (1962a)  for  nonsorbing  solutes  (Cl") , by  Bigger  and 


Nielsen  (1963)  for  sorbing  inorganic  soluees  (Mg3*) , and  by  Elrick  ec 
al.  (1966).  Kay  and  Elrick  (1967).  Davidson  ec  al.  (1968).  Davidson  and 
Chang  (1972),  and  Green  eC  al.  (1968,  1972)  for  organic  compounds 


experiments  were  all  performed  wich  packed  soil  columns,  and  while 
demonstrating  the  influence  of  soil  aggregation  on  solute  transport, 

atcempc  to  more  closely  simulate  field  conditions,  miscible  displacement 
studies  have  been  performed  using  undisturbed  soil  cores.  Early 
experiments  of  this  nature  demonstrating  nonideal  soluce  transport  of 
Cl  , which  was  attributed  to  macropore  flow,  were  reported  by  Elrick  and 
French  (1966) . Numerous  studies  of  a similar  nature  have  since  been 
reported  for  Cl"  and  ocher  nonsorbing  solutes  (see  White,  1985;  Beven 
and  Germann,  1982) . Similar  work  has  been  reported,  although  to  a much 
lesser  extent,  for  sorbing  Inorganic  solutes  (c.f.,  Souchworth,  1987; 
Jardlne  ec  al.,  1988)  and  organic  solutes  (White  et  al.,  1986; 

Southworth  ec  al.,  1987). 

Several  field  studies  of  solute  transport  under  controlled 
conditions,  where  observed  nonideality  was  attributed  to  TNE,  have  been 

Bablker,  1976).  An  early  study  of  TNE  in  field-scale  transport  of  an 
organic  solute  (herbicide  plcloram)  was  performed  by  Rao  et  al.  (1974) . 
Other  studies  for  organics  (Bowman  and  Rice,  1986;  Jury  at  al..  1986a) 


a hicontinuum  concept  may  also  be  applied  to  heterogeneous  porous 
such  as  aquifers  consisting  of  laminae  of  varying  hydraulic 


conductivities  (K) . In  this  case,  the  high-K  layet(s)  correspond(s)  Co 
the  "mobile  zone"  and  Che  low-K  layer(s)  to  the  "immobile  zone"-  Solute 
dispersion  is  affected  by  mass  transfer  between  Che  high-  and  low-K 
layers.  The  mass  transfer  mechanism  may  be  explained  as  follows: 
solute  adveccion  is  greacer  in  the  high-K  chan  in  the  low-K  layers,  and 
as  a result,  differential  solute-front  advancement  occurs  in  the 
different  layers,  which  creeces  incer-layer  concentration  gradients. 
Hence,  incer-layer  diffusive  solute  transfer  occurs,  which  results  in 
THE.  A diffusion-based  model  of  this  type  was  presented  by  Gillham  et 
al.  (1984),  while  Golcz  and  Roberts  (1989)  applied  a mass-cransfer-based 
model  to  data  obtained  from  a natural-gradient  field  experiment. 


Transport  in  fractured  media  can  also  be  modeled  with  a 
approach.  The  fractures  serve  as  zones  of  preferential  transport,  with 
the  inter-fracture  matrix  acting  as  a diffusional  sink/source  component 
(e.g.,  Grlsak  and  Pickens,  1980;  Lever  ec  al.,  1985). 

dispersive  transport  from  bulk  solution  to  the  boundary  layer  (i.e., 
adsorbed  water  surrounding  the  sorbent) ; (2)  diffusive  transport  across 
the  adsorbed  water  (i.e.,  film  diffusion);  and  (3)  pore  and/or  surface 

relative  to  the  other  two  and  may  therefore  be  designated  as  the  primary 


rate-limiting  step. 


thus  minimizing  Che  importance  of  mass  transfer  resistances  wiehin  these 
regions.  Ic  appears  that  film  diffusion  controls  the  solute  uptake  race 

1973;  Weber  and  Liu,  1980;  HcKay  et  al.,  1987).  However,  incra- 
aggregace  diffusion  appears  to  concroi  solute  uptake  for  Che  majority  of 
the  reaction  period  (Leenheer  and  Ahlrlchs,  1971;  Furusawa  and  Smith, 
1973;  Khan,  1973;  Ruchven,  1984;  Crittenden  et  al.,  1986;  HcKay  ec  al., 
1987;  Roberts  ec  al.,  1987).  Intra-aggregate  diffusion,  therefore,  is 
usually  Che  transfer  process  designated  as  Che  rate-limiting  step.  Note 
that  sorpcion  at  the  sorbent-liquid  interface  is  assumed  to  be 
instantaneous;  the  rate  at  which  the  solute  is  transported  to  and  from 


al.  (1986)  Co  assess  the  relative  contribution  of  various  mass  transfer 

Nkedl-Kizza  et  al.  (1982)  for  ’HjO,  C1‘,  and  Cm*  transport  in  packed 
columns  of  an  aggregated  soil.  Roberts  et  al.  (1987)  concluded  Chat  the 
contribution  of  internal  mass  transfer  resistance  to  dispersion  was 
three  to  ten  times  greater  than  that  of  external  mass  transfer 


sensitivity  analysis . The 


exhibited  by  nonsorbing,  i 


i well  es  sorbing,  solutes.  This  highlights 
its  responsible  for  this  nonequilibrium  is  a 

nvolved.  If  the  definition  of  sorption  is  limited  to  the  processes 
courting  at  the  solid-liquid  interface,  it  is  besc  Co  label  this  case 
s nonequilibrium  transport  as  done  above. 


soil-water  saturation.  Phenomena  similar  to  chose  discussed  above  have 
also  been  observed  in  unsaturaced  systems.  Biggar  and  Nielsen  (1960) 
performed  miscible  displacement  studies  under  unsacurated  conditions  and 
observed  "progressive  shifting  of  the  BTC  to  Che  left"  as  soil-water 


The  Influence  of  unsacurated  conditions  upon  solute-transport 

Uniform  distribution.  Nonideallcy  will  generally  Increase  when 
soil-water  content  is  reduced  from  saturation  for  a medium  that  has  a 
relatively  narrow,  uniform  distribution.  BTCs  obtained  under  saturated 
conditions  will  exhibit  a relatively  small  amount  of  dispersion  because 
of  the  uniform  distribution.  Upon  desaturation,  some  fraction  of  the 
pores  no  longer  readily  transmit  water.  These  pores  thus  become 


"Immobile"  domains  choc  acc  as  dlscrlbueed  dlffusional  sink/source 
compononcs  for  soluce.  BTCs  obcalned  under  unsacuracod  condietons  will 
exhlbic  a greaCer  degree  of  dispersion  and  may  be  shtfced  leftward  In 
comparison  Co  Che  BTC  obcalned  under  sacuraclon  because  of  Che  Impacc  of 
Che  immobile  domains.  Such  behavior,  uhere  nonidealicy  increases  as 
soil-water  concenc  Is  reduced  from  sacuraclon,  has  been  reporced  by 
several  researchers  (c.f.,  Bigger  and  Nielsen,  1960;  Krupp  and  El  rick. 
1968;  Caudec  ec  al.,  1977;  DeSmedc  and  Blerenga,  1984). 

;er  concenc  can  resulc  in  Che  effecclve 


(i.e.,  cransmiccing)  pi 
Cransporc  nonidealicy  cl 


ze  dlscribucion  oi 


nonideal,  in 


: Che  larger  soil-w 


as  soil-wacer  concenc 
> cridcal  value  of  sol 

soil-wacer  concenc  is  furcher  reduced.  Such  behavior  was  reporced  by 
Krupp  and  Elriek  (1968)  and  is  exhibiced  by  Che  daca  of  DeSmedc  and 
Blerenga  (1984). 

Nonuniform  dlscribucion.  Soluce  cransporc  under  sacuraced 
conditions  in  porous  media  characcerlzed  by  a wide  pore  size 
dlscribucion  (e.g.,  presence  of  macropores)  will  exhlbic  nonideal 
behavior,  as  discussed  above.  The  macropores  drain  upon  desacuracion, 
which  resulcs  in  a narrowing  of  Che  effecclve  pore-size  dlscribucion. 
Soluce  cransporc  nonidealicy  will  chi 


soil-wacer 


nd  French  (1966)  and  Seyfrled  ar 


It  is  likely  that,  once  the  largest  pores  are  drained,  che  porous 
medium  will  thereafter  exhibit  behavior  similar  to  that  discussed  in 
case  (1),  Further  reduction  in  soil-water  concent  will  create  ijmnobile 
regions  and  result  in  increased  nonideality.  After  the  critical  soil- 
water  content  is  reached,  further  reductions  will  eliminate  che  immobile 
regions,  resulting  in  a reduction  in  nonideality. 


d immobile  regions  ci 


explicitly,  with  the  use  of  Flck'a  law  to  describe  the 
physical  mechanism  of  diffusive  transfer; 
explicitly,  with  the  uso  of  an  empirical  first-order  mass- 
transfer  expression,  in  place  of  the  mechanistic  expression  of 
(1),  to  approximate  solute  transfer;  and 

implicitly,  with  che  use  of  an  effective  or  lumped  dispersion 
coefficient  that  includes  che  effects  of  slnk/source 
diffusion,  as  well  as  hydrodynamic  dispersion  and  axial 
diffusion.  This  lumped  parameter  replaces  the  usual 
hydrodynamic  dispersion  coefficient  (D)  in  the  advective- 
dispersive  transport  equation.  These  three  approaches  will  be 


diffusion  model  took  place 
Edeskucy  and  Amundson,  1952;  Faster 
1966).  Rao  and  co-workers  applied 
nonequilibrium  in  soils  (Rao  ec  al. 


Early  development  of  che 


le  models  to  describe 


Klzza  ec  al . , 1982) . Wu  and  Gschwend  (1986)  extended  this  approach 
study  the  sorption  kinetics  of  organic  solutes. 

negligible.  If  this  is  not  the  case,  inclusion  of  a film-transfer 


presented  by  Crittenden  ec  al.  (1986)  and  Miller  and  Webor  (1986). 
Solute  transfer  from  mobile  to  immobile  regions  Is  described  by  a t 

particle  (i.e.,  film  transfer)  followed  by  diffusion  within  Che 


The  diffusion-based  models  discussed  above  were  developed  assuming 
a spherical  geometry,  A conceptually  similar  model  employing  a 
cylindrical  geometry  was  presented  by  van  Genuchcen  at  al.  (198&)  to 

The  increased  accuracy  obtained  by  Including  the  actual  solute- 
transfer  mechanism  In  the  A-D  model  is  countered  by  a concomitant  - 
increase  in  model  complexity.  To  develop  the  diffusion  equation,  a 
geometric  description  of  the  porous  medium  is  required.  Most  models 
have  assumed,  primarily  for  mathematical  simplicity,  a spherical, 
uniformly  sized  aggregate  structure,  natural  porous  media,  however,  may 
be  far  from  this  idealization.  For  the  above  models  to  be  applicable  to 
more  than  Just  a few  ideal  cases,  aggregate  size  distributions  and 
variations  in  aggregate  shape  must  be  accounted  for.  In  an  attempt  to 
apply  the  simplified  model  to  nonideal  cases,  Rao  at  al.  (1982)  assessed 
whether  solute  diffusion  from  non-spherical  aggregates  could  be 
approximated  by  diffusion  from  equivalent  spheres.  They  showed  that 


represented  by  equivalent  spherical 


aggregates  whose  radii  are  such  that  tho  sphere  volume  is  equal  Co  Che 

diffusion  model  developed  for  spherical  aggregates  (1982)  to  slabular 
and  cylindrical  shapes  (1985a) . A form  factor  was  employed  to  adjust 
the  intra-particle  diffusion  equation  to  the  selected  shape.  It  was 
shown  that  one  aggregate  shape  could  be  approximated  with  another  by 
employing  the  same  ratio  of  external  aggregate  surface  area  to  total  bed 
porosity.  The  approximations  gave  nearly  identical  BTCs  for  short  and 
long  contact  times,  and  relatively  minor  deviations  for  intermediate 


A method  that  extends  Che  physical  diffusion  model  to  more  general 
conditions  involving  aggregates  with  non-spharical  geometries  has  been 
developed  by  van  Genuchcen  (1985).  The  method  is  based  on  Che  use  of  a 
geometry-dependent  shape  factor  (f)  that  can  be  used  Co  transform  an 


aggregate  of  a given  shape  and  size  into  an  equivalent  sphere  with 
similar  diffusion  characteristics  as  the  original  aggregate.  While  the 
shape  transformation  was  found  to  be  accurate,  it  is  limited  to 


uniformly  sized  aggregates,  and  no  provisions  are  given  to  deal  with 
mixtures  of  aggregate  sizes.  The  effect  of  various  aggregate  shapes  on 
sorption  was  assessed  by  Rasmuson  (1985a).  For  early  times,  solute 
uptake  was  nearly  identical  for  spherical,  cylindrical,  and  slabular 
aggregates.  For  larger  times,  the  uptake  race  vas  greatest  for  the 
slabular  and  lowest  for  the  spherical  aggregates.  The  slabular  and 
cylindrical  aggregates  have  solute  uptake  races  greater  chan  the 
spherical  because,  for  equivalent  characteristic  lengths  (i.e.,  equal 


greater  chan  chat  of  a sphere's  (6,  S.7,  < 
result  of  greater  uptake  races,  a delay  It 


lute  breakthrough  (i.e., 
n tailing  would  be  expecced 


.r  aggregate  systems. 


i slngl 


radii  of  each  aggregate-sire 

ixed  aggregate-sized  medium, 
validated  by  Nkedl-Kizza  et 


>r  cylindrical 

aggregate-sir 
variations,  n 

approximation  of  diffusion  in  Che  actual 
The  viability  of  this  approach  was  furthe 
al.  <1982),  where  a distribution  of  aggregate  sizes  (0.2  to  0.47  cm  it 
diameter)  and  shapes  was  well  represented  by  an  equivalent  spherical 
aggregate.  Several  forms  of  an  aggregate-size-distrlbucion  weighting 
factor  were  derived  by  Ivory  (1983). 

Combining  the  shape  transformations  of  von  Genuchten  (1985)  or 
Rasmuson  (1985a)  with  the  aggregate  size-distribution  transformation  c 
Ivory  (1981)  or  Rao  et  al.  (1982),  a porous  medium  consisting  of 
different  sizes  of  aggregates  having  various  non-spherical  shapes  can 

spherical  aggregate  of  an  equivalent  uniform  size.  In  this  manner,  r 
physical  diffusion  model  may  be  extended  to  many  nonideal,  field 


Aggiregate-size  distribution  (ASD)  was 

modeled  by  assuming  a progressively  swelling  hypothetical  aggregate 

swelling  representative  aggregate  is,  in  effect,  an  effective 
characteristic  aggregate,  whose  sire  Increases  with  time.  The  size  of 
the  effective  aggregate  is  time-dependenc  because  of  differencial  solute 
uptake  races  for  different-sized  aggregates,  Ac  early  times,  smaller 
aggregates  dominate  uptake;  therefore,  the  effective  aggregate  Is  small. 
As  time  increases,  the  larger  aggregates  become  more  important,  which 
results  in  an  increase  in  the  effective  aggregate  size.  The  effective 

solution.  Knowledge  of  the  size-distribution  function  or  of  the  average 
aggregate  size  is  not  required.  It  is  not  possible,  however,  to  derive 
an  absolute  value  for  the  intra-particle  diffusion  coefficient. 

A model  that  directly  Includes  in  the  transport  equation  provisions 
for  an  aggregate-size  distribution  has  been  developed  by  Rasmuson 
(1985b) . Any  discrete  distribution  is  allowed,  and  the  following 
parameters:  aggregate  shape,  intraaggregate  diffusion  coefficient,  film 
mass-transfer  coefficient,  and  sorption  equilibrium  constant  may  vary 
for  each  size  class.  The  size  distribution  is  accounted  for  by  the 
employment  of  the  parameter  Fib),.,  the  volume  fraction  of  aggregate 

aggregate  size-class  i to  the  total  aggregate  volume).  In  effect,  the 
influence  of  each  size  class  on  solute  transport  is  weighted  by  its 


aggregate-size  distributions  and  aggregate-shape  variations  to  be 

The  effect  of  aggregate-size  distribution  on  solute  transport  has 
been  assessed  by  Moharlr  et  al.  (1980b)  and  Rasmuson  (1985b).  ASD  was 

. C/C.-1).  The  latter  effect  is 

arly  breakthrough)  is  a result  of  the  increased  number  of  smaller-than- 

surface  area.  This  results  in  an  Initial  rate  of  solute 
uptake  greater  than  average,  which  results  in  a delayed  breakthrough. 

average-size  aggregates.  Because  of  their  increased  individual  volumes 

saturation  more  slowly  than  the  average.  Therefore,  complete 
breakthrough  is  delayed.  The  ASD  effect  on  transport  was  shown  to  be  of 

distributions,  and  ASDs  skewn  toward  smaller  aggregates 


(Moharlr  el 


in  the  above  ASD  effects,  t: 


i existence  of  ai 


when  none  is  assumed  would  resule  in  Che  physical  nonequilibrium  model 
under-predicting  tailing  and  over-predicting  early  breakthrough.  The 
effects  of  ASD  may  be  ignored  (i.e.,  Che  porous  medium  may  be  created  as 
an  assemblage  of  uniformly-sited  aggregates)  when  the  site  distribution 
is  sufficiently  narrow.  Ruthven  and  Loughlin  (1971)  determined 
mean/scandard  deviation  > 10  as  being  sufficiently  narrow.  Cooney  ec 
al.  (1983)  found  an  ASD  range  of  less  than  100  mesh  to  be  sufficiently 


cases  of  aggregates  having  different  sorpcive  properties.  For  example, 
sorbents  with  different  intraparticle  diffusion  resistances  or  different 

Ehyalsal  mass-transfer-  Recognizing  the  uncertainty  surrounding 
pore  geometry  analysis,  the  diffusion  model  may  be  simplified  by 

•C  description  of  diffusive  transfer  with  a 
s transfer  expression.  Actual  description  of  the 
•re  is  then  no  longer  required.  Solute  transfer 


kinetic  first-order 


s transfer  coefficient  and  is  assumed  to 
average  immobile-region  concentration  is 


assumed  (i.e.,  a perfectly  mixed  region).  This  assumption  is  in 


Early  devel 


physical  mass-transfer  (or  •two-region") 


presented  by  Dean  (1963)  and  Coats  and  Smith 
was  extended  to  sorbing  solutes  by  van  Genuchten  and 


Wierenga  (1976).  This  nodal  is,  in  essence,  Che  linear  driving  force 

of  diffusion  inco  spherical  particles.  This  approximation  vas  first 

Through  the  transformation  technique  mentioned  previously,  van 

Expressions  for  Che  mass-transfer  coefficient  in  terms  of  physical 

1984) , and  Parker  and  Valocchi  (1986) . With  Che  use  of  the  above 
relationships,  Che  otherwise,  empirical  mass-transfer  coefficient  can  now 
be  expressed  in  terms  of  measurable  physical  parameters.  This  allows 

transfer  model  (parameters  obtainable  only  through  curve  fitting)  may  be 


Effective  dispersion  model.  The  third  physical  nonequilibrium 

conventional  A-D  transport  model.  This  effective  dispersion  coefficient 
implicitly  accounts  for  any  diffusive  transfer  between  mobile  and 
immobile  regions.  Early  work  in  soil  science  employing  this  approach 
was  presented  by  Passioura  (1971)  and  Passioura  and  Rose  (1971). 

BTC  will  increase  to  a point  where  asymmetry  is  difficult  to  discern. 

The  effective  dispersion  model  will  fit  experimental  data  especially 


The  three  physical  models  provide  a range  in  complexity  and 
conceptual  accuracy.  It  has  been  shown  that  both  the  physical  mass- 
transfer  and  the  effective  dispersion  models  are  degenerate  cases  of  the 
physical  diffusion  model  (Raacs,  1981).  Those  approximations  will  be 


use  is  dependent  upon  Che  extant  conditions,  and  upon  the  accuracy 
provided  in  relation  to  their  relative  complexity. 

The  performances  of  the  diffusion  and  mass-transfer  models  have 
been  compared  by  Rao  et  al.  (1980b).  'Gold  (1986),  Gold  and  Roberts 
(1986),  and  Hiller  and  Haber  (1986).  In  all  four  cases,  Che  diffusion 
model  gave  similar  or  only  slightly  better  results.  Van  Cenuchten 
(1985)  assessed  Che  ability  of  the  physical  mass-transfer  model  to 
approximate  solutions  given  by  Che  physical  diffusion  model.  BTCs 
predicted  by  the  two  models  for  different  values  of  the  dimensionless 
mass-transfer  coefficient  were  compared.  The  mass-transfer  model  gave 
good  approximations  when  values  of  Che  dimensionless  mass-transfer 
coefficient  were  greater  than  10.  The  physical  mass-transfer  model 
appeared  to  give  results  virtually  identical  to  a physical  diffusion 
model  when  the  following  condition  was  met  (Parker  and  Valocchi.  1986): 
(D.  L/v  R„  b2)  > 0.3 

where  D.  is  lntra-aggreg.c.  diffusion  coefficient.  L 1,  length  of 
interest,  v is  average  pore-water  .velocity,  R,„  is  intra-aggregate 
retardation  factor,  and  b is  aggregate  radius.  Raghavan  and  Ruchven 
(1985)  found  that  the  mass-transfer  model  gave  good  approximations  of 
the  diffusion  model  for  (Da  i/ba  v)  > 10. 


From  the  Foregoing,  it  appears  char  Che  increase  in  conceptual 
accuracy  gained  with  Che  diffusion  model  may  noc  Justify  Che  additional 
mathematical  complexity  and  increased  number  of  parameters  a 
with  the  model.  Hence,  the  physical  mass-transfer  mod 


is  experimental  BTCs.  In  an 


Che  extensive  tailing  exhibited  by  sod 
attempt  to  account  for  this  tailing,  i 
proposed  (Brusseau  and  Rao,  1989).  Extensive  tailing  is  hypothesized  to 
result  from  an  additional  intraaggregate  diffusion  mechanism.  The 

This  dual  intraaggregate-poroslty  model  was  constructed  by  combining 

todels  presented  above  with  the  bldlsperse- 
il  engineering.  The  resultant  is  a porous 
ie  porosity  regions:  macropore  (interaggregate 
raaggregate  macroporosity) . and  micropore 


particle  model  used  in 
porosity) , mesopore  (intr 


developed  by  dividing  t! 
transport  is  at 


gate  porosity  model  discussed  above,  the 

ons.  *A  three-region  model  may  also  be 

i et  ml.  (1986).  Advective-dlspersive 
in  both  the  "rapidly  mobile"  and 


"moderately  mobile"  regions,  while  mass  transfer  is  assumed  to  occur 
between  the  rapidly  and  moderately  mobile  and  the  moderately  mobile  and 


The  blmodal  characterisation  of  the  mobile  (Morisawa  et  al.,  1986) 
or  Immobile  (Brusseau  and  Rao,  1989)  domain  employed  by  the  above  models 
is  in  direct  contrast  to  Che  unimodal  characterization  employed  by  the 
biconcinuum  models.  It  is  very  likely  that  a distribution  of  pore 
sizes,  rather  than  one  distinct  size,  exists  within  en  aggregated  soil, 

conditions.  It  may  be  difficult,  however,  to  quantitatively  discretize 
the  mobile  or  immobile  porosity  into  two  distinct  fractions.  In  any 

media  whose  lncraaggregate  porosity  Is  not  extremely  large.  The 
blcontlnuum  (i.e.,  two-region)  models,  therefore,  are  probably 


Obviously,  the  two-  and  three-region  models  are  simple 
approximations  of  true  porous  media.  Higher-order  models  could  be 
developed  by  progressively  Increasing  the  number  of  explicitly  defined 
regions  (i.e.,  H regions  and  N-l  mass  transfer  resistances). 

Ultimately,  this  discretization  will  lead  to  a continuous  pore-size  (or 
pore  velocity)  distribution  function.  For  example,  a transfer  function 
model  (Jury,  1982;  Jury  et  al.,  1986b)  considers  a continuous  pore-water 
velocity  distribution.  Pore  size  distribution,  estimated  using  a 
capillary  bundle  model,  has  been  utilized  in  the  past  to  model  soluce 


significant  amount  of  published  data  where  nonideal  behavior  is 
exhibited  in  the  absence  of  TNE;  for  example,  nonsorbing  solutes 
exhibiting  symmetrical  BTCs  whereas  sorbing  solutes  do  not  (cf. , Elrick 

al.,  1987;  Bouchard  et  al.,  1988;  Lee  et  al.,  1988).  Such  behavior 
suggests  that  the  nonequilibrium  is  a result  of  a sorption-related 
mechanism  (Brusseau  and  Rao,  1989).  Sorption  nonequilibrium  (SHE)  can 
result  from  either  chemical  nonequilibrium  or  lntrasorbent  diffusion 
(Brusseau  and  Rao,  1989).  Chemical  nonequilibrium  results  from  rate- 
limited  interactions  between  the  solute  and  specific  sorption  sites  of 
the  sorbent.  Nonoquilibrium  resulting  from  lntrasorbent  diffusion 
Involves  the  diffusive  mass  transfer  of  sorbate  within  the  sorbent 
matrix.  Although  both  TNE  and  lntrasorbent  diffusion  Involve  a 
diffusive  mass-transfer  mechanism,  they  are  considered  different  because 
TNE  is  predominantly  a pore-diffusion  process,  while  lntrasorbent 
diffusion  is  more  similar  to  a solid-diffusion  process.  The  results  of 
several  recent  experimental  and  theoretical  analyses  (Bouchard  at  al. 
1988;  Brusseau  and  Rao.  1989;  Nkedl-Kixxa  et  al. , 1989)  suggest  that 
sorption  nonequilibrium  for  organic  chemicals  may,  in  many  cases,  be 
controlled  by  rate-limited  interactions  between  the  sorbate  and  sorbent 


if  nonequilibrium  at 


organic  matter.  It  is  important  tc 

dispersion  was  Wilson  (1940),  who  stated  "the  width  of  a band 
(chromatographic  pulse)  may  Increase  because  the  leading  edge  of  Che 
band  migrates  coo  rapidly  on  accounc  of  a low  rate  of  adsorption,  or 
because  the  trailing  edge  of  Che  band  migrates  Coo  slowly  on  accounc  of 

sorption  kinetics  involved  simple  "one-site-  first-order  models  in 
which  Che  sorption  race  is  taken  as  a function  of  Che  concentration 
difference  between  the  sorbed  and  solution  phases  (cf . , Oddson  ec  al. 
1970) , This  approach  was  based  on  Chat  taken  by  researchers  in  chemical 
engineering  (o.g. , Lapidus  and  Amundson,  1952).  This  model  has  failed 
to  predict  experimental  data  (e.g.,  Schwartenbach  and  Wescall,  1981;  Rao 
and  Jessup,  1983;  Wu  and  Cschwend,  1986;  Boesten  and  van  der  Pas,  1988). 
Sorption  data  from  batch  experiments  have  been  found  to  exhibit  a 
two-stage  approach  to  equilibrium;  a short  Initial  phase  of  fast  uptake, 
in  extended  period  of  much  slower  uptake.  This  pattern  is 
all,  sorption  reactions  (Jaffe  and  Ferrara, 

1983;  Wu  and  Cschwend,  1986).  Generally,  roughly  50%  of  the  sorption 
occurs  within  the  first  few  minutes  to  hours,  with  the  remainder 
occurring  over  periods  of  days  or  months  (e.g.,  Karlckhoff,  1980; 

McCall  and  Agin,  1985;  Oliver,  1985;  Ball  and  Roberts,  1985;  Coates  and 


Elzerman,  1986). 


approximation  with  a bicontinuum  approach.  Models  formulated  with  tl 


conceptualization  have  been  variously  called  "two-site",  "two- 

in  these  models,  nonequilibrium  is  assumed  to  be  Che  result  of  a 
time-dependent  sorption  reaction  of  the  solute  with  the  sorbent. 

Sorption  may  be  represented  by  two  reactions,  occurring  either  in  series 
or  in  parallel,  with  the  latter  being  Che  more  common  biconcinuum 
conceptualization.  The  sorbent  is  hypothesized  as  having  two  sorption 
domains,  where  sorption  is  instantaneous  for  one  and  is  race-limited  for 

a first-order  rate  equation,  while  the  ocher  domain  is  represented  by  an 
equilibrium  isotherm  equation.  It  should  be  noted  that  as  long  as  one 

1980;  Karickhoff  and  Morris,  1985).  Also,  Che  equilibrium-controlled 
domain  can  be  replaced  such  that  both  domains  are  kinetic  controlled. 

first  made  by  Ciddings  and  Eyring  (1955). 

developed  by  Selim  ec  al . (1978)  and  Cameron  and  Klute  (1977) , while 
Selim  et  al.  (1976)  and  Karickhoff  (1980)  presented  models  where  both 

Co  be  of  Che  time-dependent  class  (Selim  ec  al.,  1976;  van  Genuchcen, 
1981).  The  two-site  model  has  generally  been  able  to  represent 


where  C - solution-phase  concentration  (M/L1) , S,  - sorbed-phase 

concentration  (M/M)  for  kinetic  sices.  Assuming  chat  the  Preundlich 
equation  is  applicable  Co  both  groups  of  sites,  the  total  sorbed 
concentration  at  equilibrium  (St)  is  given  by 


S,  - FK,C‘ 

S,  - (1-F)  K,  C“ 

The  distribution  between  Che  two  types  of  sices  is  described  by  P. 
defined  as 


P - Kfl  / K; 

i is  the  sorption  partition  coefficient  for  Che  equilibrium 
" / M")  and  K,  is  Che  sorption  partition  coefficient  for  the 
lls  (Is"  / M"). 

» following  nondlmenslonal 


!f  * <•**>»“  if  * <w».  % - 

P 0“  ■ ax  (1-4) 

a-/l)R,|f-  mjc-i-s-)  (!— 5) 


hy  defining  the  following  dimensionless  parameters: 


1-F)  K,  c; 

i/»>  K.  Cf1 


“ - kj  (1-#)  RL/v 
/»-[!  + ((■/»)  F Kj  Cf'l  / 8. 

where  kj  is  the  ness  transfer  rate  coefficient  (1/T) , and  n is  the 
Freundllch  nonlinear  sorption  isotherm  coefficient.  Note  that,  in 
contrast  to  many  earlier  models  which  have  employed  linear  Isotherms 
(e.g. . van  Genuchten.  1981;  Nkedt-Kizza  et  el.,  1984).  che  Freundlich 
isotherm  is  used  to  represent  sorption  equilibrium.  Thus,  isotherm 
nonlinearity  can  also  be  accounted  for  with  the  equations  presented 
above.  In  the  model  formulation  of  Selim  et  al.  (1976),  which  was  in 
dimensional  form,  nonlinearity  was  included  as  well. 

The  discretisation  of  the  sorbent  surface  into  two,  parallel 
classes  of  sices  differing  in  reaction  time  may  be  extended  to 

model.  The  limiting  case  would  be  a continuous  distribution  of  sices 
and  associated  rate  constants,  A model  describing  this  case  has  been 
developed  by  Vlllermaux  (1976) , where  the  site  population  is  represented 
by  the  transfer  time  distribution  (i.e..  race  constant  distribution). 

The  models  employing  Che  diffusion  equation  are  intellectually 
pleasing  because  they  employ  a truly  mechanistic  approach.  These  models 
have  been  successfully  used  for  simplified  systems.  However,  once  real 
soil  is  considered,  with  all  its  complexity,  the  usefulness  of  these 


models  becomes  questionable.  For  example,  Co  apply  these  models  Che 
size  and  shape  distribution  of  the  structures  (e.g. . aggregates)  are 
required.  Determination  of  Chese  distributions,  especially  at  the  field 

It  must  be  stressed  that  these  models  have  been  developed  Co 
represent  systems  Influenced  by  THE.  Application  of  these  models  to 
systems  influenced  by  processes  ocher  than  TOE  can  result  in  erroneous 
results.  For  example.  Miller  and  Weber  (1986)  applied  an  intra- 
aggregate  diffusiAn-based  model  to  the  transport  of  organic  solutes 
(lindane,  nitrobenzene)  through  packed  soil  columns.  The  intra- 
aggregate  diffusion  coefficients,  determined  by  optimization  of  data 
obtained  from  batch  experiments,  were  on  Che  order  of  10_u  cmz/sec. 

These  values  are  several  orders-of-magnitude  smaller  Chan  those  expected 
for  typical  soils  (e.g.,  10"*-10-7).  Such  a discrepancy  suggests  that 
the  observed  nonequllibrium  Is  the  result  of  some  process  other  than 
simple  diffusion  into  soil  aggregates . A similar  application  of  a 
diffusion-based  model  was  made  by  Hutzler  ec  al.  (1986),  »ho  modeled  the 
transport  of  chloride,  trichloroethene,  end  bromoform  through  columns  of 
packed  soil.  They  utilized  aggregate  radii 

measured.  This  again  suggests  that  a process  other  than  intra-aggregate 
diffusion  was  responsible  for  the  observed  nonequllibrium. 

The  physical  mass-transfer  (i.e. , two-region)  and  Che  two-site 
sorption  models  have  been  shown  to  be  mathematically  equivalent  for  Che 
n written  in  dimensionless  form  (van  Genuchten, 

, 1984) . Hence,  they  produce  similar  fits  Co 


s as  the  fitting  parameter; 


linear  isotherm  c.i 
1981;  Nkedi— Kizza 


experimental  data  (Rac 


fitting  models  to  experimental  data. 

the  processes  represented  by  these  parameters.  Such  models  would, 
therefore,  be  constrained  to  operation  in  a calibration  mode. 

transport  under  nonequilibrium  conditions.  Each  has  advantages  and 
disadvantages.  The  choice  of  model  is  constrained  by  available  data  and 
the  purpose  of  the  modeling  effort.  It  is  essential  to  distinguish 

understanding  of  the  processes  being  modeled. 


PURPOSE/OVERVIEW 


The  purpose  of  this  work  was  to  investigate.  I 

sorption.  An  additional  objective  was  to  develop  conceptual  and 
mathematical  models  for  the  transport  of  solutes  influenced  by 

nonequilibrium  are  investigated  in  Chapter  4,  where  the  methods 

evaluated  in  Chapter  5.  Solute  transport  at  the  field  scale  is 
investigated  in  Chapter  6;  a model  that  accounts  for  multi— source  and 
scale-dependent  nonideality  is  developed. 


EXPERIMENTAL  APPROACHES 


CHAPTER  3 

FOR  INVESTIGATING  SORPTION  NONEQUILIBRIUM 


Introduction 


K.  sorption  kinetics  of  organic  chemicals  has  recently  become  an 
area  of  research.  In  the  Investigation  and  modeling  of  sorption 


(HOC)  techniques.  The  GP  technique  vas  first  applied  to  sorption  by 
uclear  aromatics  from  freshwater  sediments.  The  GP  technique  has 


since  been  used  by  ochers  (Karickhoff  and  Morris,  1983;  Oliver,  1983; 
Coates  and  Elzerman,  1986).  To  employ  the  GP  technique,  a reactor 

purged  with  a gas  such  as  air  or  nitrogen.  The  gas  scream  scrips  the 
solute  from  the  water  phase,  thus  inducing  desorption  of  Che  sorbace 
from  the  sorbent.  The  gas  stream  is  trapped  (e.g.,  cenax)  upon  exiting 

values  for  kinetic  parameters.  This  technique  as  described  above 
provides  data  only  for  the  desorption  phase.  A modified  approach  was 
utilized  by  Wu  and  Gschwend  (1986)  to  investigate  both  adsorption  and 
desorption  phases.  The  apparatus  was  designed  to  operate  under  closed 

nondestructive  detector  is  required.  A photolonitation  detector  (PID) 

Miscible  displacement  techniques  have  been  widely  used  in  chemical 
engineering,  petroleum  engineering,  and  soil  science  to  Investigate 

of  Che  MDC  technique  to  investigate  the  nonequilibrium  sorption  or 
sorption  kinetics  of  organic  solutes  was  initiated  by  soil  scientists 
(e.g.,  Elrick  at  al.,  1966;  Kay  and  Elrlck,  1967;  Davidson  at  al.,  1968; 

were  evaluating  the  behavior  of  pesticides.  The  MDC  technique  has  since 
been  used  to  investigate  the  sorption  kinetics  of  ocher  types  of  organic 
chemicals,  such  as  chlorinated  aromatics  and  ethanes  (e.g., 

Schwarzenbach  and  Westall,  1981;  Lee  et  al.,  1988). 


There  has  been  to  dace  no  critical  evaluation  either  of  Che  general 
applicability  or  Che  equivalency  of  these  two  techniques.  This  paper 
reports  Che  results  of  experiments  that  vere  conducted  to  elucidate  the 
conditions  under  vhich  the  two  techniques  are  viable.  Various 
modifications  that  enhance  viability  are  discussed.  The  performance  of 


A disadvantage  of  the  Karickhoff  apparatus 
is  that  it  can  only  be  operated  in  the  desorption  mode.  A disadvantage 
with  Che  design  used  by  Hu  and  Cschvend  is  that  desorption  can  only  be 
studied  with  the  dilution  approach  (i.e.,  a given  volume  of  soil -water- 
solute  slurry  in  a condition  of  equilibrium  is  diluted  by  adding  a given 
volume  of  "clean"  water,  which  induces  desorption).  Another 
disadvantage  is  the  detection  limit  associated  with  the  PID.  This  would 
be  of  special  concern  for  solutes  with  very  low  solubilities;  small 
changes  at  minute  concentration  levels  may  not  be  discemable.  In  such 
cases  the  design  used  by  Karickhoff  (1980) , where  a trapping  device  is 
employed,  would  be  advantageous.  The  apparatus  employed  herein  was 
designed  to  combine  the  advantageous  features  of  Che  Karickhoff  (1980) 
and  Hu  and  Gschwend  (1986)  designs  and  to  eliminate  the  associated 

The  general  design  was  based  on  that  of  Hu  and  Cschvend  (1986).  A 
one  liter  glass  reaction  vessel  was  connected  to  a stainless  steel 


portion  (approximately  30  cmJ/min)  of  the  flow  was  diverted  through  a 
parallel  loop  to  which  a flow-through  photo  Ionization  detector  (Perkin- 
Elmer  PI-52-02)  was  attached-  The  potential  for  mass  loss  from  the 
system  (i.e.,  leakage)  was  evaluated  by  Injecting  a known  mass  of  solute 
Into  the  closed  system  and  monitoring  the  output  signal  for  decay.  A 

of  every  experiment.  A second  closed  loop  bypassing  Che  reactor  was 
utilized  Co  determine  base-line  detector  response.  The  P1D  was  housed 
in  a gas  chromatograph  (Perkin-Elmer  Sigma  2000) ; Che  detector  was 
maintained  at  a temperature  of  150°C  to  prevenc  vapor  condensation.  The 
detector  output  signal,  representing  the  gas-phase  concentration  of  the 
solute,  was  recorded  by  a Linear  (model  1200)  chare  recorder. 

Modifications  to  the  system  Included  the  following: 

1)  Three-way  valves  were  placed  in-line  so  that  the  apparatus  could  be 
operated  In  either  open  or  closed  modes.  This  allowed  the  performance 
of  an  adsorption  experiment  in  the  manner  of  Hu  and  Gschwend  (1986)  and 
Chen,  with  the  seme  slurry,  the  performance  of  a desorption  experiment 
in  the  manner  of  Karickhoff  (1980) , 2)  An  in-line  metering  valve  was 
enplaced  to  control  Che  gas  flux.  This  enhanced  Che  ability  to  optimize 
experimental  conditions,  3)  An  additional  valve  was  enplaced  to  allow 
use  of  trapping  devices  when  concentrations  were  such  chat  detection  by 


was  treated  by  carbon  filtration  Cl 


operated  In  the  open  modo. 

Materials.  The  following  chemicals  were  used:  benzene, 
chlorobenzene,  1,2-dlchlorobenzene,  1,2, 4- trichlorobenzene, 

materials  were  used  in  the  experiments.  All  soils  and  aquifer  materials 
were  passed  through  a 1 mm  sieve;  the  less-than-one  mm  size  fraction  was 
utilized  for  all  experiments. 


le  results  of  ti 


reaction  vessel,  which  was  then  sealed  and  placed  on  a wrist-action 
shaker  for  48  hours  Co  allow  the  establishment  of  equilibrium.  After 
the  equilibration  period,  Che  vessel  was  attached  Co  the  experimental 
apparatus  and  the  system  was  operated  in  the  closed  mode  to  determine 
the  equilibrium  gas-phase  concentration  of  the  solute.  With  this  value, 
an  equilibrium  sorption  constant  (K,,)  could  be  determined  given  the 
following  quantities:  gas  volume  (V,),  water  volume  (V.) . sorbent  maos 
(M,),  solute  mass,  and  Henry's  constant  (H). 

Values  for  all  of  the  quantities  mentioned  were  determined  by 

for  H,  values  of  which  were  taken  from  Ashworth  et 
a steady  output  signal  was  established,  Che  system  was 
I the  open  mode  Co  effect  gas-purge  desorption. 

c each  of  the  chemicals  by  monitoring  the  response  of  the 


operated  in 


of  Che  detector  recorded  during  a desorption  experiment 

determined  with  Che  calibration  curve.  The  detector  responses  for  all 
experiments  were  within  Che  range  covered  by  the  calibration  curves. 

Bala  analviis.  Early  attempts  to  model  sorption  kinetics  employed 
a "one-sice'  model  where  all  sorption  is  race-limited;  this  approach  has 
failed  to  adequately  represent  experimental  data  (e.g.,  Schwerzenbach 
and  Wes  tall,  1981;  Rao  and  Jessup,  1983;  Uu  and  Cschwend,  1986;  Boesten 
and  van  der  Pas,  1988).  Sorption  data  have  been  found  to  exhibit  a two- 
stage  approach  to  equilibrium:  a shore  initial  phase  of  fast  uptake  or 
release  where  roughly  50»  of  total  sorption  occurs  within  minutes  to 
hours,  followed  by  an  extended  period  of  much  slower  uptake  occurring 
, Karickhoff,  1980;  Ball  and 

1985;  Coates  and  Elzerman,  1986;  Wu  and  Cschwend,  1986).  This 
behavior  is  readily  approximated  with  a bicontinuum  approach. 

experiments  were  fit  with  a biconclnuum 
: Che  two  kinetic  parameters  F and  kj. 

The  system  is  conceptualized  as  follows: 

=.<-->  C.  < --->  s,  ^ S,  (3-1) 

where  Cf  is  Che  solute  concentration  in  the  gas  phase  (M/L3),  C_  is  the 
solute  concentration  in  Che  water-phase  (M/L3),  Sj  is  the  sorbed-phase 
concentration  in  Che  equilibrium  domain  (M/M) , S,  is  Che  sorbed-phase 
10  kinetic -controlled  domain  (M/M),  and  k,  and  k,  are 


over  periods  of  days  or  me 
Roberts,  1985;  Karickhoff 


o determine  v; 


-s  (T~l) , respectively. 


3-2.) 


the  forward  and  reverse  sorption  kinetic  conscai 
S,  - F K,  C. 

dSj/dc  - k,  S,  - k,  S2;  S2  - (1-F)  Kp  C.|t  . (3-2b 

The  cocal  mass  of  solute  in  the  system  (Kj)  is  given  by 

Mt  - (V,  C.)  + (V.  C,/H)  + (M,  F It,  Cp/H)  + (M,  S2)  (3-3) 

The  mass  flux  of  solute  from  the  reactor  is  given  by 

dH,/dt  - -q  C.  - R,  (dCp/dt)  + M,  (dSj/dc)  (3.4) 

where  R,  - V,  + (yj H)  +(M,  F Kp/  H) . Equation  3-4  may  be  rewritten  in 
nondlmensional  form  and  is  a second  order  linear  ordinary  differential 
equation,  which  can  be  solved  analytically. 

The  model  required  the  data  to  be  in  terms  of  cumulative  mass  flux 
of  solute  out  of  the  system.  This  quantity  was  obtained  by  trapezoidal 
integration  of  the  concentration-time  data.  Knowledge  of  the  following 
parameters  are  required  to  run  Che  model:  Kp,  H,  gas  flux  (q),  VB,  v 
H,,  mass  of  solute,  kj,  F,  and  Ms  (fraction  of  total  solute  mass  removed 
from  reactor).  The  value  for  Kp  was  determined  as  discussed  above, 
values  for  H were  taken  from  Ashworth  et  al.  (1988);  all  other  parameter 
values,  except  for  k*,  F.  and  M,.  were  measured.  Values  for  Hf  were 
obtained  by  dividing  the  total  quantity  of  mass  removed  from  Che  reactor 

the  baseline  value  (i.e. , "zero"  concentration)  was  reached,  M, 
represents  the  fraction  of  mass  lose  during  the  experiment,  including 
the  equilibration  step.  Considering  the  nature  of  the  chemicals  and  the 
relatively  short  time,  involved,  this  loss  is  most  likely  a result  of 
volatilization  from  Che  reactor  during  the  equilibration  period  and 


he  techniques  employed  for  the  miscible 
displacement  experiments  were  similar  to  those  utilised  by  Nkedl-Kl 
et  ml.  (1987)  and  Lee  at  al.  (1988). 

Altex/Beckman  preparative  chromatography  column  (No.  252-18)  made  of 
precision-bore  borosillcate  gloss,  with  an  Internal  diameter  of  2.5  cm. 
An  adjustable  shaft  was  employed  such  that  the  length  of  the  packed  soil 
bed  could  be  varied  up  to  25  cm.  Bed  supports  on  both  ends  were  woven 
FEP  Teflon  diffusion  mesh  In  contact  with  a porous  Teflon  filter  disk  of 
pore  site  30  to  60  pm.  The  diffusion  mesh  enhanced  radial  distribution 
of  Influent  solution  and  also  minimised  dispersion  at  the  effluent  end. 
The  column  apparatus  was  designed  to  have  a minimum  void  volume  In  the 
end  plates.  Two  single-piston  HPLC  pumps  (Gilson  Medical  Electronics, 
Model  302)  were  connected  to  the  column  with  a Rheodyne  switching  valve 
(Model  7060)  to  facilitate  switching  between  solutions  with  and  without 
the  solute  of  interest.  The  system  was  designed  so  that  the  solute 
contacted  only  stainless  steel  or  Teflon. 

Effluent  from  the  column  was  directed  through  a flow-through, 
variable  wave-length  UV  detector  (Cllson  Bolochrome)  to  allow  continuous 


monitoring  of  the  solution-phase  concentration  of  the  solute,  vhich  i 
recorded  on  a strip  chart  recorder  (Linear,  Model  1200).  Previous 
miscible  displacement  studies  have  usually  employed  collection  of 
tn  analyzed  by  GC  or  HPLC,  For 

te  potential  for  loss  1 


effluent  fractions,  which  vi 

volatilization  during  the  fraction  collection  at 
use  of  che  flow-through  detection  method  eliminates  ci 
also  significantly  simplifies  the  analytical  process.  However,  this 
method  is  limited  to  a specific  sec  of  chemicals. 

To  compare  che  performance  of  the  flow-through  and  effluent 
fraction  detection  methods,  a miscible  displacement  experiment  was 
performed  wherein  an  automated  fraction  collector  (1SC0,  Hodel  ISIS)  was 
connected  in  series  to  the  uv  detector.  Samples  collected  with  the 
fraction  collector  were  analyzed  by  HPLC  (Waters  Associates,  Model  450 
OT  detector;  Waters  Associates,  WISP  710B  automated  sampler).  In  this 

detection  modes.  A nonvolatile  chemical,  atrazine,  was  selected  so  chat 
any  observed  differences  in  BTCs  obtained  with  che  two  methods  would  be 
a result  of  factors  other  than  volatilization.  The  results  of  this 
experiment  are  shown  in  Figure  3-1.  There  is  virtually  no  difference 
between  the  two  BTCs,  from  which  it  may  be  concluded  that  the  two 
detection  methods  yield  comparable  results. 

HaterUls.  The  following  chemicals  were  employed  for  the  miscible 
displacement  experiments:  triclated  water,  benzene,  chlorobenzene,  1.3- 
dichlorobenzene,  1,2, 4 -trichlorobenzene,  cecrachloroethene, 
tr ichloroe thene , and  atrazine.  Eustis  soil  was  used  for  these 
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Pore  Volumes 


seeps  with  Che  dry  material  to  establish  uniform  bul 
completion  of  packing,  electrolyte  soil 


established. 


periodically  throughout  the  experiments.  When  the  column  was  i 
it  was  carefully  disassembled  and  mass  balances  on  pore  water  i 

id  dead  volumes  were  used  to  determine  tht 


Saturated  stock  sole 
quantities  in  excess 


lectrolyte  for  all  solutions  was  CaCl2  (0.01  N> . 
is  were  prepared  for  each  solute  by  placing 
solubility  in  contact  with  the  electrolyte 


solution  and  stirring  the  solution  for  at  least  68  hours.  Solutions 
were  filter  sterilized  with  0.45  pm  filters  prior  to  use.  Aliquots  from 
the  stock  solution  were  diluted  with  electrolyte  solution  until  the 
desired  concentration  was  achieved.  The  C0  values  for  the  experimental 
solutions  were  generally  50%  of  the  aqueous  solubility. 

A typical  experiment  involved  the  continuous  injection  of  a pulse 


point  of  C/C0  - 1 was  reached.  The  solute  pulse  was  then  displaced  with 
0.01  N CaCl3  solution  until  the  point  of  C/CQ  - 0.  Constancy  of  flow 
vas  periodically  checked  by  collecting  effluent  fractions.  These 


injection 


measured  flow  races  were  checked  against  Che  average  floi 


Data  analysis.  The  results  of  the  miscible  displacement 
experiments  were  analyzed  with  a bicontinuum  model  to  attain  equivalency 

Models  of  this  type,  where  sorption  is  assumed  to  be  instantaneous  for 


nonequilibrium,  and  incraj 
nonequilibrium  sorption  (I 


tedium  and  race-limited  for  cl 
i et  al.  (1976),  van  Genuchcen  and  Uierenga 
e (1977).  These  models  can  represent  each 


ac'/ap  + (flR-i)  ac*/ap  + (i-«r  as*/ap  - u/p)  aV/ax*  - ac'/ax  (3-5) 
(1-«R  as‘/ap  - e>  (C‘  - S*>  (3-6) 


c*  - C/C0  (3-7.) 
s‘  - Sz/(l-F)Kp  (3-7c) 
p - vt/L  (3-7e) 
X-x/h  (3-7g) 


P - vL/D  (3-7b) 
R - 1+  (p/8) K„  (3-7d) 
8 - (1  + P(p/«)Kp)/R  (3-7f) 
u - k,(l-a)RVv  (3-7h) 


and  where  D is  the  dispersion  coefficient  (L*/T) , v is  the  average  pore- 
water  velocity  (L/T) , x is  distance  (L) , L is  column  length  (L) , p Is 
bulk  density  (M/Ls) , and  8 is  volumetric  soil-water  content.  Note  that 
the  kinetic  parameters  k,  (first  order  rate  constant)  and  F (fraction  of 


To  run  the  model,  knowledge  of  the  following  pi 
required:  the  Peclet  number  P,  which  represents  the  dispersive-flux 

the  effect  of  sorption  on  transport;  the  fraction  of  instantaneous 
retardation,  t> ; the  Damkohler  number  w,  which  is  a ratio  of  hydrodynamii 
residence  time  to  characteristic  time  of  the  sorption  "reaction";  and 

normalized,  absolute  temporal  moment  of  a dirac  pulse.  The  general 


"i  - of  C*P  «?/  oT  C-  dp  (3-8) 

Application  of  eq  3-8  to  a BTC  will  not  directly  yield  a value  for  R for 
a non-Dirac  Input  pulse,  i.e.,  when  T0  is  finite.  A correction  factor 
of  0.5T0  (Valocchi,  1985)  was  applied  to  quantities  determined  with  eq 
3-8  to  calculate  the  R value.  The  value  of  T„  vas  known  from 


-squares  optimization  program  CF1TIM  (van 
o determine  values  for  Che  two  unknowns. 


Remits,  ana  Pis; 


The  results  of  gas-purge  di 

muck  soil  and  four  solutes , benzene  (BNZ) , chlorobenzene  (CB; 

.,2,4-trichlorobenzene  (TCB) , are 

The  relationship  between  k. . values  of  which  > 

energy-relationship  technique,  in  the  manner  of  that  to  be  discussed  in 
Chapter  4 of  this  dissertation.  The  results  are  shown  in  Figure  3-3. 

relationship  has  been  observed  by  others  (Karickhoff,  1980;  Karickhoff 
and  Morris,  1985)  and,  for  the  types  of  systems  Investigated  herein,  is 
suggestive  of  an  intrasorbent-diffusion  process  (Karickhoff  and  Morris, 

chlorobenzenes  by  Karickhoff  and  Morris  (1985),  Oliver  (198S),  and  Wu 
and  Gschwend  (1986).  To  compare  the  results  obtained  by  these 
researchers  to  those  obtained  herein,  values  for  k?  and  Kp  obtained  from 
Che  research  reports  are  plotted  in  Figure  3-3.  A combined  regression 
of  the  literature  and  experimental  data  yields  an  r2  of  0.99.  Kith  such 

herein  are  comparable  to  those  obtained  by  Che  ocher  researchers.  The 
GP  technique  would,  therefore,  seem  to  hove  associated  with  it  a high 
degree  of  reproducibility  and  consistency. 
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Time  (min) 


i tine  than  TCE. 


ir  desorption  of  tetrachloroethene 


Figure  3-2, 
requiring  n< 

optimization  were  54  hr-1  and  0.30  for  TCE  and  6 hr*1  and  0. 

TCE  and  PCE  through  a Eustis- 
•e  performed  with  TCE  and  six  ocher  soils 
and  aquifer  materials  in  addition  Co  Euscls  soil.  The  values 
determined  from  these  experiments  are  compared  in  Figure  3-4  to  the 
regression  equation  reported  in  chapter  four  of  this  dissertation  for 

The  two  question  marks  in  Figure  3-4  represent  the  indeterminant 
nature  of  the  results  obtained  for  the  two  aquifer  materials.  The 
characteristic  reaction  time  of  TCE  sorption  to  these  two  sorbents  was 
fast  enough  such  that  Che  desorption  curve  could  not  be  distinguished 

parameters  determined  from  a CP  experiment,  there  must  be  a significant 


%6o| 


o 


greater  than  the  characteristic  tine  for  gas-liquid  mass  transfer 
Given  the  inverse  relationship  that  appears  to  exist  between  kj  at 


the  condition  under  which  Td  » Tfl  and 
technique  would  be  viable,  is  strongly 
the  specific  solute/sorbent  pair.  Base 
Figures  3—3  and  3-4,  it  appears  that  tfc 
systems  having  K,  values  ranging  from  o 
becomes  unreliable  for  Kp  values  in  the 


the  sorptivity  of 


A tritiated  water  BTC  was  obtained  for  t 
hydrodynamic  character  and  to  determine  a val- 
5).  This  BTC  was  symmetrical  and  sigmoidal  i: 
with  pore-water  velocity,  vl 
nonequilibrium  (Brusseau  and 
dlchlorobenrene  (DCB3) , and 
asymmetry  of  the  BTCs;  this. 

nonequilibrium  process  (Brus 
Values  for  kj  and  Kp  obi 
Figure  3-6,  along  with 

between  kj  and  Kp  i 


o establish  the 

h suggests  the  absence  of  physical 

B are  presented  in  Figure  3-5.  Note  the 
n conjunction  with  the  symmetrical  BTCs 
Is  highly  suggestive  of  a sorption-related 

'e  plotted  in 


s-purge  experiments.  The  similarity  between  the  ti 
hat  the  two  experimental  techniques  produced 
It  is  also  evident  that  an  inverse  relationship 
displayed  by  the  MDC  dace. 
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attempting  to  use  kinetic-parameter  values  determined  by  gas-purge  to 
predict  a BTC  obtained  for  the  same  soil/solute  pair.  The  BTC  obtained 


:he  Eusels  soil  column  is 
produced  by  Che  biconcinuum 


from  miscible  displacement  of  P< 
presented  in  Figure  3-7.  The  Vi 

o (see  eq.  3-7f  and  3-7h) . The 
obtained  from  miscible  displacement  of  TCE  through  the  Eustls  column 


(0.25).  which  is  within  Che  range  where  Che  GP  technique  may  be 
unreliable.  These  two  predictions  to  my  knowledge  represent,  along  vlth 
that  reported  in  Chapter  5 of  this  dissertation,  the  first  successful 
cases  where  nonequilibrium  BTCs  of  organic  solutes  were  predicted  with 
all  parameter  values  being  obtained  independently  (i.e.,  without  curve 
fitting) . 


yield  comparable  results,  based  on  tl 
3-7  and  3-8.  This  would  suggest  thar 
same  physico-chemical  process.  This 
model,  since  it  was  used  successfully  in  data  analysis  for  boch  methods, 
provides  a valid  representation  of  the  process  responsible  for  the 
observed  nonequilibrium.  whatever  it  may  be.  The  parameter  kj.  rather 


nd  miscible  displacement  techniques 
also  suggests  that  the  biconcinuum 


o 
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chan  being  jusc  a fitting  function,  may  thus  be  considered  as  being 
represencacive  of  Che  nonequilibrium  mechanism. 

useful  for  low-sorptivity  systems.  Its  efficacy,  however,  is  greatly 
reduced  for  systems  comprised  of  highly  sorptive  chemicals  or  of  soils 
containing  high  levels  of  clay  and/or  organic  matter.  As  the  sorpeivity 
of  the  solute  Increases,  boch  time  constraints  and  problems  such  as 
sorption  to  Che  apparatus  become  of  increasing  concern.  These  problems 
can  be  ameliorated  through  modification  of  the  experimental  method.  One 
such  approach  Involves  Che  use  of  miscible  organic  cosolvents . The 
presence  of  cosolvenc  results  in  reduced  values  of  Kp  (Rao  ec  al.,  1985; 
Nkedi-Kizza  ec  al.,  1985;  Fu  and  Luthy,  1986b;  Nkedi-Kizza  et  al.,  1987) 
and  in  Increased  values  of  k,  (Nkedi-Kizza  et  al.,  1989).  Bith  these 
effects,  the  use  of  cosolvent  can  significantly  reduce  the  experimental 


highly  sorptive  chemicals.  It  will  be  shown  in  Chapter  4 of  this 
dissertation  that  values  for  kj  in  aqueous  systems  can  be  estimated 
successfully  by  extrapolation  from  experiments  performed  with  mixed- 

The  determination  of  kinetic  parameter  values  from  column 
experiments  is  predicated  upon  the  ability  of  the  mathematical  model  to 
successfully  simulate  the  experimental  data.  Confidence  in  the 
robustness  of  the  parameter  values  so  determined  is  attained  only  with  a 
unique  solution  (i.e.,  when  one  suite  of  parameter  values  provides  a 
solution  chat  is  significantly  better  than  all  others).  For  cases 
wherein  a system  is  near  equilibrium  or  under  extreme  nonequilibrium. 


attainment  of  a unique  solution  may  prove  difficult.  A modified 
miscible-displacement  technique,  involving  flow  interruption,  that 
potential  for  achieving  unique  solutions,  and  thus 

of  optimized  values  of  kinetic  parameters,  ui] 
chapter.  This  method  would 


be  presented  in  the  next  section  of 
appear  to  be  especially  useful  for  systems  comprised  of  aquifer 
materials , as  a low  degree  of  nonequilibrium  may  often  be  expect 
such  systems  based  on  their  relatively  low  sorptivity. 


Two  experimental  methods  for  investigating  nonequilibrium  sorption 
of  organic  chemicals  have  been  evaluated.  In  planning  to  use  either 
method  it  is  important  to  consider  their  range  of  applicability.  The 
viability  of  the  gas-purge  technique  is  controlled  by  the  vapor  pressure 
of  the  chemical  and  the  sorptivity  of  the  solute/sorbent  combination. 

The  method  will  obviously  be  applicable  only  for  chemicals  whose  vapor 
pressure  is  such  that  gas-phase  concentrations  are  of  sufficient 
magnitude  to  be  detected.  The  viable  rang,  of  vapor  pressure  can  be 
extended  through  the  use  of  trapping  devices.  The  viability  of  the 
technique  becomes  questionable  when  K,  values  are  less  than  0.5  ml/g. 
Below  this  Kp  value,  the  characteristic  time  of  sorption  is 
Indistinguishable  from  that  of  gas-liquid  m,,  transfer.  The  visbility 
of  the  miscible  displacement  technique  is  controlled  by  the  sorptivity 
of  the  solute/sorbent  combination  and  by  the  texture  of  the  sorbent. 

This  technique  is  especially  useful  for  systems  comprised  of  chemical, 
of  low  to  moderate  sorptivity  and  of  sandy  soils  or  aquifer  materials. 


The  viable  range  may  be  extended  wich  modifications  such  as  mixed- 
solvent  solutions  or  flow  interruption.  Considering  the  results 
discussed  above,  the  gas-purge  technique  would  appear  to  be  the  method 


especially  for  those  wich  soils  containing  high  clay  or  organic  matter 
contents.  The  miscible  displacement  technique  would  be  Che  method  of 
choice  for  systems  of  low  and  moderate  sorpeivity. 

The  time  required  for  desorption  was  observed  co  be  a function  of 
the  solute;  as  the  sorpeivity  increased,  Che  required  time  also 
Increased.  An  inverse  relationship  between  k2  and  Kp  was  observed  for 
the  results  obealned  from  both  experimental  techniques.  This,  in 
combination  wich  the  observation  of  asymmetrical  BTCa  for  the  sorbing 
solutes  and  symmetrical  BTCs  for  the  nonsorbing  solute,  suggests  that 
Che  observed  nonequilibrium  was  a result  of  on  Intrasorbent-dlffuslon 


EIb"  IntemiBtion; A HetW  for  InvescUfltina.gQrptton 


investigate 


re  is  a definite  need  to 


practice  ofeen  Involves  curve-fitting  experimental  data  with  various 
models  and  selecting  the  model  that  exhibits  Che  best  fit.  However,  the 
successful  simulation  of  a data  sec  by  a mathematical  model  is  a 
necessary,  but  not  sufficient,  condition  for  establishing  the  validity 


conceptual  model  are 

nonequilibrlum  for  the  case  of  linear  sorption  isotherm  (van  Genuchten, 
Rao,  1989;  Lee  et  al.,  1988).  Obviously,  there  is  need  for  experimental 

cannot  be  readily  discerned  when  analyzing  data  obtained  with  the 
typical  column  experiment.  In  these  instances  a technique  chat  has  a 
greater  sensitivity  to  nonequilibrlum  would  be  very  advantageous.  When 

; are  required  for  the  kinetic  terms.  Under  certain 
current  curve-fitting  approach  may  not  produce  a unique 
tt  of  values  for  a given  suite  of  parameters.  To  Improve  the  efficacy 
f parameter  determination,  this  ability  must  be  enhanced, 

nsicivity,  it  also  provides  more  reliable  determinations  of  kinetic 


possible 


technique  to  discriminate  between  the  mechanism(s)  responsible  for 
sorption  nonequilibrium.  The  conceptual  and  theoretical  basis  for  the 
method  is  presented,  and  its  application  is  illustrated. 

Conceptual  Basis  for  the  Method 

and  Myers,  1947;  Kressman  and  Kitchener.  1949)  which  has  been  employed 
in  chemical  engineering  to  distinguish  between  systems  chat  are 

Helfferlch  (1962) , this  is  the  best  technique  for  differentiating 

was  just  prior  to  interruption.  If,  however,  film  diffusion  is  rate 
controlling,  the  interruption  will  have  no  observable  effect  upon  the 

Weber  and  Morris  (1963)  employed  the  interruption  method  to 

carbon.  Zogorsky  ec  al.  (1976)  employed  Che  method  Co  investigate 
sorption  kinetics  of  2 . 4-dichlorophenol  on  activated  carbon  and 


that  noe  only  can  differentiate  between  filn  and  intra-particle 
diffusion,  but  also  intra-organlc  matter  diffusion  and  rate-limited 

in  a batch  system;  hence  the  name—  the  flow  interruption  method.  The 
method,  then,  involves  displacing  at  a steady  flow  rate  a solution 
containing  the  solute  of  interest  into  a water-saturated  soil  column 
stopping  the  flow  for  some  period  of  time,  and  then  reinitiating  flow. 
Either  temporal  changes  in  the  effluent  concentrations  (i.e., 
breakthrough  curves  (BTC))  or  concentration  profiles  within  the  column 
may  be  monitored.  The  former  is  experimentally  moro  feasible  and  is 
also  theoretically  more  appealing  since  the  temporal  first  moment  (i.e. , 
center  of  mass,  equivalent  to  retardation  factor)  is  independent  of 
kinetic  constraints,  whereas  the  corresponding  spatial  moment  is  not 

The  flow-interruption  technique  is  similar  in  operation  to  the  so- 
called  "arrested-flow"  chromatography  technique,  first  described  by  Knox 
and  Mclaren  (1964).  This  technique,  however,  was  developed  to  measure 
gas-diffusion  coefficients  and  not  to  investigate  sorption  kinetics. 

The  technique  is  employed  by  stopping  the  flow  of  carrier  gas  after  the 
pulse  has  reached  the  midregion  of  the  column,  and  then  recommencing 
flow  after  some  period  of  time.  During  the  interruption,  solute  spreads 
only  by  molecular  diffusion.  In  the  absence  of  sorption  or  chemical 
reaction,  the  diffusive  and  dispersive  processes  are  linear,  and 
The  Influence  of  hydrodynamic  dispersion  and 


therefore  additive. 


persion. 


entrance/exit  boundary  effects  can  be  subtracted 
concribucion  of  molecular  diffusion  Co  solute  di 
or  second  moment  of  the  BTC  (l.e.,  dispersion)  varies  with  the  length  of 
interruption.  Hence  the  slope  of  a plot  of  variance  versus  interruption 
tine  Is  proportional  to  the  molecular  diffusion  coefficient.  The  method 
was  recently  revived  by  McCoy  and  Moffat  (1986)  and  Park  ec.  ml.  (1987). 


flow  interruption  by  defining  a modified  pore 


be  modified  to  simulate 
volume  parameter: 


following  sorption  nonequilibrium  mechanisms—  physical  nonequilibrium, 
chemical  nonequilibrium,  or  diffuslonal  mass  transfer  of  sorbace  within 
mineral  or  organic  components  of  the  sorbent. 

boundary  condition  and  an  initial  condition  of  C*-0  using  a Crank- 
Ntcolson  finite  difference  method.  Parameter  estimation  was 
accomplished  by  coupling  the  numerical  model  with  a non-linear 
least-squares  optimization  program  (NUS)  (Meeter  and  Wolfe,  1968). 


When  sorption  nonequilibrium  exists,  flow  Interruption  on  the 
frontal  limb  of  a BTC  will  result  in  a drop  in  the  effluent  solute 
concentration  (except  when  film  diffusion  controls  the  race  of 
sorption).  An  example  of  a BTC  obtained  from  a simulated  flow- 


interruption  experiment 


Figure  3-9  correspond  Co  the  pre-  and  posc-interrupc  effluent 
concentrations,  respectively.  The  dip  in  the  BTC  is  a result  of 

thus  reducing  the  in  situ  aqueous-phase  concentration.  The  solute 


is,  therefore,  lower  chan  Che  concentration  Just  prior  to  interruption. 

breakthrough  front  of  the  solute  pulse  Injected  upon  recommencing  flow. 
To  maximize  Che  magnitude  of  the  resultant  drop  in  effluent 

le  interruption  should  occur  within  a specific  range  of 

lemlcal  as  well  as  Che  sorptive  nacure  of  the  sorbent. 
Longitudinal  diffusion  of  the  solute  in  the  resident  solution  can  occur 
during  the  quiescent  period  as  a result  of  concentration  gradients  along 

diffusion  during  conditions  of  no  flow  was  demonstrated  by  Bigger  and 
Nielsen  (1962b)  and  Fukui  and  Kacsurayama  (1975) . Molecular  diffusion 
will  redistribute  the  solute  from  the  region  near  Che  influent  end  (X-0) 
toward  the  region  near  Che  effluent  end  (X-l),  thus  increasing  the 
post-interruption  effluent  concentration.  Because  such  solute 
redistribution  opposes  the  effects  of  race-limited  sorption,  the  true 
magnitude  of  the  concentration  drop  can  be  masked.  The  lower  the  C* 


Pore  Volumes 


greater  the  resultant  solute  redistribution.  Thus,  to  minimize  the 
effects  of  longitudinal  molecular  diffusion,  the  flow  Interruption 

solution  concentration  gradient  is  small. 

A flow  Interruption  experiment  with  tritleted  water  was  performed 
to  demonstrate  the  effect  of  longitudinal  diffusion  during  the 
interruption  period.  Displacement  of  tritiated  water  through  a column 
packed  with  Eustis  fine  sand  (a  sandy  soil  with  organic  carbon  concent 
of  0.56%)  was  stopped  at  the  Initial  stage  of  solute  breakthrough 
(C‘-0.02).  The  interruption  period  was  approximately  21  hours.  Upon 
recommencing  flow,  the  Initial  value  of  C*  was  0.25.  For  a BTC  measured 
under  similar  experimental  conditions  but  with  no  flow  Interruption,  the 
relative  effluent  concentration  at  the  equlvelent  pore  volume  was 
approximately  0.08.  Solute  redistribution  resulting  from  longitudinal 
diffusion,  then,  caused  C‘  Co  change  by  approximately  0.17. 

The  effect  of  interrupting  flow  at  a low  C*  value  for  a sorbing 
solute  is  demonstrated  in  Figure  3-10;  note  that  all  parameter  values 
are  identical  to  chose  for  Figure  3-9,  except  for  the  time  at  which  the 
interruption  was  initiated.  The  magnitude  of  the  concentration  drop  is 
negligible,  compared  Co  that  observed  in  Figure  3-9.  At  lower  c', 
decreeses  in  solution  concentration  resulting  from  continued  sorption 
are  negated  by  an  Increase  in  concentration  resulting  from  the 
occurrence  of  longitudinal  diffusion  during  the  quiescent  period.  The 
magnitude  of  the  drop  obcalned  at  low  C'  is  such  Chat  it  could  not 
readily  be  discerned  from  eh.  scatter  typical  to  experimental  data. 


Pore  Volumes 


Interrupting  flow  at  a higher  C‘  value,  however,  nay  not  allow  the 
development  of  a significant  concentration  drop  when  near-equilibrium 
conditions  exist.  For  such  systems,  the  interruption  muse  be  initiated 
at  somewhat  lower  C*  values  to  maximize  conditions  for  elicitation  of 
the  concentration  drop.  Hence,  the  opposing  effects  of  the  degree  of 
nonequilibrium  and  of  longitudinal  diffusion  during  the  quiescent  period 
must  be  considered  in  determining  the  optimal  C"  value  for  initiating 
flow  interruption. 

It  should  be  noted  that  the  interruption  can  also  be  initiated  on 
the  distal  portion  of  the  BTC.  In  this  case,  an  increase,  rather  than  a 
drop,  will  occur  in  the  effluent  concentration  when  nonequilibrium  is 
present.  To  optimize  the  magnitude  of  this  spike,  the  interruption 
should  occur  at  lower  c'  values.  In  addition  to  the  single- 
perturbation,  stop-flow  approach  discussed  so  far,  flow  variation 
(period  of  decreased  or  Increased  flow  velocity)  and  multi-point 
perturbation  (multiple  periods  of  flow  interruption  or  variation) 
approaches  may  also  be  employed. 

One  problem  associated  with  BTCs  (non  Interrupted)  is  the 
relatively  poor  resolution  that  exists  in  the  regions  near  C*  - 0 and  1. 
As  the  C‘  value  approaches  the  limiting  value  of  one  (or  zero) , the 
changes  occurring  in  C‘  become  very  small.  When  the  magnitude  of  the 

:o  eh«  analytical  sensitivity,  resolution  is  lost.  In 
it  at  which  C*  actually  reaches  one  (or  zero)  is 
a result  of  this  poor  resolution,  BTCs  that  appear  to 
e.  , to  have  reached  C*  - 0 or  1)  may  actually  have 
.n  example  of  a BTC  that  appears  to  have  reached  c‘-l 


tailing  present. 


sensitivity.  The  lncerrupeion  period  allows  the  effects  of  sorption 
nonequilibrium  Co  become  more  apparent,  thus  enhancing  che  ability  Co 
delineate  and  investigate  conditions  of  sorption  nonequilibrium. 

Problems  may  exist  when  NILS  optimization  methods  are  used  to 
estimate  che  model  parameters  for  non-interrupced  BTCs  Influenced  by 
nonequilibrium.  This  is  a result  of  two  factors:  (1)  che  curve  fitting 
is  weighted  cowards  che  central  portion  of  the  BTC  (e.g.,  C*  - 0.2  Co 


0.8)  because  che  sum  of  the  squares  of  che  deviations  to  be  minimized  li 
taken  in  Che  vertical  direction.  This  portion  of  che  BTC  is  not  as 
sensitive  Co  nonequilibrium  effects  as  che  lower  (early  breakthrough) 
and  upper  (tailing)  portions,  which  are  given  less  weight;  (2)  The 
portions  of  che  BTC  most  sensitive  to  nonequilibrium  are  also  those 
portions  most  difficult  to  determine  experimentally.  As  a result  of 
these  two  factors,  curve-fitting  of  non-interrupced  BTCs  can,  in  some 
cases,  produce  ambiguous  parameter-value  determinations.  This  is 
demonstrated  in  Figure  3-12,  where  it  is  seen  that  different 


same  parameter  values  used  to  produce  Figure  3-12  were  also  used  to 
simulate  the  flow-interrupted  BTCs  shown  in  Figure  3-13.  The  presence 
of  the  concentration  drop  in  BTCs  obtained  with  che  flow-interruption 
method  supplies  additional  ■information"  on  sorption  nonequilibrium, 
thus  enabling  more  reliable  determinations  of  kinetic  parameter  values. 
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Race— limited  sorption  reactions  are  not  the  only  processes  that 
can  cause  a concentration  drop.  Transformation  reactions  (biological  or 
chemical)  occurring  during  Che  interruption  period  would  also  result  in 
a reduction  in  in  situ  solute  concentrations.  These  reductions  would 
increase  Che  magnitude  of  Che  observed  drop  in  effluent  concentration 


Process  investigation 

exhibits  a concentration  drop,  physical  nonequlllbrlum  (i.e.,  diffusive 
mass  transfer  between  "mobile"  and  "immobile"  regions)  is  presumed  to  be 

(2)  Perform  flow  interruption  with  Che  sorbing  solute  of  interest. 

physical  nonequilibrium  is  present,  the  kinetic  parameters  determined 
for  the  nonsorbing  solute  may  be  used  to  simulate  the  BTC  of  the  sorbing 


predicted  well,  physical  nonequilibrium 


responsible  for  sorption  nonequilibrium, 
race-limited  sorbace-organic  matter  il 

interaction  between  Che  sorbate  and  the  mineral  components  is  suggested 

(4)  Perform  a column  experiment  without  flow  interruption  to 


r partially  <or  totally) 
s reduced  (or  eliminated) , 


d apparently  be  cl 


e significance  < 


insignificant.  This 
diffusion 


ssenc  the  contributions  of  all 

it  is  good,  film  diffusion  is 

generally  insignificant  relative  to  intra-particle 
r many  systems  of  interest  (Leenheer  and  Ahlrichs,  1971; 


; Crittenden  el 


>;  Roberts  et  si.,  1987). 


Application  of 


with  unpublished  experimental  data  of  Rao  and  Davidson  (1979) , who 
obtained  BTCs  from  miscible  displacement  of  tritiated  water  as  a 


techniques  w< 


•-D  has  a solubility  of  890  mg/1  2S°C  (Verschueren.  19B3) , a log  K„ 
2.81  20°C  (Verschueren,  1983),  and  a pK.  of  2.64  (Worthing,  1983). 

:s  for  these  solutes  were  obtained  for  displacement  through 
;er-saturated  columns  of  Eustls  and  Webster  soils.  Analytical 
lermination  of  3H  and  ;'c  activity  in  the  effluent  soludon  fractions 
* accomplished  using  liquid  scintillation  counting.  Experimental 
similar  to  those  employed  by  Rao  et  al.  (1979). 
m displacement  of  tritiated  water  through  the  Eustis 
in  Figure  3-14.  Flow  was  interrupted  at  c"  - 0.75  such 
n changes  resulting  from  molecular  diffusion  were 
likely  to  be  insignificant.  Rote  Che  absence  of  a drop  in  the  BTC 
which  suggests  that  physical  nonequilibrium  was  not  present.  Absence  of 
mobile-immobile  soil-water  regions  can  be  expected  for  this  system 
(sandy  soil,  grain  size  less  chan  1 mm,  saturated  conditions).  The  BTC 
for  displacement  of  2,4-D  amine  through  the  Eustls  column,  however,  does 
exhibit  a significant  drop  (Figure  3-15).  The  combinations  of  Figures 
3-14  and  3-15  suggests  the  existence  of  a sorption-related 


nonequilibrium  mechanism.  A BTC  for  2,4-D  amine  was  also  obtained 
(Figure  3-16)  for  a column  packed  with  Eustls  soil  that  had  been  treated 
with  hydrogen  peroxide  (H,0;)  to  partially  remove  organic  matter 


(organic 


).1»). 
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Comparing  Figures  3-15  and  3-16,  a substantial  decrease  in  Che  magnitude 
of  the  effluent  concencracion  drop  was  observed  for  che  experiment 
involving  HaO,-created  Eustis  soil.  From  this,  it  may  be  inferred  that 

organic  matter.  Considering  che  porous,  polymeric  structure  attributed 
to  soil  organic  matter  (e.g. , Khan,  1978;  Burchill  et  al.,  1981),  the 
nonequilibrium  mechanism  may  Involve  intra-organic  matter  diffusion. 

This  hypothesis  was  also  proposed  by  Khan  (1973),  who  attributed  the 
kinetic  sorptive  behavior  exhibited  by  2,4-D  to  the  diffusion  of  sorbate 


two-site  sorption  model,  using  values 

a chat  the  fit  was  obtained  with  fixed, 
calculate  R„  were  determined  from  batch 


for  the  kinetic 
in  Figures  3-15 

experiment.  The  fit  obtained  with  the  R„  value  optimized  was 
essentially  the  same  as  that  obtained  with  R„  fixed.  This  supports  the 
validity  of  using  independently-determined  K,  values  for  simulating 
column  BTC  data;  the  same  conclusion  has  been  reached  by,  among  others, 
Lee  et  al.  (1988). 

The  BTC  for  displacement  of  2,4-D  amine  through  a water-saturated 
column  of  Webster  soil  is  given  in  Figure  3-17,  along  with  the  solaced 


(untreated) — 0 . li 


larger  chan  chat 


ss  determined  from  che  model  simulations  were;  Eustis 
!;  Eustis  (created)— 0.44;  Webster— 0.03.  Note  chat  che 
! experiment  employing  Eustis  soil  created  with  H;02  was 
obtained  from  the  experiment  using  the  non-treated 
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two  soils,  the  Webster  value  is  small* 
transfer  employed  in  the  IOHD  model. 


5 the  kz  values  obtained  for  the 

(OC)  (Webster  0C-  3.9%;  Eustis 
■e  conceptualization  of  mass 
Ihe  mass  transfer  coefficient  (k2) 

diffusive  path  length  is,  in  turn,  correlated  to  the  quantity  of  organic 
matter  associated  with  a sorbent.  Thus,  a decrease  in  kj  is  expected 

also  been  observed  by  Karickhoff  (1980) , Karickhof  f and  Morris  (1985) , 
Walters  and  Guiseppie-Elie  (1988),  and  Brusseau  and  Rao  (1989). 

and  Aylmore,  1980);  however,  1 am  not  aware  of  any  previous  application 
demonstrated  with  miscible  displacement  of  2,4-D  through  saturated  soil 

suggested  to  be  a result  of  incra-organlc  matter  diffusion. 
Experimentally-determined  k2  values  show  an  Inverse  relationship  to 
organic  matter  concent,  as  would  be  expected. 


The  flow-interruption  technique  should  prove  useful  in  identifying 
Che  presence  of  nonequilibrium  and  in  evaluating  the  conditions  under 
which  it  is  important.  The  method  should  also  provide  improved 

responsible  for  sorption  nonequilibrium.  Although  the  method  is 
indirect,  its  significance  is  not  lessened,  considering  the  present  lack 
of  direct,  molecular-scale  techniques. 

In  addition  to  its  use  as  a process-investigation  technique,  there 

restoration  of  contaminated  aquifers  by  flushing  ( i . e . , contaminant 

such  is  Che  case,  the  concentration  of  solute  in  the  groundwater  will, 
after  an  initial  rapid  decline,  exhibit  a long,  asymptotic  approach  to 

efficiency  of  flushing,  pulsed-flushlng,  rather  than  continuous  pumping, 
times  as  necessary.  In  a similar  vein,  the  pulsed-flow  approach  may 


applied  eo  soil-venting,  Che  relatively  new  technique  developed 
re  volatile  chemicals  from  the  unsacuraced  zone. 


CHAPTER  4 

INVESTIGATION  OF  SORPTION  NONEQUILIBRIUM 


Sorption  Kinetics 


Introduction 

As  discussed  earlier,  sorption  nonequilibrium  may  result  from 
chemical  nonequilibrium  (ONE)  or  Intrasorbent  diffusion.  Given  the 
nature  of  HOCs,  chemical  nonequlllbrium  is  not  likely  to  be  an  Important 
process  (Chiou  et  al..  1979;  Karlckhoff,  1981;  Brusseau  and  Rao,  1989). 
Thus,  intrasorbent  diffusion  may  be  the  predominant  nonequilibrium 
process  for  HOCs.  Intrasorbent  diffusion  may  be  separated  Into  two 
classes,  lntraorganlc  matter  diffusion  (IOHD)  and  diffusion  into  mineral 
particles  (1MD) , according  to  the  two  general  types  of  sorbents. 


mineral  particles  is  not  likely  to  be  important  under  most  conditions. 
This  is  supported  by  the  work  of  Lafleur  (LaFleur  et  al.,  1978;  LaFleur 
1979a;  1979b),  who  reported  that  sorption  for  a large  number  of 
pesticides  on  inorganic  substrates  was  independent  of  time.  This 
indicated  that  sorption  seemed  to  be  restricted  to  external  surface 
sites,  and  that  sorption  on  internal  surfaces,  if  present,  was 
negligible.  Conversely,  sorption  to  organic  matter-containing  sorbents 
exhibited  a significant  long-term  approach  to  equilibrium.  Similar 


behavior  has  been  reported  by  Siegel  ec  el.  (1951),  Harris  and  Barren 
(1964),  Moyer  et  al.  (1972),  Salesman  ec  al.  (1972),  and  Grover  (1974). 

Incramineral  diffusion  may,  however,  be  important  for  ocher 
chemical  classes.  Por  example,  the  long-term  sorption-desorption 
behavior  exhibited  by  Fluridone  on  chree  California  soils  was  attributed 
to  diffusive  transfer  to  incerlamellar  regions  of  the  clay  fraction 

Diffusion  within  the  matrix  of  sorbent  organic  matter  has  been 
proposed  as  a possible  mechanism  responsible  for  race-limited 
sorption/desorpeion  by  several  authors  (Hamaker  et  al.,  1966;  Hayes, 

1970;  Leenheer  and  Ahlrichs,  1971;  Bums  ec  al.,  1973*  Khan  1973- 
LaFleur  et  al.,  1978;  Moreale  and  van  Bladal.  1979;  Kenaga  and  Goring, 
1980;  Rao  and  Jessup,  1983;  Freeman  and  Cheung,  1981;  Boescen  and  Van 
Der  Pas,  1983;  Isaackson  and  Frink,  1984;  Karickhoff  and  Morris,  1985; 
Oliver.  198S;  Bouchard  et  al. . 1988;  Leo  et  al.,  1988;  Brusseau  and  Rao, 


Partitioning  into  soil  organic  matter  has  been  proposed  as  the 
primary  sorption  mechanism  for  HOCs  (Chiou  et  al.,  1979;  Moreale  and  v 
Bladel,  1979;  Freeman  and  Cheung,  1981;  Wershaw,  1986).  Soil  organic 


matter  may  be  pictured  as  a three-dimensional  network  of  randomly 
oriented  polymer  chains  and  as  having  a relatively  open,  flexible 
structure  perforated  with  voids  (Degens  and  Hopper,  1976;  Khan,  1978; 
Schnlczer,  1978;  Burchlll,  1981).  Humic  materials  in  soils  are  pictured 
as  membrane-like  aggregates  in  a model  proposed  by  Wershaw  (1986).  With 
such  a "porous"  structure  attributed  to  organic  matter,  the  existence  of 
an  intraorganic  matter-diffusion  process  is  definitely  possible.  If 


sorption  nonequilibrium  Is  caused  1 
is  investigated  in  this  section. 


■ IOMD,  the  structural  n 


Materials  und. Methods 

The  literature  was  searched  for  data  on  the  nonequilibrium 

materials.  Values  for  the  partition  coefficient  (Kp)  and  the  sorption 
rate  coefficient  (k2)  were  desired.  To  facilitate  data  analysis,  the 
linear  free-energy  relationship  (LFER)  technique  was  employed.  This 
correlation  analysis  technique  is  widely  used  in  chemistry  to  evaluate 

linear  relationship  occurs  when  the  logarithms  of  t 
a series  of  compounds  are  plotted  against  the  logarithms  of  their 

All  of  the  data  were  obtained  with  one  of  three  experimental 

batch-rate  and  the  gas-purge  techniques  had  been  analyzed  with  the  use 
of  a blcontinuum  model;  in  these  cases,  the  appropriate  values  of  Kp  and 

(1985)  presented  batch-rate  data  and  reported  the  time  to  reach  50%  of 
complete  desorption  (Ts0);  for  this  data,  a first-order  race  constant 
s 0.693/TM. 


been  quantitatively  analyzed.  In  these  cases,  a blconcinuujs  model  was 
employed  Co  determine  Che  appropriate  sorption  nonequilibrium 
parameters.  The  nonlinear  least-squares  optimization  program  CFITIM 
(van  Genuchcen,  1981)  was  used  to  fit  the  model  to  measured  breakthrough 

The  organic  chemicals  for  which  sorption  data  were  compiled  fall 
into  two  general  classes — hydrophobic  organic  chemicals  (HOCs)  and 
nonhydrophobic  organic  chemicals  (NHOCs).  HOCs  are  chose  whose  sorption 
is  considered  to  be  predominated  by  hydrophobic  sorption.  Hydrophobic 
sorption  is  generally  characterized  by  linear,  reversible,  and  singular 

Because  of  their  more  reactive  nature,  NHOCs  are  prone  to  specific 

The  paired  Kp  and  kj  values  used  for  the  analyses  presented  herein 


n Tables  4-1  and  4-2. 


re  divided  into  two  ty| 
nd  NHOCs  (Type  II;  see  Table  4-2). 
f experiment  used,  whether  the  data  wi 
s or  by  the  current  author,  and  the  si 


Results  and  pl?cug}ion 

Identification  of  mechanisms  responsible  for  sorption 
nonequilibrium  is  not  possible  by  curve-fitting  of  experimental  data 


Kp  (ml/g)  kj  (he  l)  Chemical  Exp.  Tech.  Simulation  Reference 


tetrachlorobenzene 

trichloroechene 

crichloroethene 

phenanthrene 

phenanthrene 

phenanthrene 

phenanthrene 

naphthalene 

naphthalene 


3 or  4-ethyl toluene*  i 


trimethylbenzene 


B - batch;  GP  - gas-purge-desorption;  C - column;  SA  - source  authors;  PA  - 
+ - value  used  in  Figure  4-2 


Kp  (Typo  II  Chemicals) . 


being  investigated  so  that  mechanisms  which 
may  be  identified.  With  these  Improbable  mi 
number  of  probable  mechanisms  is  reduced.  In  some  cases,  it  may  be 

■ . The  discussion  is  structured  by  type  of  nonequilibrium 
cder  to  facilitate  the  elimination  of  improbable 


Physical  noneouilibrium.  As  discussed  by  Brusseau  and  Rao  (1989), 

Klute  (1977),  Lee  et  al.  (1988),  Hiller  and  Weber  (1986),  Patterson  ec 
al.  (1987),  Rao  et  al.  (1979),  and  Wood  et  al.  (1987).  In  addition,  the 
particle  size  of  the  sorbents  used  in  all  these  studies  was  less  than  1 
mm.  It  can  be  readily  shown  that,  for  Che  tortuosities  typical  of 
natural  soils  and  sediments,  pore  diffusion  within  particles  of  this 
size  is  fast  enough  such  that  TNE  will  be  relatively  insignificant.  The 
BTC  for  the  non-sorbing  solute  of  Souchworch  et  al.  (1987)  'exhibits  some 
asymmetry.  However,  these  authors  successfully  simulated  the  BTC  with 

that  might  have  been  present  was  relatively  insignificant.  The  particle 


either  non-existent  or  of  negligible  importance  for  all  of  the  data 
utilized  herein.  Hence,  TNE  may  be  eliminated  as  a potential 


Because  of  their  chemical  s 

this  hypothesis  will  be  evaluated  in  light  of  the  discussions  presented 

CNE . The  literature  is  replete  with  investigations  of  specific 
interactions  between  these  types  of  chemicals  and  soils  and  sediments. 
Intrasorbent  diffusion.  Considering  the  nature  and  sorptive 

successfully  competes  for  polar  surfaces) , IOHD  is  a likely  potential 
nonequilibrium  mechanism.  As  previously  mentioned,  XOMD  has  been 

Conversely,  IMD  is  an  unlikely  mechanism  also  because  of  Che  reasons 

exhibiting  this  type  of  behavior  have  been  reported  by  several  authors 
(Siegel  ec  al.,  1951;  Harris  and  Uarren,  1984;  Moyer  ec  al.,  1972; 
Saltzman  et  al. , 1972;  Grover,  1974;  LaFleur  et  al.,  1978;  LaFleur, 


IOMD  may  be  a possible  nonequilibrium  m* 
nonequilibrium  for  Che  analyzed  data:  a)  HOGS — intraorganic  mad 


The  daca  in  Tables  4-1  and  4-2  are  presented  in  Figure  4-1.  Three 

linear,  validating  the  LFER  approach.  We  will  first  discuss  the  second 
feature,  returning  thereafter  to  the  first. 


: IOMD.  sorption  nonequilibrium 
is  the  result  of  race-limited  mass  transfer  of  sorbate  from  the  exterior 
surface  of  Che  organic  matter  into  the  interior  matrix.  Such 

associated  with  the  sorbent  matrix,  and  Che  diffusion  path  length.  For 
HOCs,  Kp  is  influenced  predominantly  by  two  components,  Kpc  (partition 
coefficient  normalized  to  organic  carbon)  and  foc  (fraction  of  organic 
carbon) . Koe  and  diffusivity  are  both  functions  of  the  size  and 


l°g  K 


relationship  between  k.  and  Kp  observed  in  Figure  4-1  is  expected  for 
Karlckhoff  (1980),  Karlckhoff  and  Morris  (1985),  and  Nkedi-Kizza  et  al. 


nonequilibrium  mechanism  for  Che  HOC  data  analyzed.  This,  in 
combination  with  Che  inverse  relationship  discussed  above,  is  strong 
evidence  for  Che  existence  of  IOMD  and  for  it  being  the  cause  of  Che 
nonequilibrium  exhibited  by  the  HOC  data.  To  further  test  the  validity 

in  Figure  4-2  (the  data  used  in  Figure  4-2  are  noted  in  Table  4-1).  A 
difference  of  three  and  one  half  orders-of-magnitude  exists  between  the 

nonequilibrium.  As  mentioned  above,  TNE  Is  a pore-diffusion  process. 

the  diffusing  molecule,  benzene  and  HCB  would  be  expected  to  have 
different  D values.  The  mass  transfer  coefficient  (a)  for'  the  TNE  model 

(Brusseau  and  Bao,  1989  and  references  cited  within): 


/ <r  t !) 


(4-1) 


where  t is  i shape  coefficient,  T is  tortuosity,  i is  Che  characteristic 
diffusion  half  length  (L) , D0  is  diffusion  coefficient  in  bulk  water 


'*60| 


(L!/T).  and  D - (D0/r) . Hence,  Che  Influence  of  Che  sice  of  Che  soluce 

soluces  resulc  only  from  differences  in  D„  values.  Therefore,  Che 
differences  in  Che  a values  for  bencene  and  HCB  can  be  evaluaced  by 

The  D0  values  for  bencene  and  HCB,  calculated  with  Che  Hayduk  and 
Laudie  correlation  (1974),  are  plotted  in  Figure  4-2  (noce  chat  Che 
values  are  scaled  to  Che  existing  ordinate  values  and  chac  a line  is 
used  in  place  of  Che  point  values  Co  facilicace  comparison) . IC  is 

che  slope  is,  of  itself,  support  for  Che  exiscence  of  IOMD. 

Noce  Chac  Che  data  in  Figure  4-2  Include  five  differenc  sorbents. 
The  high  correlation  (r2  - 0.99)  between  k2  and  Kp  suggests  chac  Che 

race— limited  process.  This  type  of  behavior  would  not  be  expected  for  a 

expected,  however,  for  a system  involving  a physical  process  such  as 
of  CNE  and  Che  existence  of  IOMD. 

(PAHs)  are  presented  in  Figure  4-3  (che  data  used  in  Figure  4-3  are 
greater  chan  che  slope  chac  would  resulc  from  a TNE-concrolled  system. 


for  Che  chlorinated 


in  Figure  4-3  represent  five  different  sorbents;  as  for  Che  chlorinated 
benzene  dace,  the  high  correlation  (r2  - 0.97)  supports  the 
improbability  of  ONE  involvement. 

The  slope  and  intercept  values  for  the  chlorinated  benzene  series 
and  the.  PAH  series  are  essentially  identical.  The  fact  that  two  series 
of  solutes  of  different  chemical  nature  exhibit  identical  behavior  is 


evocative  of  a physical  process, 
for  CHE- influenced  systems.  Give 


behavior  would  not  be  expected 
sea  presented  in  Figures  4-2  and 


relationship  between  log  k2  and  log  K,  exhibited  in  Figure  4-1  can  be 
quantified  with  linear  regression  analysis.  The  procedure  most  commonly 
employed,  bssed  on  Che  least-squares  model,  assumes  chat  the  independent 
variable  has  no  associated  error  (i.e.,  it  is  known  absolutely).  This 
assumption  is  invalid  in  our  case  since  measured  values  of  Kp  do  have  an 
associated  error.  Hence,  the  analytical  procedure  must  account  for 

Such  a technique,  called  the  geometric  mean  functional  (( 
method,  was  discussed  by  Halfon  (1985).  ' 
the  sum  of  the  products  of  the  vertical  ai 

sum  of  the  squares  of  the  vertical  deviations. 

The  following  regression  equation  relating  log  k2  ai 
the  Type-I  data  was  found  with  the  CHF  method: 


')  regression 
is  GKF  regression  minimizes 
l horizontal  deviations  of 
in  regression  minimizes  the 


d log  Kp  fc 


n empirical  relation  to 

or  Kp,  which  can  be 

9;  Karickhoff  et  ml.,  1979; 


- 0.301  - 0.668  log  K,  (r2  - 0.95)  (4-2) 

exhibited  by  the  Type-I  data  is  remarkable,  especially 
considering  that  the  data  were  obtained  by  many  Investigators  and  with 
different  experimental  techniques. 

Equation  4-2  should  prov* 

(solubility,  Kpp,  molecular  size,  e 
of  Che  sorbent  (e.g.,  Chiou  et  al. 

Sablj  1c , 1987)  . This  provides  the  capability  to  evaluate  Che 
nonequilibrlum-sorpcion  potential  of  a given  system  (i.e.,  sorbate, 
sorbent,  solvent).  Such  information  would  be  useful  in  the  selection  of 

Equation  4-2  also  provides  a means  Co  independently  determine 
for  dispersion  and  retardation  factor  values  (see  above) , would  allow 


of  chemical 


kinetic  behavior. 


Influe 


noted  above,  it  is  readily  apparent  Chat  Che  daca  in  Figure  4-1  cluster 
in  two  discrete  groups  according  to  chemical  type.  The  kj  values  for 

those  of  the  Type  I daca.  Thus,  Che  Type  II  chemicals  appear  to  be 
influenced  by  another  race-limiced  mechanism,  instead  of,  or  in  addition 
to,  the  I OHO  mechanism  affecting  Che  Type  I chemicals.  Considering  the 
nature  of  the  chemicals  comprising  the  Type  II  group  (ionizable  or  high 

possible.  In  fact,  it  has  been  suggested  chat  specific  chemical 
Interactions  can  occur  between  functional  groups  of  soil  organic  matter 
and  the  triazine  herbicide  atrazine  (Hayes,  1970;  Li  and  Felbeck,  1972), 
the  carboxylic  acid  herbicide  picloram  (Nearpass,  1976),  and 
chlorophenols  (Boyd,  1982;  Boyd  and  Mikesell,  1987).  Thus,  the  sorption 
nonequilibrium  exhibited  by  the  Type  II  chemicals  may  be  the  result  of 
CNC  or  a combination  of  CNE  and  I0HD. 

The  influence  of  the  chemical  nature  of  the  sorbate  molecule  on  the 
kinetics  of  interaction  with  organic  matter  can  be  further  evaluated  by 

2,4-D,  and  pencachlorophenol  (FCP),  is  controlled  by  the  molecular  form 
whon  pH  is  less  than  their  pK,  (Bailey  et  al.,  1968;  Ue  et  al.,  1988). 

conditions.  For  example,  based  on  experimental  results,  physical 


adsorpcion  (i.e.,  hydrophobic  sorpcion)  has  been  suggested  as  Che 
primary  sorpcion  mechanism  under  low  pH  condicions  (pH  less  chan  pK,) 
for  2,4-D  (Khan,  1973 ; 1974) . picloram  (Nearpass,  1976;  Bailey  eC  al,, 
1968),  and  PCP  (Lee  ec  al.,  1988).  Sorpcion  race  conscancs  for  chese 
ionizable  compounds  would  be  expecced  co  fall  near  Che  Type  I group  when 
cheir  behavior  is  controlled  by  Che  molecular  form  (i.e.,  when  Che  pH  is 
less  chan  the  pK,).  Conversely,  race  conscancs  obtained  in  systems  with 
pH  values  much  greacer  chan  pK.  would  be  expecced  co  fall  wichin  Che 
Type  II  group. 


1989))  ac  two  different  pH  domains  (Table  6-2),  one  below  Che  pK,  and 
the  ocher  much  larger,  are  ploCCed  In  Figure  4-4.  For  comparison,  the 


regression  equations  obtained  for  the  daca  presented  in 

fall  near  Che  Type  I group  as  expecced.  Also  as  expected,  the  rate 
constant  obtained  under  high  pH  condicions  is  wichin  the  Type  II  group. 
Th°  h;  values  for  ionizable  chemicals  such  as  these  would  be  expecced  Cl 


be  influenced  by  the  degree  of  ionizacion  and,  hence,  Co  depend  upon  Che 
pK,  - pH  relationship. 


It  should  be  noted  Chat  the  experlmenc  reported  by  Patterson  and 
Llebscher  (1987)  was  performed  with  a very  low  ionic  strength.  For 
organic  acids  such  as  pentachlorophenol . the  magnitude  of  Che  ionic 
strength  may  be  expecced  co  Influence  their  sorption  dynamics.  In  fact, 
Lee  (1989)  observed  the  sorption-kinetic  behavior  of  pentachlorophenol 
in  a high  pH.  relatively  high  ionic  strength  (0.01  normal)  system  to 


attributed  to  the  formation 


of  me  cal  cation/phenolate  Ion  pairs.  The  behavior  exhibited  In  Figure 
4-4  is  strong  evidence  that  che  Type  I chemicals  are  influenced  only  by 
10HD.  whereas  the  Type  II  chemicals  are  influenced  by  CUE  or  by  CUE  and 


Given  the  wide  range  in  physicochemical  properties  exhibited  by  NHOCs, 

widely.  In  particular,  it  is  expected  that  k,  values  for  NHOCs  would  be 
dependent  upon  the  nature  of  specific  sorbate-sorbent  interactions.  The 


kj  values  for  NHOCs  might 

that  Type  II  chemicals  wot 
rather  than  a single  ono  s 
relative  positions  of  the 


expected,  therefore,  to  exhibit  greater 
oose  for  HOCs.  As  such,  it  is  more  likely 
exhibit  a series  of  roughly  parallel  LFERs, 
is  shown  in  Figures  4-1  and  4-4.  The 
SRs  (l.e. , the  intercept  values)  would,  of 

interactions.  The  intercept  value  of  che  LFER  reflects  the  degree  to 

would,  therefore,  be  indicative  of  the  specific-interaction  potential  of 
the  system.  The  smaller  the  specific-interaction  potential  of  a NHOC 

LFER  would  be  to  that  of  the  HOCs.  The 


relative  paucity  of  kj  values  for  Type  II  chemicals  limits  the 
discussion  in  this  paragraph  to  reasoned  speculation.  Additional  data 
are  required  to  substantiate  the  postulated  behavior. 

So  for,  the  behavior  of  cationic  organic  species  (COCs)  has  not 
been  discussed.  Although  no  data  were  found  for  such  chemicals,  it 


might  be  speculated,  t 


u exchange. 


chac  k~  values  for  COCs  would  be  larger  chan  chose  for  HOCs  (i.e.,  the 
LFER  for  COCs  would  ploc  above  Che  HOC  LFER) . Cation  exchange  has  been 
demonscraeed  to  bo  an  Instantaneous  physical  process  involving  no 
chemical  bond  formation  (e.g. , Helfferich  (1962)),  Therefore,  ion 
diffusion  co/from  the  exchanger,  rather  than  Che  actual  exchange 
"reaction"  itself,  may  be  considered  as  the  rate-limiting  step  for  ion 
exchange  (Helfferich,  1962).  Considering  the  nature  of  soil  organic 
matter,  it  is  likely  that  a relatively  large  proportion  of  cation 
exchange  sites  associated  with  organic  matter  are  located  near  the 
soil/solution  interface.  Furthermore,  based  on  the  double  layer  theory, 
it  may  be  anticipated  that  Che  electroneutrality  requirements  for  ion 
exchange  can  be  satisfied  by  ions  residing  within  the  electrical  field 
of  the  sorbent.  As  a result,  the  effective  diffusion  path  lengch  for 
COCs  might  be  expected  to  be  shorter  than  that  for  HOCs.  The  kj  values 
for  COCs  would  therefore  be  larger,  since  kj  is  a function,  in  part,  of 
the  diffusion  path  length.  Because  of  the  nature  of  cation  exchange  and 
the  conformation  of  organic  matter,  as  described  above,  Che  ion  exchange 
of  cations  on  organic  matter  may  be  considered  as  essentially  a film- 
diffusion  process.  In  other  words,  any  rate-limited  behavior  exhibited 
in  such  systems  may  usually  be  attributed  to  film  diffusion.  This  is 
supported  by  the  work  presented  by  Bunzl  (Bunzl,  1974,  Bunxl  et  al., 
1976).  They  investigated  the  kinetics  of  ion  exchange  of  metal  cations 
and  concluded  chat  Che  race-limiting  step  was 
tsearch  on  NHOCs  and  COCs  is  required  Co  further 


film  diffusion, 
elucidate  what 


A preliminary  investigation  ol 

and  methods  employed  were  the  same 
the  appendix  for  a listing  of  the  : 

species,  considering  that  cation  exchange  is  expected  ti 
performed  in  the  manner  described  in  Chapter  3.  Properties  of  the 


:he  behavior  of  COCs  was  performed 
le  values  for  k2  and  Kp  obtained 


The  BTC  obtained  for  4sCa  is  shown  in  Figure  4-5,  along  with  the 

relative  pore  volumes  (i.e.,  R/p)  to  facilitate  comparison  of  relative 
nonideal  behavior.  It  is  apparent  that  while  45Ca  exhibits  some 
nonideal  behavior  (a  slight  tail) , this  nonideality  is  much  less  chan 
that  exhibited  by  naphthalene.  In  other  words,  the  sorption  of 


uojjejjuaouoo  SAiiB|ea 


Relative  Pore  Volumes  (R/p) 


cation  exchange  of  *sCa  is  much  saa 


is  suggested  t 


that  associated 


relatively  rapid  rate  of  ion  exchange  of  metal  cations  in  organic  natter 
was  also  reported  by  Bunzl  at  al.  (1976).  Such  behavior,  where  Che 
cation  (45Ca)  exhibited  a lesser  degree  of  nonequilibrium  chan  did  the 
HOCs,  was  also  observed  by  Nkedi-Kixza  et  al.  (1989). 

The  influence  of  nolecular  structure  of  Che  sorbate  on  sorption 

for  the  BTCs  are  provided  in  the  appendix, 

in  Figure  4-6.  It  is  apparent  that  each  series  is  described  by  a 
slightly  different  LFER.  However,  the  differences  are  so  slight  that 

Kp,  r2  - 0.96).  This  suggests  that  differences  in  behavior  resulting 

the  rate-limited  behavior  of  these  HOCs. 

The  influence  of  substitution  pattern  on  rate-limited  sorption  was 

relationship  between  k.  and  Kp  (see  Figure  4-7) . The  influence  of  a 


log  K 


toluene,  p-xylene,  1 . 3 . 5-crime thy lbenzene  series.  It  appears  ehac  che 
presence  of  che  functional  group  does  not  affect  the  nature  of  the  k.  - 
Kp  relationship  (see  Figure  4-7).  The  influence  of  multi-carbon 
function  group  was  evaluated  using  che  echylbenzene,  n-propylbenzene , u- 
bucylbenzene  series,  it  appears  chat  the  presence  of  a chain  of  three 
or  more  carbons  alters  the  k2  - Kp  relationship.  For  a process 
Involving  diffusion  within  an  organic  polymer,  this  behavior  may  be 
explained  as  resulting  from  the  additional  resistance  to  diffusion  that 


exists  for  the  bulkier  molecules.  In  othor  words,  che  presence  of  che 
alkyl  chain  hinders  the  diffusional  mass  transfer  of  n-propylbentene  and 
n-butylbenzene.  The  influence  of  molecular  structure  on  diffusion  in 
polymers  has  been  investigated  in  the  field  of  polymer  physics. 


The  influence  of  che  chemical  nature  of  the  sorbate  molecule  on  th 
rate  of  ics  interaction  with  natural  sorbents  was  recognized.  Strong 
indirect  evidence  was  presented  for  the  existence  of  the  intraorganic 
matter  diffusion  process  and  for  it  being  che  cause  of  the  sorption 
nonequilibrium  exhibited  by  the  HOC  data  analyzed.  An  equation  was 
presented  that  may  be  used  to  estimate  approximate  values  of  the 
sorption  rate  constant  for  HOCs.  This  provides  a means  to  evaluate 
nonequilibrium  potential  and  to  attempt  to  predict  nonequilibrium 
behavior.  Evidence  was  also  presented  chat  suggested  that  NHOCs  may  be 
influenced  by  some  combination  of  two  nonequilibrium  mechanisms, 
chemical  nonequilibrium  and  intraorganic  matter  diffusion.  It  thus 
appears  that,  in  many  cases,  the 


sorption  kinetics  of  organic  chemicals 


Such  processes  could  have  several  important  effects,  such  as 
reducing  the  effective  travel  time  of  a solute.  Increasing  the  time 
required  to  "flush"  a contaminated  aquifer,  causing  che  nonsingularity 
observed  in  many  sorption  experiments . and  being  responsible  for  the 
formation  of  so-called  "bound  residues". 

The  sorption  kinetics  of  lonizable  chemicals  appears  to  be 
significantly  influenced  by  pH.  Additional  work  is  required  to  exploro 


The  influence  of  organic 
of  hydrophobic  organic 


distribution  (i 


solubility  of  HOCs  in  ] 
with  increasing  volume 
The  log-linear  model  hi 
descriptor  of  che  influen 

A model  describing  ch 
equilibrium  sorption 


is  sorption  and  transport 

esearch.  The  Impact  of  cosolvents  on  phase 
ity,  sorption)  of  HOCs  has  been  che  major 

in  of  cosolvent  (f„)  in  a binary  mixture. 

..  1988;  Yalkowsky  and  Roseman,  1981). 
feet  of  cosolvent  upon  che  sorption  of 
1.  (1985).  The  relationship  between  the 


log  K„  - c f„  * d (4-3) 

where  K^  Is  Che  sorpcion  equilibrium  consCanc  for  Che  mlxed-solvenc 

che  cosolvency  power  of  Che  cosolvenc  expressed  via  Che  slope  of  Che  HOC 
solubillzacion  profile;  a is  an  empirical  conscanc  chac  represents  any 

solubility;  and  d-  log  . where  Kp*  is  che  equilibrium  conscanc  for 
sorpcion  from  water.  The  model  has  been  shown  Co  be  valid  for  boch 
sorpcion  (Nkedi-Kizza  ec  al..  1985 ; 1987 ; Fu  and  Luchy,  1986b)  and 
transport  (Nkedi-Kizza  ec  al.,  1987}  of  HOCs. 

The  influence  of  cosolvencs  on  che  kineclcs  of  sorpcion  of  HOCs, 
however,  has  received  very  liccle  accencion.  Nkedi-Kizza  et  al.  (1987) 

Increasing  fraction  of  cosolvenc  (fe)  for  displacemenc  of  che  herbicides 

between  BTC  asymmetry,  which  was  accribuced  by  che  authors  Co  sorpcion 

che  presence  of  a cosolvenc.  Nkedi-Kizza  ec  al.  (1987;1989)  suggested 

sorbenc  organic  matter  was  responsible  for  che  observed  behavior. 
Walters  and  Guiseppl-Elie  (1989)  reported  chac  sorpcion  equilibrium  of 
dioxins  by  soils  from  wacer/mechanol  mixtures  was  attained  in  shorter 

cosolvencs  on  Che  race  of  sorpcion  of  HOCs  has  been  qualicacive  in 
nature.  To  invescigace  and  evaluate  che  relacionship  between  sorpcion 


appears  Co  also  have  an  impact  on  Che  race  of  sorption.  The  existence 
of  a relationship  between  these  two  phenomena,  and  the  nature  of  the 
relationship  if  one  exists,  are  Che  central  questions  to  be  investigated 
herein.  The  critical  component  to  any  such  relationship  would  be  che 
existence  of  a relationship  between  che  sorption  race  constant  (kj)  and 
Kp.  Such  a relationship  was  presented  earlier  in  Chapter  4,  where  kp 
was  found  to  be  a log-log  linear  inverse  function  of  K.,: 

log  k,  - a log  K,,  + b (4-4) 


hr  and  Kp  has  units  of  ml  g*1.  Race-limited  interaction  between 
sorbent  organic  matter  and  the  sorbets , involving  diffusive  mass 
transfer  wlchin  the  organic  matter  matrix,  was  postulated  as  being 
responsible  for  the  observed  relationship. 

sorption  of  HOCs  may  be  developed  by  coupling  che  relationships  that 
have  been  previously  reported  between  kj  and  Kp,  and  between  Kp  and  fe. 
Note  chat  the  units  for  Kp  in  eq  4-3  and  4-4  differ.  To  attain 

log  kj  - a log  Kp’  + b'  (4_5 


where  b*  - b + a log  V;  V Is  Che  molar  volume  <ml  mol-1)  of  the  solvent; 

Substitution  of  eq  4-3  into  eq  4-5  results  in  the  following: 

log  - * f5  + r (4-6) 

where  lt„  is  the  sorption  rate  constant  (hr-1)  for  the  mixed-solvent 


A log-linear  functionality  between  k,  and  f=  is  expected  for  HOCs, 
based  on  eq  4-6.  Such  an  experimentally-determined  relationship  for  k; 
versus  fe  could  be  extrapolated  to  estimate  a k;  value  for  an  aqueous 
system.  If  the  functionality  between  k,  and  K,  were  Co  hold  for  mixed- 
solvent  systems,  che  regression  equation  they  reported  should  be 

mixed-solvent  systems  could  then  be  estimated  independently. 


Miscible  displacement  experiments  were  performed  with  naphthalene 
using  columns  packed  with  a sandy  soil  (Eustis  fine  send,  organic  carbon 
content  (OC)  - 0,39%).  The  experimental  apparatus  and  methodology  were 
similar  to  chose  employed  by  Hkedl-Kizza  et  ml.  (1987)  and  Lee  ec  al. 
(1988).  BTCs  for  phenanchrene,  naphthalene,  and  diuron,  obtained  by 
miscible  displacement  wich  mixtures  of  water  and  methanol  through 
columns  packed  with  Eustis  soil  (OC  - 0.72»)  as  reported  by  Hood  et  al. 
(1987),  were  used  as  an  additional  source  of  data.  A bicontlnuum  model, 
where  sorption  is  assumed  to  consist  of  both  an  equilibrium  domain  and  a 
race-limited  domain,  was  used  to  simulate  the  measured  BTCs  A 


Co  characteristic  time  of  sorption) . t 


M (van  Genuchcen, 
the  fraction  of 


contains  k; . Values  for  all  ocher  model  parameters  were  obtained  from 
independent  sources:  Peclet  number  (i.e..  dispersivity,  aL)  from  Che  BTC 
of  triciated  water,  a nonsorbing  solute;  retardation  factor  (R) , which 

sorbing-soluce  BTC  (c.f.,  Valocchi,  1985);  and  pulse  width  (T)  from  the 
experiment  itself.  Resultant  values  for  kp  snd  Kp  were  analyzed  Co 


Results  and  Piswsaion 

Evaluation  of  results 

The  BTCs  for  naphthalene  in  aqueous  solutions  (f0  - 0)  and  binary 
mixtures  (fc  - 0.1,  0.2,  and  0.3)  are  presented  in  Figure  4-8. 


with  increasing  fe  is  readily  apparent.  The  BTC  for  triciated  water 
(not  shown)  was  highly  symmetrical  and  was  described  with  a large  P 

exhibited  by  the  naphthalene  data  is  the  result  of  a sorption-related. 

To  test  Che  validity  of  eq  4-6,  log  k?  was  regressed  with  fc;  these 
results  are  presented  in  Figure  4-9.  The  log-linear  relationship 
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predicted  by  eq  4-6  appears  to  be  a good  dascripcor  of  tb 
Equation  4-6  was  further  tested  with  the  naphthalene  and 
data  obtained  from  Wood  et  al.  (1987).  Results  of  regressing  log  k2 
fc  are  presented  In  Figure  4-10.  Again,  the  log-linear  relationship 


eq  4-6.  The  question  then  becomi 
organic  matter  and  the  sorbate  hi 


that  the  rate  of  sorption  of 

,s  log-linear,  as  predicted  by 
t Is  the  mechanistic  explanation 
ted  Interactions  between  sorbent 
suggested  as  a primary  cause 


, 1988)  support 


experiments  (Nkedl-Klzta  et  al.,  1989;  8ouc 
this  proposal.  As  discussed  earlier,  soli  organic  matter  has  been  shown 
to  be  a complex  polymeric  system.  The  permeation  of  small  molecules 
through  polymers  Is  a function  of  several  factors  Including,  the  size, 
shape,  and  substitution  pattern  of  the  molecule  and  the  nature  and 
conformation  of  the  polymer  (Rogers,  1985;  Stannett,  1968).  The 
conformation  of  a polymer  is  In  turn  Influenced  by  the  nature  of  the 
local  environment;  for  example,  a polymer  will  expand  or. contract  In 
response  to  changes  In  solvent  polarity.  The  permeation  of  small 
molecules  la  affected  by  this  contractlon/expanslon.  The  conformation 
of  a hydrophobic  polymer  will  be  In  a relatively  condensed  state  when 
water  Is  the  sole  solvent.  A reduction  In  solvent  polarity,  as  achieved 
by  addition  of  an  organic  solvent  such  as  methanol,  will  result  In 
subsequent  expansion  or  -swelling-  of  the  polymer.  Concomitantly,  the 


expected  chat  0.  Che  fraccion  of  instantaneous  retardation,  would  also 
be  influenced  by  Che  presence  of  cosolvent.  Polymer  expansion,  by 

the  relationship  between  0 and  fc.  regression  analyses  were  performed 
for  two  published  data  secs.  A linear  relationship  between  0 and  fc  was 
found  for  Che  naphthalene  (r2  - 0.90)  and  phenanthrene  (r2  - 0.54)  data 
of  Hood  ec  al.  (1987)  and  for  the  acrazlne  (r2  - 0.96)  and  diuron  (r2  - 


6uaa».i)l«tton  I1 


Predicted  k„  (hr'1)  Measured  « Difference  Sorbate 


very  well  with  Che  experimentally-determined  value.  This  comparison  was 
also  performed  for  the  dluron  data  presented  by  Nkedi-Kizza  et  al. 
(1989),  who  obtained  hj  values  for  fe  values  of  0.0,  0.03,  0.1,  and  0.3. 
Again,  the  predicted  k„  value  comperes  very  well  with  the  measured 


die  results  presented  above,  it  appears  thac  reasonable 
,,  values  may  be  obtained  by  extrapolation  from 
experiments  performed  with  mixed  solvents.  Furthermore,  it  appears  ctla 


! should  bo  amenable  to  independent 
prediction.  To  test  this,  values  of  i and  t were  estimated  and  compared 
with  values  obtained  through  experimentation.  Recall  chat  ^ Is  a 
function  of  o,  «,  and  a,  the  slope  of  the  log  k;-log  K,  functionality. 


A regression  equation,  relating  o 


o octanol/water  partition 


and  Luchy,  1986b),  a value 


The  parameter  t is  a function  of  the  intercept  and  slope  of  the  log 
ka-log  Kp  relationship  and  the  intercept  of  the  log  Kp-fc  relationship 
(i.e.,  log  Kp.) . The  value  of  0.301,  reported  in  section  4-1  was  used 


Che  log  kj-log  Kp  incercepc. 


Karlckhoff  ec  al.  (1979),  relacing  et 
parcicion  coefficient  (K0,)  to  r.cw . wa 


e organic  carbon-normalized 

5 used  to  estimate  values  for  Koc , 
from  which  1C*  values  were  calculated  (i.e.,  Kp,  - Kp.  OC) . 

The  predicted  and  c values  for  the  naphthalene  data  are  2.06 
and  -0-.20,  respectively  (an  experimentally-determined  value  of  2,300 
(Karlckhoff,  1981)  was  used  for  K„.).  These  values  compare  relatively 

and  naphthalene  data  of  Wood  et  al.  (1987),  using  an  experimental  Kp, 
value  of  37,000  (Karlckhoff,  1981)  for  phenanthrene . The  predicted 

0.37  for  naphthalene.  Again,  these  values  compare  well  with  the 
experimentally-determined  values  presented  in  Figure  4-10.  It  should  be 

Deviations  from  ideal  behavior 

upon  the  Kp-fc  functionality.  The  cosolvency  effects  of  weaker 

il  generally  follow  the  log-linear  model; 


cosolvents  si 


however,  the  effects  of  stronger 
observed  to  deviate  from  log-linearity  (Rao  ai 


Yalkowsky , 19B7).  Hence.  Che  k,-f,  functionality  may  be  expected  Co 
deviate  from  log-linearity  for  systems  Involving  stronger  cosolvents. 
The  slope  of  che  k2-fc  functionality,  i , is  dependent  upon  the 
solute-cosolvenc-sorbent  combination  specific  to  a given  system.  For 
example,  d would  be  expected  to  Increase  with  decreasing  cosolvent 
polarity.  Thus,  stronger  cosolvencs  would  induce  greacer  effects; 

nonideal  behavior.  Therefore,  the  use  of  weaker  cosolvents,  such  as 
likely  to  follow  ideal  behavior  is  preferred, 
far,  has  centered  on  HOCs,  The  influence  of 
>f  sorption  of  HOCs  containing  "reactive" 

pacific  interactions  with  a sorbent,  cosolvents 
ace  of  sorption  of  these  chemicals  in 
erved  for  HOCs  not  containing  reactive  functional 
tionship  between  kj  and  fc  (r!  - 0.83)  for  a 
tern  was  reported  by  Nkedi-Klzza  ec  ml.  (1989).  To 
of  eq  4-6  for  the  diuron  data  of  (Nkedi-Klzza  et 
regressed  with  ft;  the  resultant  r1  value  of  0.77 
srted  for  the  k;-f„  regression.  The  hypothesis 
l diuron  data  reporced  by  Wood  ec  al.  (1987). 


potential  to  undergo 

comparison  to 

diuron/Eustis 
ovaluate  the  validity 
al.,  1989),  log  kj 


was  furthor 

better  descriptor  of  the  k2-ft  relationship 
vs  0.87,  respectively).  These  results  supp, 
sorption  kinetics  of  HOCs  containing  react!' 


n log-linear  (r2  - 0.93 
the  hypothesis  that  the 


containing 


reactive  moieties. 

The  influence  of  cosolvent  on  the  k;-fc  relationship  for  HOCs 
containing  functional  groups  vas  further  investigated  by  analyzing  the 
atrazine  data  reported  by  Nkedi-Kizza  et  al.  (1989).  No  relationship 

of  regression  analysis  (r*  - 0.02).  This  suggests  tl 

presence  of  methanol.  Interestingly,  as  reported  ab 
linear  relationship  with  f0  for  atrazine,  just  as  it  did  for  the  "non- 

increases  with  Increasing  fe,  the  race  of  sorption  of  atrazine  is 

that  observed  for  the  "non-reaccive"  HOCs.  Such  behavior  would  be 

and  McIntosh,  1971;  Sullivan  and  Felbeck  (1968). 


led  on  Che  results 
Sac  Che  mechanism 
c affected  by  the 
ove,  & exhibited  a 


The  results  presented  above  suggest  chat  data  on  sorption  kinetics 
collected  in  mixed  solvents  can  be  extrapolated  to  predict  the  kinetics 
of  HOC  sorption  from  aqueous  solutions.  Mixed  solvents  may  be  utilized 
to  facilitate  experimentation  with  organic  chemicals,  especially  chose 

sorption  to  apparatus  surfaces.  Obviously,  there  will  be  some 


uncertainty  associated  with  parameter  values  estimated  by  extrapolation. 
However,  the  magnitude  of  this  uncertainty  may  be  no  greater,  and  in 

values  obtained  from  direct  experimentation  with  very  hydrophobic 

ramifications  with  respect  to  modeling  the  transport  and  fate  of  organic 

when  modeling  solute  transport  (e.g. ; Valocchi,  1985).  If  sorption 
equilibrium  is  attained  more  rapidly  in  systems  containing  a cosolvent. 


Summary 

A model  describing  the  influence  of  organic  cosolvents  on  the 
sorption  kinetics  of  organic  chemicals  was  presented.  This  model 
predicts  the  sorption  rate  constant  of  hydrophobic  organic  chemicals  to 
be  a log-linear  function  of  the  fraction  of  cosolvent  present  in  the 
system.  The  experimental  data  presented  appear  to  substantiate  the 
validity  of  the  proposed  model.  The  mechanism  responsible  for  this 
phenomenon  was  postulated  to  involve  changes  in  conformation  of  the 
organic  matter  matrix,  which  are  induced  by  changes  in  solvent  polarity 
that  result  from  the  addition  of  a cosolvent.  The  diffusive  mass 
transfer  of  sorbate  within  the  organic  matter,  the  apparent  cause  of  the 

conformation  of  the  organic  matter  matrix.  Given  this  mechanism,  the 


fraction  of  instantaneous  retardation  (i.e.,  sorption)  would  be  expected 
to  increase  as  the  fraction  of  cosolvent  increases.  This  was  indeed  the 

Values  of  the  regression  coefficients  for  the  relationship  observed 
between  sorption  rate  constant  and  fraction  cosolvent  for  the 
naphthalene  and  phenanthrene  data  were  predicted  with  simple  empirical 

with  mixed-solvent  data.  This  supports  the  validity  of  extrapolating 
results  obtained  from  mixed-solvent  systems  to  aqueous  systems.  The 

valid  for  mixed-solvent  systems.  Thus,  the  relationship  represented  by 

Deviations  from  log-linear  behavior  are  expected  under  certain 
conditions.  For  example,  the  sorption  kinetics  of  chemicals  that  have  a 
greater  potential  to  specifically  interact  with  the  sorbent  might  be 

data  sets,  two  Involving  the  pesticide  dluron  and  one  the  pesticide 


MODELING  SOLUTE  TRANSPORT  INFLUENCED  BY 
MULTIPROCESS  NONEQUILIBRIUM 


reported  chec  suggests  thet  the  transport  of  solutes  con  be  Influenced 
by  tronsport-releted  ond  by  sorpclon-relsted  nonequilibriuin  processes. 
It  Is  logical  to  assume  that  in  some  cases,  especially  in  the  field, 
solute  transport  will  be  constrained  by  both  of  these  nonequilibrium 
processes.  When  investigating  solute  transport  it  is  imperative  that 
the  potential  for  multiple  sources  of  nonequilibrium  be  recognized  and 
evaluated. 

Nonequilibrium  is  often  conceptually  modeled  using  a bicontinuum 
approach,  wherein  the  porous  medium  Is  discretized  into  two  domains. 
The  so-called  "two-compartment" . "two-region",  or  "two-site" 
mathematical  models  (von  Genuchten  and  WIerenga,  1976;  Selim  et  al. 
1975;  Cameron  and  Elute,  1977;  Karickhoff,  1980)  chat  result  from  this 
approach  may  be  used  to  represent  all  three  nonequllibrium  processes 


Rao.  1989).  Such  models  have  b 

and  WIerenga,  1977;  Karickhoff 
Coates  and  Elzerman,  1986;  Mill 


assumption,  implicit  or  explicit,  that  a single  process  ii 
Cor  the  observed  nonequilibrium.  As  discussed  al 
to  be  situations  where  nonequilibrium  is  caused  by  more  chan  one 
process.  A model  chat  accounts  Cor  multiple  sources  oC  nonequilibrium 
would  be  required  to  accurately  represent  such  systems.  For  example, 


nonequilibrium  processes  would  yield  global  or  lumped  kinetic  terms; 
i.e.,  parameters  chat  represent  the  combined  eCCects  oC  all  processes 
contributing  to  nonequilibrium.  Conversely,  the  application  oC  a model 
designed  Cor  multi-process  nonequilibrium  would  yield  discrete  kinetic 
terms,  each  representing  individual  nonequilibrium  processes. 
Mathematical  models  are  sometimes  used  Co  assist  in  the  evaluation  o C 
the  conceptual  models  which  they  represent.  For  mathematical  models  tc 
be  employed  in  this  Fashion,  each  oC  the  processes  oC  concern  should  be 
represented  by  a discrete  term.  The  use  oC  global  or  lumped  terms, 
while  satisCaccory  Cor  purposes  oC  calibration  or  Fitting,  is  not 
conducive  Co  process  investigation.  Hence,  a model  that  accounts  Cor 
multiple  sources  oC  nonequilibrium  s 
use  in  process-level  investigations 
An  ultimate  goal  oC  modeling  is 
independently.  Such  capability  means  chi 
or  Fitting  parameter  values  is  abrogated 


>e  oC  greater 


'at  the  necessity  Cor  calibrating 


C system  response  to  speciCied  perturbations  are  possible. 
Global  parameters  ore  not  amenable  to  independent  determinations, 
because  they  represent  multiple  processes.  They  are,  in  essence 


fleeing  parameters  unique  Co  ehe  specific  systems  and  sec 
noe  particularly  useful  for  independent  prediction  of  solute  transport 

explicitly  accounts  for  multiple  sources  of  nonequilibriua  Influencing 
solute  transport  in  porous  medis.  The  conceptual  and  mathematical 
framework  of  Che  model  is  described.  Model  parameters  and  their 
determination  are  discussed.  A sensitivity  analysis  is  performed  to 

on  system  behavior.  Finally,  the  model  is  applied  to  several  published 
data  sets  to  illustrate  its  use  and  to  evaluate  ehe  validity  of  the 
conceptualization  upon  which  it  is  based. 


h transport-related 


nonequilibrium.  The  sorption  dynamics  of  such  systems  are  represented 
by  a model  consisting  of  two  serially-arranged  biooncinuums  coupled  in 
parallel.  A schematic  of  the  model  conceptualisation  is  shown  in  Figure 
5-1.  This  conceptualization  results  in  discretisation  of  the  porous 
medium  into  four  sorption  domains,  where  instantaneous  sorption  occurs 
in  the  first  domain  and  rate-limited  sorption  occurs  in  the  other  three. 
An  example  system,  where  nonequilibrium  may  be  caused  by  PNE  and 


intraorganic  matter  diffusion,  wi 
appreciable  organic  naccer  concei 

Che  aggregates  (sorpcion  domain  : 


coupled  by  diffusive  mi 


coupled  by  diffusive  mi 


d be  an  aggregaced  soil  with  an 

he  inceraggregace  (or  mobile)  region 
er  sicuated  on  external  surfaces  of 
Sorption  domains  1 and  2 are 

e (intraaggregace)  regions  are 

(i.e,,  TNE) . Transfer  between  the  immobile  region  and  sorpcion  domain  3 
is  instantaneous,  as  it  was  between  the  mobile  region  and  domain  1.  The 
coupling  between  sorption  domains  3 ond  4 is  similar  to  that  between 
domains  1 and  2.  Note  that,  although  mass  transfer  between  the  immobile 
region  and  domain  3 is  instantaneous,  accessibility  of  sorpcion  "sites" 
in  domain  3 is  rate  limited  by  mobile-immobile  mass  transfer.  Also, 
note  that  sorpcion  in  domain  4 is  constrained  by  two  mechanisms,  TNE  and 
intraorganic  matter  diffusion. 

The  TNE  component  in  the  MPNE  model  con  be  used  to  represent  any 
transport-related  nonequilibrium  process.  For  example,  in  the  case  of 
macroporous  soil,  the  macropores  can  be  designated  as  the  mobile  region 
and  the  bulk  soil  matrix  as  the  immobile  region.  The  sorption-related 
nonequilibrium  component  may  represent  intraorganic  matter  diffusion 
intramineral  diffusion,  or  chemical  nonequilibrium.  For  chemical 
nonequilibrium.  Che  process  should  be  one  chat  can  be  modeled  with  a 
serial  representation  (e.g.,  "physisorption"  which  is  generally  fast, 
n rate-limited). 


followed  by  "chemisorption”,  w 


alternate  conceptualization 


3e  formulated  by  linking  domains  1 
and  domains  3 and  u in  parallel  t 


and  chemical-related  processes  cause  nonequilibrium  and  where  Che 
chemical  nonequilibrium  process  is  best  modeled  with  a parallel 
representation  (e.g.,  sorption  by  two  different  types  of  sites,  where 
sorption  is  rate-limited  at  one  and  instantaneous  at  the  other).  The 
formulation  illustrated  in  Figure  5-1  was  selected  because  incraorganic 
matter  diffusion  was  considered  to  be  the  more  predominant 
nonequilibrium  process  for  HOCs,  in  comparison  to  chemical 


Ac  equilibrium,  sorption 


respectively.  The  Cotal  s« 


■ . Sr  (mess  sorbsce/cocel 


comprising  che  mobile  region. 


Si  - « (Su  + S^)  + (l-£)  (: 
where  f is  che  mass  fraction  of  sorb 
Ac  equilibrium, 

S,  - f K.  C.  + (1-f)  K„  C„ 

Recognizing  that  C„  - c„  at  equilibrium,  eq.  (3)  reduces  to 
sr  - <s.i  Xp  - If  «.  * (1-f)  Kin) 
where  Kp  represents  the  weighted-average  equilibrium  sorption  ct 


Sorption  dynamics  in  both  che  mobile  and  immobile  regions  are 
described  as  follows: 

asM 

TF  - X.1  S.1  - Ha  S^  - Iqa  [(1-F.)  K.  C.  - S^l  (5-5c) 

HT  ‘ k‘-i  Swi  * k»!  s»!  ■ *ua  ( (1-Fu.)  K„  C,.  - Sum]  (5-5d) 
where  k.s,  Iq^,  klBl , and  kioC  are  first-order  sorption  kinetic 
coefficients  (1/T). 

The  one-dimensional,  advecclve-dispersive , solute-transport 
equation,  formulated  in  terms  of  the  MPNE  model,  is  written: 


ssa 


asH 


1-Op 


~ ' *•  D ~ i ST 


volumetric 


and  Immobile  regions,  respectively;  p is  the  bulk  density  (H  L‘J)  of  the 

hydrodynamic  dispersion  coefficient  in  rha  mobile  region  (L’/T) : q is 
Darcy  flux  (L/T);  and  v - q/f  is  average  pore-water  velocity  (L/T).  It 
is  assumed  that  advective-dlspersive  mass  transport  occurs  only  in  the 
mobile  region.  All  solute  mass  in  the  immobile  region  is  subject  to 
mass  transfer  between  Che  mobile  and  immobile  regions: 


e a is  a first-order  mass  cransfer  coefficient  (1/T). 

Equation  5-6  may  be  rewritten  by  substituting  in  eqs  5-5  and  5- 
SC. 

<».  + CP  p.  5T  + fp  Km  I (1-P.)  K.  C.  - S*)  ♦ a (C,  - C„)  - 


fp  F„  SC,'  fp 

+ — W'TTr  V «V  - S.*)|  + w (€.•  - C„‘) 
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ax  (5-9) 


by  defining  the  following  dimensionless 
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vher«  L is  the  length  of  interest, 
Che  encice  porous  medium,  and  C0  is 
paramecers  described  in  equations  ! 
Damkohler  numbers  representing  the 


it  concentrecion.  Note  chat  tl 


three  nonoquilibrlum  terms  included 


id  may  be  rewritten  in  dimensionless  c< 


IT  ‘ - s-'> 

- kj»*(C,.*  - S„‘) 

Equations  5-9  and  5-10  may  be  simplified  to 
SC.’ 

*1  ~sT  * *■»  (C“"  - S*‘>  ♦“<«*-  ?„*)  - 

1 aC-' 


(9-12) 

(5-13) 


(5-14) 

(5-15) 


R=!  - U-F„>  K. 


<1-0* 

— (1-F„> 

From  eqs.  (5-4)  and  (5-16),  note  chat  Che  cocal  recardaclon 

R-R«  + ^tRuu  + Rua-Il  + <*/* 


(5-17) 


h - -j- 


(5-18a) 


fraction  of  recardacion  thac  is  accribuced  c 


is  based  on  Che  linear  driving-force  approximation,  firsc  presented  by 
Gleuckauf  and  Coates  (1947) . An  alternative  representation  for  PNE 
could  be  obtained  with  the  use  of  diffusion  equations  based  on  Fick's 
law.  A great  deal  of  work  utilizing  the  diffusion-based  approach  has 


was  used  by  Rao  and  covorkers  Co  describe  solute  transport  in  aggregated 
media  (Rao  et  al.,  1980a, b;  Wkedi-Kizza  et  el.,  1982)  and  more  recently 
by  several  others  (e.g.,  Crittenden  et  al.,  1986;  Hutzler  et  al.,  1986; 
Roberts  et  al.,  1987;  Miller  and  Weber,  1988).  Use  of  Che  diffusion- 
based  approach  requires  detailed  knowledge  of  the  physical  structure  of 
the  porous  medium  and  is  complicated  by  several  factors,  including 
complex  structural  shapes,  size-distributions  of  the  structures,  end  the 
presence  of  poorly-defined  structures.  These  and  ocher  factors 
constrain  the  number  of  cases  to  which  diffusion-based  models  can  be 
successfully  applied.  For  this  reason,  the  mass-transfer  approximation 
was  used.  Several  investigators  have  shown  that  models  utilizing  the 

obtained  from  diffusion-based  modols  (c.f.,  Rao  et  al.,  1980b;  Raghavan 
and  Ruchven,  1985;  van  Genuchcen,  1985;  Parker  and  Valocchi,  1986). 

A computer  program  developed  by  Jessup  (1988)  was  used  to  solve  the 
pertinent  governing  equations  (5-12  to  5-15)  under  the  following  initial 

C-’  (X'°>  - cui'  (*.<»  - 0 (5— 19a) 

S-‘  (X'0)  - Si.’  <*.0)  - 0 (5— 19b) 


(5-19d) 


*C.‘  (X.p) 


5—l9e) 


Performance  of  che  numerical  solution  was  checked  againsc  an  analytical 
solution  (Brenner,  1962)  for  several  simplified  cases  (Jessup,  1988). 

The  computer  program  developed  by  Jessup  (1988)  was  based  on  Che 
finite-difference  method  and  utilized  the  Crank-Nlcholson  technique  with 
equal  weighting  for  the  time  component  of  the  difference  equation.  The 

equation.  The  discretisation  (i.e.,  number  of  grids)  for  che  spatial 
component  was  optimised  based  on  che  value  of  the  Peclet  number;  the 
number  of  grids  ranged  from  20  co  50.  The  discretization  for  che 
temporal  component  (i.e.,  pore  volumes)  was  optimized  based  on  Che 
values  of  several  parameters  of  the  governing  equations  (Jessup,  1988). 

is  being  employed.  As  has  been  widely  discussed  in  che  literature,  this 
type  of  boundary  condition  is  more  appropriate  chan  a concentration-type 


Parker,  1986).  It  is  assumed  for  this  model's  formulation  that  sampled 
effluent  represents  mobile-region  solution,  or  in  other  words, 


represents  flux-averaged  concentrations.  The  boundary  condition  Imposed 
in  eq.  5-19e  is  that  employed  by  Brenner  (1962)  and  allows  the 
determination  of  flux-averaged  concentrations  at  X - 1.  It  has  been 
shown  that  simulations  produced  with  a solution  employing  this  boundary 
condition  are  similar  Co  those  produced  with  a solution  employing  che 
finite  column  with  no  boundary  effects  accruing 


Co  Che  endplate,  foe  P groacer  chan  approxlmacely  5 (van  Genuchcen  and 
Parker,  1984). 


be  Che  nose  difficulc  Cask  In  modeling.  One  mechod  used  roucinely  is 
Che  curve-ficcing  of  model  simulaclons  Co  experlmencal  dace,  wherein  Che 

coupled  Co  a nonlinear,  lease-squares,  opcimizacion  program  (UHAUS; 

fashion.  A model's  concepcual  basis  Is  cruly  cesced  and  valldaCed  only 
when  Che  model  can  simulaCe  a sec  of  experlmencal  daCa  wlch  all 
paramecer  values  decermined  independencly  of  Che  daca  being  simulaced. 

A general  scheme  for  obcainlng  Che  paramecer  values  required  for  che 
MPNE  model  is  presenCed  below. 


cransporc.  Hydrodynamic  dispersion  is  generally  considered  Co  be 
comprised  of  cwo  componcncs,  molecular  diffusion  and  mechanical 
dispersion.  The  concribucion  of  diffusion  Co  hydrodynamic  dispersion  is 
generally  insignificanc,  especially  in  Che  range  of  velocicies  cyplcal 
co  many  cases  of  concern.  As  a resulc,  all  soluces  experience 
essenclally  che  same  dispersive  cransporc.  The  P value  decermined  from 
che  breakchrough  curve  of  a nonsorbing  soluce  may  cherefore  be  used  for 
sorbing  soluces.  In  using  che  MPNE  model  ic  is  assumed  choc  PNE  is 
presenc.  Under  such  condiclons.  che  BTC  of  a nonsorbing  soluce  will 


1989).  A bicontinuum  model  a! 


•i  m 


i been  much  discussion  in  Che  liceracure  on  Che  validicy  of 

constrained  Co  be  as  similar  as  possible  Co  chose 
calculated  from  sorption  isotherms  may  be  successfully  used  for  column 

1989;  Lee  et  al.,  1988).  Other  means  of  determining  R values  Include: 

(1)  measuring  the  area  above  Che  BTC  for  a frontal  inpuc  (c.f.,  Roberts 
et  al.,  1980;  van  Cenuchten  and  Parker,  1984;  Nkedi-Klzza  ec  al.,  1987); 

(2)  determining  che  first  cemporal  moment  for  a pulse  input  (c.f., 
Kucera,  1965;  Valocchi,  1985);  and  (3)  che  use  of  various  ci 
equations  relating  K,  to  properties  of  che  sorbent  am 
Chiou  ec  al.,  1979;  Lyman,  1982;  Karlckhoff,  1984). 


nonequilibrium.  When  TNE  ! 
will  exhibit  nonideality  ar 

determined  by  fitting  che  t 


present,  the  BTC  of  a nonsorbing 


a pore-diffusion  process,  the  u 


differences 


diffusivieies 


The  firsc-order  mass-transfer  coefficient  (a)  for  Che  PNE  term  ci 
following  form  (Brusseau  and  Rao,  1989  and  references  cited  therein): 


where  D0  is  Che  diffusion 
factor:  a is  a geometric 
diffusion  length.  For  po 

the  only  solute-specific 

determined  from  that  of  a 
equation: 


5 (5-21) 

where  is  the  first-order  mass-transfer  coefficient  calculated  from 
Che  u value  determined  for  the  nonsorbing  solute  and  where 

1 ~ W®0o  (5-22) 

(s  and  n representing  the  diffusion  coefficients  for  the  sorbing  and 
nonsorbing  solutes . respectively) . Note  chat  the  corrected  value  of  u 
for  a sorbing  solute  will  usually  be  smaller  chan  the  u value  for  a 
nonsorbing  solute,  since  D„,  < and  T < 1;  i.e.,  the  diffusion 
coefficients  for  sorbing  solutes  are  generally  smaller  than  chose  of 
solutes  employed  as  tracers  (e.g. , triclated  water,  chloride). 


(SNE)  Co  the  cocal  nonequi librium  Influencing  soluCe  transport, 
independent  determination  of  sorption-kinetic  parameter  values  has 

values  is  an  independent  experiment  for  sorption  kinetics,  such  as  the 
batch-rate  or  the  gas-purge  desorption  methods. 

Another  approach  would  be  the  use  of  an  empirical  correlation 
equation  that  relates  the  sorption-rate  coefficient  to  the  equilibrium 
sorption  coefficient,  as  presented  in  section  4-1.  The  regression 
equation  reported  in  Figure  4-1  for  the  Type-I  data  may  be  used  to 
estimate  sorption  kinetic  coefficients  £kz)  for  HOCs.  These  values  Chen 

regression  equation  reported  for  the  Type-11  data.  However,  the  k; 


values  estimated 


Genuchten,  1981; 


1 have  greater  uncertainty  because,  a. 
e likely  exhibit  a series  of  LFERs  ra 


noted  that  the  definition  of  k*  as  given  in  eqs.  5-llh 
equivalent  to  the  definition  used  for  the  equivalent 
in  the  two-site  model,  which  is  k2I*(l-0)/v  (c.f.,  von 
Lee  et  al.,  1988).  The  definition  of  k*  can  be 


modified  in  Che  following  manner: 
k*  - k'/Rfft  + ft) 


(5-23) 


Co  atcain  equivalency  with  Chat  of  Che  two-sice  model.  This 


determine  independently.  0, 
recardaclon  (i.e..  sorption). 


could  be  used  Co  escimace  pl. 


lardacion,  are  Che  most  difficult  Co 
represents  Che  fraction  of  inscancaneous 
For  an  aggregated  soil,  Che  amount  of 
ill  time  interval  in  a batch  experiment 
To  independently  determine  values  for 


F's  could  be  determined  by  fitting  a bicontinuum  model  to  che  results 
physical  properties  of  che  porous  medium. 


Sensitivity  Analysis 

transport  equation,  in  comparison  to  chat  imposed  by  che  models 
employing  the  biconclnuum  nonequilibrium  conceptualization  and 
especially  the  local  equilibrium  assumption  (LEA).  It  is  important  to 
identify  the  conditions  under  which  the  use  of  the  MPNE  model  is 

adequate.  This  endeavor  is,  in  essence,  one  of  identifying  the 
conditions  under  which  che  LEA  is  valid.  Consider,  for  example,  a case 
where  nonequilibrium  is  caused  by  a transport-related  and  a sorption- 


related  process,  when  extant  conditions  are  such  that  the  LEA  is  valid 
for  one  of  Che  two  nonequilibrium  processes,  Che  MPNE  model  reduces  Co 

processes , a LEA-based  model  would  be  appropriate . 

Identification  of  the  conditions  under  which  the  LEA  is  valid  has 
recently  become  a focus  of  attention.  A component  of  this  research  has 
bean  the  development  of  criteria  that  may  be  used  for  evaluating  LEA 
validity.  The  Damkohler  number  has  been  the  most  widely  used  criterion, 
which  is  not  surprising  given  chat  the  Damkholer  number  is  a ratio  of 
hydraulic  residence  time  to  reaction  time  and.  as  such,  characterises 
the  degree  of  nonequilibrium.  The  use  of  Damkohler  numbers  as  LEA- 
validity  criteria  has  been  discussed  by.  among  others,  Helnyfc  at  al. 
(1983),  Rao  and  Jessup  (1983),  Jennings  and  Kirkner  (1984),  Valocchi 
(198S),  and  Bahr  and  Rubin  (1987).  The  use  of  Damkohler  numbers  as  In- 
validity criteria  has  been  criticised  by  Valocchi  (1985)  on  Che  basis 
thac  parameters  ocher  than  Damkohler  numbers  can  influence  LEA  validity. 

To  quantitatively  evaluate  LEA  validity  one  must  have,  in  addition 
to  a criterion,  a means  to  measure  deviations  from  "ideal"  behavior 
(i.e.,  that  exhibited  under  conditions  where  LEA  is  valid).  A method 
based  on  comparison  of  time  moments  of  the  solute  BTCs  was  presented  by 
Valocchi  (1985)  and  Parker  and  Valocchi  (1986).  A quantitative  error 
ss  defined  as 

- O^s  - lO/U£„  (5.24) 

is  the  n-th  central  time  moment  and  the  superscripts  K and  E 
° the  "'’"equilibrium  and  LEA  models,  respectively.  Inspection  of 


eq  5-24  shows  Chat  Eq  will  increase  as  departure 
local  equilibrium  increases  and  char  E„  - 0 when 
The  degree  of  "spreading"  or  dispersion  is  c 


central  temporal 
Lta  pulse,  was  given 


by  Parker  and  Valocchl  (1986) 

U2  - (2R*/P)  + (T!0/12) 

In  eq  5-25,  the  dispersion  Influencing  solute  transport  is  seen  to  be 
comprised  of  two  components,  hydrodynamic  dispersion  and  that  related 
a finite  input  pulse.  The  equation  for  the  second  central  temporal 
moment  for  the  MPNE  model  may  be  written  in  a manner  similar  to  that 


(5-25) 


employed  by  Valocchl  <1985)  for  the  bicontinuum  model: 

0,  - <2rVp>  + <2R**V«)  + <2R‘nVk-.)  ♦ (2RV.A-J  + <TVl2)  (5-26) 


in  the  MPNE  model,  along  with  those  related  to  hydrodynamic  dispersion 
and  finite  pulse,  contribute  to  observed  dispersion.  This  linear 
additivity,  where  total  dispersion  is  the  sum  of  the  contributions  of 
each  dispersion-causing  component,  has  long  been  recognised  in 


chromatography  (c.f. . Giddings,  1965). 

The  quantitative  error  criterion  for  the  MPNE  model  is 
Ej  - PiflVw  ♦ /lVk"„  ♦ flVk'i.] 

when  the  Influence  of  T0  is  neglected.  Inspection  of  eq  5-27 
LEA  validity,  as  measured  by  B,,  is  a function  of  P,  the  0's 
representing  the  three  rate-limited  sorption  domains,  and  the 


Damkohler  numbers.  This  concurs  with  the  analyses  presented  by  Valocchl 
(1985)  and  Parker  and  Valocchl  (1986).  Not.  that  the  retardation 


explicitly  influenced  by  R 


region  i)  Co  R,  which  i 


b recio  of  Rj  (recardacion  associaced  with 
iel  parameters  on  Che  validity  of  the  LEA  cai 


is  directly 


assessed  with  eq  5-27.  As  evident  from  eq  5- 
Valocchi  (1985) , Che  magnitude  of  deviation  i 
proportional  to  P.  This  functionality  is  a result  of  the  manner  in 

5-27  is  an  absolute  measure.  It  can  be  shown  that  LEA  validity  is.  in 
fact,  independent  of  P if  a relative  approach  is  employed.  Consider, 
for  example,  a case  where  P - 10,  u - 1,  - 0.25,  and  k*a  - k\B  - 10; 

Che  associated  Ej  value  calculated  with  eq  5-27  is  0.875.  This  value 

expected,  E;  is  larger  for  Che  caso  where  the  additional,  sorption- 
related  parameter  contributes  to  nonequilibrium.  The  deviation  between 
the  two  Ej  values  is  40%.  Now  consider  two  additional  cases  where  all 
parameter  values  are  Identical  to  the  previous  two  cases  except  for  the 
value  of  P,  which  is  now  set  equal  to  100.  The  resultant  Ea  values  are 

respectively.  The  magnitude  of  the  Ej  values  for  P - 100  are  greater 


larger  values  i. 


o smaller  values. 


validity  of 


the  0's  representing  the  three  nonequilibrium  domains 


numbers  may  be  expected  to  r 
than  the  range  expected  for 
the  parameter  chat  has  Che  g 


three  0's  is  0.1  to  0.5.  The  Darakohler 
i from  0.01  to  100,  which  is  much  greater 

:est  potential  for  variation  when 
Hence,  in  the  sensitivity  analyses 

he  MPNE  model  for  k°  is  presented  in 
i discussion,  it  is  assumed  that  k*B  - 


The  contribution  of  sorption  nonequilibrium  to  total  nonequilibrium  is 
relatively  insignificant  under  conditions  where  k°  is  greater  than  10, 
and  becomes  virtually  nonexistent  at  values  approaching  100.  A 
bicontinuum  model  representing  only  the  influence  of  PNE  would  be 
adequate  in  these  situations.  The  k°  value  at  which  sorption 
nonequilibrium  becomes  relatively  insignificant  is  a function  of  the 
relative  sizes  of  che  sorption  domains  (i.e.,  the  0 values);  however, 
this  functionality  exhibits  low  sensitivity.  Of  course,  if  0,  * 0,  - 0, 


uoj|ej|uaouo3  eAj|E|ey 


the  model  becomes  totally  Insensitive  to  k°.  We  are,  however,  concerned 
with  cases  where  ^ " A > 0 (i,e.,  when  the  MPNE  model  Is  applicable). 

When  k°  - 0.01,  Che  MPNE  simulation  Is  similar  to  a blconclnuum 
simulation  wherein  Che  value  for  R is  reduced  by  one-half.  For  k° 

<0.01,  the  rate  of  mass  transfer  is  so  slow  that  the  sorption 
nonequilibrium  component  has  minimal  effect  on  solute  transport. 
Sorption  domains  2 and  4 are,  in  essence,  switched  off  and  domains  1 


n of  the  solute  Is  reduced  by  fiz  + 

s (see  Figure  5-1) , its  influence  on  solute  transport  is 
expected  Co  be  dependent,  in  part,  upon  Che  value  of  w.  From  the 
results  of  Che  sensitivity  analysis  lc  appears  chat  for  u equal  to 


that,  for  cases  where  Che  system  is  relatively  insensitive  to  k.a* , 
setting  kla'  equal  to  k_*  should  not  adversely  affect  performance  of  the 
model.  Hence,  under  these  conditions  the  number  of  parameters  for  which 

A sensitivity  analysis  for  u is  presented  in  Figure  5-3.  Similarly 


sorption  nonequilibrium,  TNE  becomes  relatively  insignificant 


UOIJEJJU30U03  aAi)E|oy 


1 representing  only 


the  contribution  of  sorption  nonequilibrium  (not  shown) . As  was  the 
case  for  k , Che  w value  at  which  TNE  becomes  relatively  insignificant 
is  only  slightly  sensitive  to  p.  Ac  the  lower  limiting  value  of  u - 
0.01,  the  behavior  is  again  similar  to  that  observed  for  SNE.  In  this 
case,  however,  the  race  of  mass  transfer  between  the  mobile  and  immobile 
regions  is  so  slow  that  domains  3 and  4 are  effectively  switched  off. 
leaving  1 and  2 as  the  predominant  active  domains. 

Close  scrutiny  of  Figures  5-2  and  5-3  shows  some  differences 
between  the  two  cases.  For  values  of  a less  chan  approximately  10, 

simulations  where  u was  fixed  at  1 and  k° 


deviations  e: 


i 2.  Figure  5-3) 
u decreases . XI 
the  status  of  Che  immobile  ) 
interaction  between  the  immt 
influenced  by  k°.  whereas  it 
1 when  k"  equals  0.01,  tl 
participation.  However. 


The  degree  of  divergence  is  observed  to 
e deviations  result  from  differences  in 

•t  is  influenced  by  u.  For  example,  for  case 
immobile  pore  water  remains  in  active 
>r  case  2 when  u equals  0.01,  the  immobile 
d from  significant  participation, 
of  a nonsorbing  solute  is  often  used  to  assess  the 
d relative  importance  of  TNE.  This  approach  may  result  in 


is  conclusions.  Consider,  for  example,  a 
is  obtained  by  fitting  a bicontlnuum  mode 


151 

(see  eq.  5-22)  for  a particular  sorbing  solute,  the  u value  for  that 
solute  would  be  2.5.  The  conclusion  chat  TNE  is  insignificant  nay  not 

maybe  Used  for  a sorbing  solute  after  correcting  for  differences  in 

which  to  evaluate  the  relative  significance  of  PNE  for  sorbing  solutes. 

The  four  limiting  cases  for  the  MPNE  model  with  respecC  to  the 
nonequilibrium  parameters  (i.e.,  Damkohler  numbers)  are  given  in  Table 

produced  by  the  MPNE  model  for  the  limiting  cases  are  compared  in  Figure 
5-4  to  simulations  produced  by  a model  employing  the  LEA.  Note  chat 
values  of  100  and  0.01  were  used  in  place  of  “>  and  0,  respectively.  The 

model  for  all  four  cases.  For  example,  the  simulation  produced  by  Che 
reduces  to  Che  LEA  model  tt 


4,  in  addition  c< 
true  limicing-cai 


analyses  presented  in  Figures  5-2  and  5-3, 

Lues  of  **  and  0,  respectively.  These  values,  then, 


significance  f 


1.  Limiting 


MPNE  Model 


R.  - R (fit*  ft) 

R.  - R (ft  + ft) 


■e  equivalent  retardation  factor  fc 
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and  of  Figure  5-4  was  to 
MPNE  model  reduces  Co  Che 
le  LEA  model,  respecCively.  The  magnicude  of  Che 
is  used  as  Che  cricerion  for  chis  evaluadon.  As 
le  Damkohler  number  has  been  used  by  ochers  as  an  LEA- 
validicy  erlcerlon.  Resales  of  chelr  analyses  showed  Che  llmldng 
Damkohler  number  value  for  LEA  validicy  Co  be  approximacely  100  (Melnyk 
ec  al.,  1983;  Jennings  and  K&rkner,  1984;  Valocchi,  1985;  Bahr  and 
Rubin,  1987;  Lassey,  1988).  These  resulcs  are  In  agreemenc  wlch  chose 
obealned  herein.  Jennings  and  Kirkner  (1984)  reporced  a Damkohler 
number  of  10  as  being  sufficienC  Co  produce  a "reasonably  good 


from  che  analyses  are  applicable,  in  a relacive  sense,  Co  any  P and  R 
values,  based  on  Che  discussion  presenced  aC  Che  beginning  of  Chis 


illuscraced  in  Figure  5-4.  For  limicing  cases  2,  3,  and  4,  Che  vercical 
porcion  of  Che  BTCs  of  Che  MPNE  model  are  essenclally  coincidenc  with 


whereas  Che  LEA-model  BTCs  do  n 


The  magnicude  of  chis  calling  is 
bove  Che  MPNE-model  BTCs  are  Che  same 
1.  The  close  proximicy  of  che  cail 


•e  generally  Indeterminant  in  these  cases  and 
such  values  result  in  apparent  retardation  factors 
chat  are  lower  than  the  true  value.  These  apparent  retardation  factors 
often  exhibit  a velocity  dependence,  which  confounds  efforts  to  analyte 
sorption  dynamics.  The  use  of  apparent  retardation  factors  thus  may 
lead  to'  various  problems,  especially  when  the  nature  of  these  apparent 

alternative  experimental  techniques  are  required.  One  such  technique, 
called  flow  interruption,  was  presented  in  section  3-5.  It  should  be 


Figures  5-2,  : 


f nonequilibrium  on  solute  transport  is  illustrated  in 
and  5-4.  In  general,  a BTC  is  rotated  clockwise  and 
is  shifted  leftward  under  conditions  of  nonequilibrium,  thus  becoming 
asymmetrical  and  skewed.  When  neither  nonequilibrium  process  is 

model  (limiting  case  1) . As  one  process  becomes  increasingly 
significant  (while  the  other  remains  insignificant),  the  BTC  is 
displaced  leftward  as  predicted  by  a bicontinuum  model.  As  the  second 

simulated  by  a bicontinuum  model. 


To  illustrate  the  use  of  the  MPf 
validity  of  the  model's  conceptual  bj 
several  published  data  sets.  Ail  of 
were  obtained  are  aggregated  soils  wl 


s.  the  model  was  applied  tc 
e systems  from  which  Che  di 
moderate  organic  carbon 


concents.  All  of  the  experiments 
conditions  of  steady  water  flow  u: 

Several  approaches  were  used  Co  obtain  independent  values  for  the  mi 
parameters.  In  some  cases,  the  approach  employed  may  not  represent 


possible  given 
st  alternative 


most  desirable;  however,  the 

the  constraints  of  limited  data.  In  these  cases,  the 
approach  is  selected.  It  is  hoped  that  the  examples 
will  serve  as  a useful  illustration  of  Che  various  strategies  available 


1.  (1977).  They  presented  BTCs  obtained  from 
miscible  displacement  of  the  herbicide  2,4,5-T  through  columns  packed 
with  an  aggregated  soil  (Glendale  loam,  organic  matter  content  - 0.8»). 
The  purpose  of  those  experiments  was  to  evaluate  Che  performance  of  a 
two-region  model  developed  by  van  Genuchten  and  Wlerenga  (1976).  The 
simulations  resulting  from  a three-parameter  (P.  $,  u)  optimization  of 
the  data  compared  well  with  the  experimental  BTCs,  although  not  as  well 
as  did  the  optimized  simulations  of  a nonsorbing  solute  (tritiated 
water)  reported  in  a companion  paper  (van  Genuchten  and  Wlerenga.  1977). 


e obtained  from  the  triciated-wacer  data.  Since  Che 

additional  nonequilibrium  process  is  suggested. 
This  was  noted  by  van  Genuchcen  et  al.  (1977)  and  was  attributed  to 

The  first  data  set  selected  for  use  was  from  Experiment  1-6  (v  - 

re  linearised  Kp  value  obtained  by  vi 
1 batch  isotherm  experiment,  in  which  it  1 
nonlinear  (n  - 0.79,  where  n is  the  Freundlich  exponent). 

nonequilibrium.  Hence,  the  w value  obtained  from  optimization  of  Che 

which  were  Chen  used  to  calculate  ~l  - 0.3  using  eq.  5-22.  0 then,  is 
the  only  remaining  parameter  to  be  determined.  A simplifying  assumption 
of  f - d (Nkedi-Kizza  et  al.,  1986)  was  made  in  order  to  estimate  a 

fitting  the  bicontinuum  model  to  the  BTC  of  a nonsorbing 


0.65  cm/hr,  aggregate 

triciated-water 
(1977)  from  a bs 


breakthrough)  a 


generally  valid,  1c  serves  as  a useful  flrsC  approximation.  The  value 
0 value  estimated  by  assuming  f - 0.93  is  0.93,  where  0 - {0  + f(R-l))/R 
With  the  procedure  used  above,  values  for  all  model  parameters  have 

le  bicontinuum  model  is  presented  in  Figure  5-5.  The 
*C  exhibits  a significant  leftward  shift  (i.e.,  early 
.n  asymmetrical  profile  in  comparison  to  the  model 
prediction.  The  early  breakthrough  and  asymmetry  exhibited  by  the  data 
suggests  the  presence  of  an  additional  nonequilibrium  mechanism.  If 
such  an  additional  mechanism  were  present,  the  MPNE  model  should  provide 
an  improved  simulation  of  the  data.  The  values  for  P,  R,  and  w were 
obtained  for  the  MPNE  model  in  Che  same  manner  as  described  for  the 

was  assumed  that  k‘a  - k‘lra  - k*,  an  assumption  chat  would  seem 
equation  for  the  Type-II  chemicals  was  used,  given  the  chemical  nature 
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i 03  (recall  chat  Che  foe 
r be  used  CO  calculate 

Che  remaining 


The  remaining  unknown  parameters  are  . fi2i  e 

values  for  these  parameters  if  F,  f,  and  i ar 
same  assumption  that  f - d as  used  above,  F i 

the  work  of  Karickhoff  (1980)  and  Karickhoff  and  Morris  (198S)  who 
evaluated  sorption  kinetics  for  several  HOC/soil  systems  with  a 

bicontlnuua  model  shows  that  the  MPNE  model  provides  a much-improved 
one  nonequilibrium  process  is  affecting  the  transport  of 

A multi-parameter  (£, , . fi2,  u,  k“)  optimization  of  the  data  was 

performed  with  the  MPNE  model.  The  optimized  values  for  u and  k«  were 
0.49  and  6,  respectively.  The  value  of  the  optimized  a is  the  same  as 

to  correct  for  differences  in  diffusivities.  The  equivalency  of  the 
calculated  and  optimized  u values  supports  the  validity  of  obtaining  an 

factor-of-chree  larger  chan  the  k2  value  estimated  from  the  k vs  Kp 
regression  equation,  a difference  chat  is  within  the  scatter  of  the  data 


fortuitous,  it  may  also  support  the  suggestion  that  the  sorption-related 
process  involved  in  this  system  is  comprised  of  a rate-limited 
interaction  betveen  2,4, 5-T  and  sorbent  organic  matter.  The  sorption  of 

Anderson,  1974).  If  sorption  is  controlled  by  organic  matter,  it  is 
probable  that  the  sorption-kinetics  mechanism  would  involve  organic 

nonequilibrium  sorption  exhibited  by  2,4-D  (Khan,  1973;  Section  4-5),  a 


from  optimization  of  oi 


aggregate  diameter  s 6 am) 
and  u were  independently  ol 
the  independent  predictions  o: 


*1-  (W77) 


parameters  obtained 


equal  to  k',  four  parameters  require  optimization  (/),,  0Z,  pa,  k").  The 
4-parameter  optimized  simulation  of  experiment  3-3  is  presented  in 
Figure  5-6.  The  k2  value  calculated  from  k‘  is  similar  to  the  kj  value 
calculated  from  the  optimized  k*  value  obcained  from  experiment  1-4  (14% 
deviation).  The  closeness  of  these  two  values  to  that  predicted  with 
the  k2  vs  Kp  regression  equation  is  again  supportive  of  rate-limited 
sorbate-organlc  matter  interaction  as  the  nonequilibrium  mechanism. 

The  parameters  determined  for  experiment  3-3  were  used  to  predict 
the  BTC  obtained  from  a different  column  experiment  performed  under 


aggregate  diameter 
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experiment  3-3.  The  values  for  u and  k*  were  adjusted  to  account  for 
the  slight  differences  in  velocities  between  the  two  experiments.  p and 
R were  obtained  from  the  fit  of  the  bicontlnuum  model  to  the  trltlated 

independent  prediction  is  presented  in  Figure  5-7a.  The  parameters 

experiment  1-4  (see  Figure  5-7b),  a system  with  a velocity  and  an 
aggregate  site  one-third  that  of  experiment  3-3.  Mote  that,  again, 

k*  reflect  adjustment  for  differences  in  velocities.  That  the  same 
values  for  the  nonequilibrium  parameters  (i.e.,  n,  k2 , fit)  can  be  used 
to  predict  two  experiments  with  different  velocities  and  aggregate  sites 
suggests  that  the  nonequilibrium  processes  are  relatively  insensitive  to 

values  associated  with  these  two  experiments.  While  k*  would  be 
expected  to  be  Independent  of  aggregate  size,  a would  be  expected  to  be 
dependent  (see  eq.  5-20).  The  exhibited  insensitivity  vould  thus 
suggest  that  a sorption-related  (i.e.,  kj-cons trained)  process  is 


Hornsby  (1972) . He  Investigated  the  influence  of  aggregate  size  on 
miscible  displacement  of  the  herbicide  fluometuron  through  columns 
pocked  with  Norge  loam  (organic  matter  content  - 1.7%).  Aggregates  of 


prepared  by  Creating 


solution  of  Krilium,  a soil  conditioner.  BTCs  obtained  for  chloride, 

behavior  consistent  with  that  of  a PNE-influenced  system.  Conversely, 
the  BTC  for  Che  column  packed  with  the  smallest  aggregates  (0.25-0.5  mm 

.1-aggregate  * column . 

in  Figure  5-9.  The  nonlinear,  least-squared,  optimization  program 

The  value  for  P was  determined  from  the  chloride  data  and  R was 
determined  by  moment  analysis;  i.e.,  using  the  first,  normalized, 
absolute  temporal  moment.  The  nonequilibrium  exhibited  by  the  small- 
aggregate  data  would  seem  Co  result  from  a sorption-related  process, 
given  chat  the  comparable  chloride  data  exhibited  no  nonequilibrium. 


present  in  the  large-aggregace/chloride  system  and  if  SHE  is  present  for 
the  small-aggregate/fluometuron  system,  it  is  logical  to  assume  chat  the 
fluomecuron  BTC  for  the  large-aggregate  column  is  influenced  by  both 
transport-  and  sorption-related  nonequilibrium.  Hence,  the  need  for  the 

Parameter  values  for  the  MPNE  model  were  obtained  in  the  following 
manner.  The  Peclet  number  was  obtained  by  fitting  a bicontinuum  model 
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co  Che  chloride  BTC  obcained  from  c 
value  was  al90  obcained  from  Che  sa 
factor  (eq.  5-22)  was  calculaced  an 
discussed  previously.  The  value  fo 
decermined  by  moment  analysis.  A v 
llh)  using  the  kj  value  d 


e large-aggregate  column 
e chloride  data; 


»e  fluomecuron/ small-aggregate 


F was  obtained  from  the  fluomecuron/small- 
A value  for  d was  determined  from  Che  BTC  o. 

he  simplifying  assumpcion  of  f 


determined  independently.  The  independent  prediction  ol 


ITC  obcained  from  t: 
lodel  prediccion.  ' 


le  simulated  BTC 


large  aggregates  is  shown 
a is  closely  matched  by  the 
displaced  leftward  of 


the  data.  This  deviation  is  primarily  a reflection  of  the  value  used 
for  Pi  which,  in  turn,  is  a reflection  of  the  assumpcion  made  to 
determine  a value  for  f.  Note  that,  given  the  low  value  of  u.  the 
assumption  of  kj’  - k.‘  most  likely  did  not  influence  the  predicted 


The  final  data  set  utilized  is  chat  of 
Although  their  experiments  were  performed  m< 
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data  set  represents  one  oC  the  most  complete  sets  collected  to  date. 
They  Investigated  the  effect  of  flow  velocity  upon  transport  of  the 
pesticide  lindane  through  columns  packed  wich  an  aggregated  soil 


simulate  lindane  BTCs  with  che  LEA  model 
to  the*  existence  of  intraaggregate  diffusion  (i 
of  flow  rate  on  the  transport  of  chloride  is  si 
BTC  for  the  higher-velocity  (v  - 2.48  cm 

tailing  and  clockwise  rotation  exhibited  by  the  chloride  d. 


s inability 
he  influence 


If  chloride  transport  was  affected  by  PHE  at  the  higher  velocity, 
lindane  transport  should  be  as  well.  Based  on  the  analyses  of  the 
previous  experiments,  it  might  be  expected  that  a sorption-related 


test  this  hypothesis,  parameter  values  for  the  MPNB  model  were  obtained 
independently  as  follows.  The  Peclet  number  and  w were  obtained  from  a 
fit  of  the  bicontinuum  model  to  the  high-velocity  chloride  data;  the  t> 
value  was  corrected  for  differences  In  diffusion  coefficients  as 
described  above. 


(1967) . To  obtain  values  for 

-kinetics  experiment  was  fit  w: 
5-12).  The  value  determined  f< 


>r  was  calculated  using  a Kp  value 
experiment  reported  by  Kay  and  Elrick 
i sorption-kinetic  parameters  F and  k2, 
.ck  (1967)  from  a batch  sorption 

kj  was  used  with  eq.  5-llh  to  calculate 


calculated 


is  determined 


predicted  by 


5-18,  employing  the  assumption  that  £ - 4,  and 
from  the  chloride  data. 

to  the  high-velocity  lindane  BTC  in  Figure  5- 
the  model  describes  the  data  relatively  veil.  Some  of  the  small 
deviation  is  probably  a result  of  an  accidental  change  in  flow  rate  tl 
was  reported  by  Kay  and  Elrlck  (1967)  to  have  occurred  during  the 
experiment.  The  parameter  k°  defined  in  eq.  5-23  can  be  used  to 
evaluate  LEA  validity  for  the  sorption-related  nonequilibrium  process. 
Using  values  given  in  Figure  5-13  for  the  required  parameters,  a valui 
of  52  is  calculated  for  k#.  Recall  from  the  sensitivity  analysis 
section  that  a value  of  100  vas  determined  as  the  limit  at  vhich  the  I 
vas  valid,  and  that  the  MPNE  model  reduces  to  the  bicontinuum  model  it 
this  situation.  Given  the  large  calculated  k°  value,  the  bicontinuum 
model  should  provide  a prediction  very  similar  to  that  provided  by  the 
MPSE  model,  vhich  in  fact  it  does  (BTC  not  shovn).  This  experiment  is 
good  example  of  a case  vhere,  although  tvo  nonequilibrium  processes  ve 
in  operation,  experimental  conditions  vere  such  that  only  one  of  the  t 
had  a significant  impact  on  solute  transport.  If,  however,  another 
experiment  had  been  performed  at  some  larger  velocity  (e.g.,  greater 
than  25  cm/hr) , it  might  be  expected  that  the  bicontinuum  model  vould 
longer  be  able  to  adequately  predict  the  outcome. 


A model  vas  presented  that  explicitly  accounts  for  multiple  si 
of  nonequilibrium.  The  performance  of  this  model  vas  compared  to  i 
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ts  in  additional  terms.  Obviously,  if  values  for  th 

perform  better  simply  because  of  the  greater  numbei 
parameters.  Evaluation  of  model  performance  under 


these  additional 
ideal  condition 


lists  when  the  number  of  fitting  parameters 
, independent  prediction.  Hie  results 


Cher,  it  vas  an  attempt  to  improve  the  description 
of  solute  transport  in  natural  systems  by  introducing  a more  complete 
representation  of  transport/sorption  dynamics.  The  additional 
parameters  introduced  with  the  MPNE  model  represent  actual 
physicochemical  processes  and,  thus,  are  not  simply  "adjustment 
factors".  The  ability  to  obtain  values  for  the  kinetic  parameters  of 

parameters  represent  actual  physicochemical  processes. 

Of  course,  additional  parameters  could  be  produced  with  further 
discretisation  of  the  porous  medium.  This  process  could  be  extended  to 
the  limiting  case  of  n domains  and  n-1  mass-transfer  processes.  Models 
based  on  such  an  approach  may  be  able  to  simulate  experimental  data 
well;  however,  the  number  of  parameters  involved  would  greatly  exceed 


present  capability  to  independently  evaluate  the  processes 


Che  MPNE  model,  Okham's  razor  was  ucllized  In  balancing  the  number  of 

suggest  that  the  relationship  may  serve  useful  In  the  investigation  of 
sorption  nonequilibrium.  The  MPNE  model  was  successfully  used  to 
provide  completely— independent  predictions  of  published  data.  These,  to 
our  knowledge,  are  the  first  such  successful  predictions  for  organic 
solutes.  The  success  of  the  MPNE  model  in  independently  predicting  Che 
data  provides  strong  support  for  the  validity  of  the  conceptualization 

Analyses  of  che  published  daca  would  appear  to  suggest  that  under 

nonequilibrium.  This  potential  would  be  especially  significant  at  the 
field  scale.  Hence,  when  investigating  or  modeling  solute  tl 
is  important  to  consider  the  possible  existence  of  multiple  s 


The  local  equilibrium  assumption  (LEA)  and  Che  assumption  of 
spatial  homogeneity  in  aquifer  properties  are  two  primary  components  of 
many  models  used  to  predict  the  transport  of  contaminants  In 
groundwater.  The  validity  of  these  assumptions  has  recently  become  a 
topic  of  considerable  interest,  especially  in  light  of  data  showing 
deviations  from  that  predicted  with  the  simple  models.  Several 
processes  can  cause  this  nonideal  behavior,  including  physical 
nonequilibrium,  chemical  nonequilibrium,  intrasorbent  diffusion  spatial 
variability  in  sorption  coefficient,  and  spatial  variability  in 
hydraulic  conductivity.  While  a significant  amount  of  research  on 
transport  nonideality  has  been  performed  at  the  laboratory  scale,  little 


Solute  transport  can  be  investigated  employing  either  the  Eulerian  or 
the  Lagrangian  approach.  The  nature  of  the  transport  nonideality 
displayed  is.  in  part,  dependent  upon  which  of  these  two  approaches  is 
employed.  For  example,  it  has  been  well-documented  that  the  temporal 
first  moment  of  a solute  pulse  is  independent  of  kinetic  constraints. 
Conversely,  the  spatial  first  moment  is  influenced  by  kinetic  processes 
(Goltz  and  Roberts.  1987;  Valocchi.  1988).  The  characterization  of 
transport  nonideality,  may,  therefore,  depend  upon  the 


approach 


Involved  either 


second  temporal  moment  or  the  first  spatial  moment.  Under  "ideal" 
conditions  (e.g.,  LEA  valid,  homogeneous  domain)  solute  transport  vill 
bo  characterized,  for  example,  by  constant  retardation  factor  (R)  values 
and  minimal  dispersion.  The  operation  of  nonideal  phenomena  can  result 
in  increased  dispersion  (larger  second  temporal  moment)  and 
time-dependent  R values  (time-variant  first  spatial  moment). 

Solute  transport  may  be  investigated  at  several  different  spatial 
scales,  including  pore-scale,  macroscopic  (e.g.,  laboratory  columns, 
small  portions  of  an  aquifer),  megascopic  (that  associated  with  larger 
portions  of  an  aquifer),  and  gigascopic  (that  associated  with  multiple 
units  or  formations) . The  degree  of  dispersion  exhibited  by  a solute 
pulse  at  any  given  scale  is  a function  of  all  the  processes  contributing 
to  dispersion  at  chat  scale  and  at  all  smaller  scales.  Thus,  the 
relative  importance  of  the  various  phenomena  that  may  contribute  to 
dispersion  is  dependent  upon  the  scale  of  interest.  Applying  the 
adveccive-dispersive  transport  equation  to  a solute  front  at  a given 
scale  results  in  the  determination  of  a dispersion  coefficient  or 
dispersivity  value  (o)  chat  is  associated  with  that  particular  scale  and 
chat  is  representative  of  all  processes  contributing  to  the  dispersion. 
If  the  dispersivity  associated  with  hydrodynamic  dispersion  is  taken  as 


may  be  considered  as  effective  a's.  Generally,  these  effective  m's 
Increase  in  magnitude  with  increasing  scale.  The  Interplay  between  the 
various  length  scales  of  the  system  of  interest  (e.g.,  characteristic 
lengths  associated  with  the  nonidealicy  phenomena,  length  scale  of 


device)  greedy  Influences  Che 


problem, 

observed 


Mechanisms 

se  below,  may  also  cause  nonidealley 


a)  Hydrodynamic  dispersion-  dispersion  in  this  case  results  from 
nonuniform  velocity  distributions  at  the  intra-pore  and  Inter-pore 
scales,  the  tortuous  nature  of  thB  flow  paths,  and  the  existence  of 


solute  concentration  gradients.  Increased  dispersion  can  result,  for 


example,  when  the  porous  medium  has  a wide  pore-size  distribution. 

b)  Physical  nonequilibrium-  physical  nonequilibrium  (PNE)  results 
from  a heterogeneous  flow  domain  ac  the  macroscopic  level.  Regions  of 
lower  hydraulic  conductivity  (K)  within  the  flow  domain  act  as 
sink/source  components,  with  rate-limited  dlffuslonal  mass  transfer 
between  these  regions  and  the  flow  domain  causing  dispersion  of  the 
solute  front.  These  sink/source  regions  can  take  various  forms, 


including  the  internal  porosi 
matrix  of  fractured  media,  an 
typically  found  in  aquifers, 

c)  Chemical  nonequilibrl 


K micro-layers  or  laminae 


sorbent  surface  Is  rate-limited.  As  discussed  by  Brusseau  at 


(1989) , CNE  is  generally  not  in  important  nonequilibrium  mechanism  for 
HOCs  and,  hence,  may  usually  be  excluded  from  consideration. 

d)  Incrasorbenc  diffusion-  nonequilibrium  associated  with 
intrasorbent  diffusion  (1SD)  results  from  the  diffuslonal  mass  transfer 
of  sorbate  from  the  sorbent  surface  into  the  interior  matrix  of  the 
sorbent.  Considering  organic  matter  and  mineral  surfaces  as  the  two 

intraorganic  matter  diffusion  (IOMD)  or  intramineral  diffusion  (IMD). 
IOMD  has  been  proposed  as  a cause  for  sorption  nonequilibrium  by  several 
:e  section  9-1).  IMD  in  systems  involving  clay  minerals 


Involving  typical  aquifer  materials  (i.e.,  quart!  sand)  has  only 
recently  received  attention  (cf . . Ball  and  Roberts,  1987).  Although 
both  PHE  and  ISD  involve  a diffuslonal  mass  transfer  mechanism,  they  may 
be  considered  as  two  different  processes  because  1)  P(JE  Influences  both 
sorbing  and  nonsorbing  solutes,  whereas  ISD  impacts  only  sorbing 
solutes,  and  2)  the  mass  transfer  mechanism  is  quite  different  for  the 
two  processes. 


e)  Vertical  variations  in  sorption  coefficient-  The  existence  of 
vertical  variations  in  sorption  coefficient  (Kp)  can  cause  increased 
dispersion  of  sorbing  solutes.  Consider,  for  example,  a porous  medium 
comprised  of  several  layers,  each  with  a different  Kp,  but  all  with  the 
same  K.  Solute  is  injected  simultaneously  into  all  layers.  As 
transport  proceeds,  differencial  advancement  of  the  solute  fronts  in 
each  layer  will  occur  because  of  the  variations  in  retardation.  This 
differential-front  advancement  will  lead  to  the  observation  of  increased 


dispersion  ac  a depth-integra*  -d  sampling  poinc  down-gradient  of  the 
injection  point.  In  real  situations,  spatial  variations  in  both  K and 
Kp  would  be  expected.  The  existence,  and  nature  of,  correlations 
between  the  two  properties  would  be  important  for  such  cases. 

This  process  has  received  very  little  attention  to  dote.  Smith  at 
Schwartz  (1981)  evaluated  che  Influence  of  s pat  tally-variable  cation 
exchange  capacity  (i.e.,  variable  retardation)  negatively  correlated 


with  K on  solute  transport  in  a simulated  heterogeneous  aquifer.  They 
found  dispersion  to  Increase  as  the  contrast  in  retardation  increased. 
However,  they  concluded  Choc  this  effect  would  be  of  secondary 
importance  in  comparison  to  the  effect  of  spatial  variability  in  K.  The 
influence  of  a spatially-varying  R on  solute  transport  in  a homogeneous 
aquifer  was  investigated  with  model  simulations  by  Sovlch  and  Dougherty 
(1987).  Effective  dispersion  coefficients  exhibited  scale-dependent, 
asymptotic,  oscillatory,  or  negative  values,  depending  on  the 
configuration  of  che  retardation  field.  This  mechanism  was  proposed  as 
the  cause  of  the  enhanced  dispersion  exhibited  by  lithium,  in  comparison 
to  bromide,  in  a recent  natural  gradient  experiment  at  Cape  Cod,  HA 
(Garabedlan  ec.  al.,  1988). 

f)  Macrodispersion-  macrodispersion  results  from  spatial  variations 


of  K in  the  flow  domain  at  the  megascopic  scale  (i.e.,  heterogeneous 
aquifer).  This  process  has  been  the  focus  of  much  research.  Note  that 
macrodispersion  is  similar  to  the  physical  nonequilibrium  that  results 
from  the  existence  of  laminae  of  varying  K,  that  was  discussed  above  in 
(a).  In  fact,  the  mechanism  producing  the  two  phenomena  is  the  same, 


only  difference  being 


differentiated  by  the  relati*  spatial  resolution  of  the  sampling  device 
employed  to  observe  solute  concentrations.  The  Increased  dispersion 

(i.e.,  a monitoring  well  screened  over  a substantial  portion  of  an 
aquifer)  is  considered  macrodispersion.  Conversely,  the  increased 
dispersion  observed  at  the  scale  of  a single,  discrete  sone  within  an 
aquifer  is  considered  the  result  of  PNE.  Obviously,  the  use  of 
discrete,  multilevel  sampling  devices  is  necessary  to  observe  this  PNE. 
Such  devices  have  an  effective  vertical  sampling  sone  of  approximately  8 
to  10  cm  (cf . , Mackay  et  al.,  1986a),'  Hence,  PNE-relaeed  dispersion  in 
aquifer  systems  may  be  defined  as  that  observed  at  the  centimeter  scale 

megascopic  scale. 

The  difference  between  the  two  phenomena  may  be  illustrated  with  the 
data  reported  by  Pickens  and  Grisak  (1981).  A dispersivity  of  0.035  cm, 
representing  the  contribution  of  hydrodynamic  dispersion,  was  measured 
from  miscible  displacement  through  columns  packed  with  aquifer  material. 
At  the  scale  of  a discrete  level  within  the  aquifer,  a dispersivity  of 
0.7  cm  was  measured  for  data  collected  with  multi-level  point  sampling 
wf  « field  tracer  experiment.  The  difference  between  this  value  and  the 
laboratory  value  represents  the  effect  of  PNE.  The  -full-aquifer- 
dispersivity  of  50  cm,  derived  from  a small  two-well  tracer  test  with 
depth-integrated  sampling,  represents  the  effect  of  macrodispersion. 

The  a values  obtained  at  different  scales  may  be  compared  to  assess 


dispersion  mechanisms. 


Grisak  (1981) , and  DeSmedt  and  Ulerenga  (1984)  show  than  PNE—  Induced 
nonequilibrium  can  resulc  in  roughly  an  order-of-magnitude  increase  in 
dispersivity.  Data  from  Pickens  et  al.  (1981),  Hutxler  et  al.  (1986), 
and  Lee  et  al.  (1988)  show  that  sorption-related  nonequilibrium  (e.g., 

dispersivity.  This  order-of-magnitude  effect  is  also  demonstrated  by 
dispersivity  values  determined  at  different  scales  for  two  aquifer 
systems,  Borden  and  Chalk  River.  An  a value  of  0.05  cm,  representing 

addition  to  hydrodynamic  dispersion,  were  obtained  by  fitting  the  A-D 
model  to  the  BTCs  of  two  HOCs  obtained  from  miscible  displacement 
through  Borden  sand,  a values  of  approximately  5 cm  were  determined  for 
BTCs  obtained  at  individual  levels  for  a field-scale  experiment 

dispersivity  represents  the  PNE-inducod  dispersion  associated  with  the 
effects  of  the  micro-layers  that  exist  in  the  aquifer,  as  well  as 
hydrodynamic  dispersion  and  sorption  nonequilibrium-induced  dispersion 
(see,  for  example,  Sudicky  (1986)  for  a discussion  of  the  heterogeneity 
of  the  Borden  aquifer).  Similar  o values  for  the  various  processes  were 
obtained  by  Pickens  and  Grisak  (1981)  and  Pickens  et  al.  (1981).  The 
order-of-magnitude  impact  associated  with  the  macroscopic  dispersion 


processes  (e.g. , PNE,  ISD  Ct  ) could  be  significant  for  systems  where 

vertically-averaged  samples  are  taken,  these  order-of-magnitude  effects 
are  likely  to  be  small  compared  to  the  impact  of  macrodispersion.  For 
example,  "macroscopic"  dispersivlty  values  of  approximately  50  cm, 

natural-gradient  solute  transport  experiment  was  performed,  were 
obtained  by  Sudicky  (1986)  and  Freyberg  (1986) . Full-aquifer 
macrodispersivity  values  of  30-60  m were  determined  from  analyzing  the 
transport  of  "bomb"  tritium  in  the  Borden  aquifer  (Egboka  et  al.,  1983). 

relative  effect  of  macroscopic  vs.  megascopic  dispersion  processes. 

, ISD,  and  ONE  at  the  field-scale  is  supported  by  a 

(1988).  His  analysis  indicated  that  "kinetic  effects  may  not  be 
important  for  modeling  depth-averaged  plume  behavior  in  highly 
heterogeneous  systems."  A similar  conclusion  was  reached  by  Bretz  et 
al.  (1986).  Based  on  the  extrapolation  of  experimental  results  obtained 

processes  at  the  field  scale  is  also  supported  by  the  work  of  Holz  and 


tracer  tests  they  have  evaluated  results  from  vertical  distributions  of 


Eleld-Scale  Transport  Nonldealltv 
A recently  completed  natural-gradient  field  experiment  in  a sand 
aquifer  at  the  Canadian  Forces  Bose,  Borden.  Ontario  (Mackay  et  al., 
1986a;  Roberts  et  al.,  1986;  Curtis  et  al.,  1986;  Sudicky,  1986; 
Freyberg,  1986;  Colts  and  Roberts,  1986)  provides  an  excellent  case 
study  of  organic  contaminant  transport  at  the  field  scale.  As  this  is 
the  first  comprehensive  field  experiment  specifically  designed  to 
eveluate  processes  affecting  the  transport  of  organic  solutes  in  an 
aquifer  system,  valuable  Information  would  be  expected  from  the  results. 
An  evaluation  of  the  results  with  respect  to  sorption  nonldeality  and 
dispersion  is  presented  below.  Data  obtained  with  an  Eulerlan  approach 
will  be  discussed  first,  followed  by  a discussion  of  data  obtained  with 
the  Lagrangian  approach. 


Eul«rian  Apprcach 

Breakthrough  curves  (BTCs)  for  nonsorbing  (chloride)  and  sorbing 
(carbon  tetrachloride,  bromoform,  and  tetrachloroethene  (PCE))  solutes 
obtained  at  a specific  x,y,z  coordinate  were  reported  by  Coltz  and 
Roberts  (1986).  The  authors  fitted  two  bicontinuum  models,  representing 
PNE,  to  the  data.  Inspection  of  their  Figures  6 a,  b,  c and  7 a,  b,  c 
shows  that,  while  the  bicontinuum  models  provided  a generally  good  fie 
to  the  data,  the  tailing  exhibited  by  the  BTCs  of  the  sorbing  solutes 
was  more  extensive  than  that  predicted  by  the  models.  Conversely,  the 
models  did  predict  very  well  the  tailing  of  the  nonsorbing  solute,  which 


. chat  of  the  sorbing  solutes.  This  behavior. 

additional,  sorption-related,  nonequilibriun  mechanism.  Fitting  a model 
that  includes  a single  nonequilibrium  term  Co  data  Influenced  by  two  or 

or  global  kinetic  parameter.  While  this  approach  may  be  adequate  for 
estimating  Che  BTC,  it  is  not  conducive  to  process-level  investigation 

To  further  investigate  the  apparent  nonequilibrium  exhibited  by  the 

developed  in  section  5 was  used  to  simulate  the  BTC  of  PCE  (see  Figure 
6-1) . This  model  accounts  for  the  presence  of  two  nonequilibrium 
mechanisms,  one  associated  with  PNE,  and  a second  that  may  represent  one 

below).  The  u value,  which  is  the  nondimensional  rate  parameter  (i.e., 
Damkholer  number)  representing  Che  effect  of  PNE,  was  determined  by 
fitting  a bicontinuum  model  Co  the  Cl"  BTC.  This  u value  was  adjusted 

Cl*  and  PCE.  The  retardation  factor  value  was  calculated  using  a 
partition  coefficient  determined  from  a batch  isotherm  reported  by 
Curtis  et  al.  (1986).  The  Peclet  number  (P) . representing  hydrodynamic 
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utilizing  Borden  aquifer  material  reported  by  Lee  et  al.  (1988).  The 

associated  with  instantaneous  sorption),  are  fit  values.  The  values  for 
ond  Pk  (fractions  of  rate-limited  sorption)  were  arbitrarily  sec 
equal  such  that  Che  simulation  is,  in  effect,  a two-parameter  fit  (the 
sum  of  the  four  0’s  equals  one). 

From  the  k*  value  of  3,  a dimensional  first-order  rate  constant 
(kj)  of  0.0019  hr-1  is  calculated.  This  value  can  be  compared  to  values 
that  might  be  expected  from  various  sorption-related  nonequilibrium 

0.0019  hr-1,  which  suggests  that  10MD  is  not  an  important  mechanism  at 
the  scale  of  this  experiment.  Curtis  (1984)  reported  a relaxation  time 

aquifer  material  and  PCE.  Using  this  value,  kj  values  of  roughly  0.1  Co 

values  may  represent  the  combined  effects  of  IMD  (diffusion  of  solute 
into  microporous  sand  grains,  as  proposed  by  Ball  and  Roberts  (1987)) 
and  10MD.  Comparing  these  kj  values  to  0.0019  hr-1  shows  that  the 
former  are  one  to  two  orders-of-magnitude  larger.  This  suggests  that, 
while  the  processes  represented  by  the  batch  experiment  (e.g.,  1SD)  may 
contribute  to  the  nonideality  exhibited  by  the  PCE  BTC,  this 
contribution  is  not  significant  enough  to  explain  all  of  the 


nonideali cy . Thus, 
involved. 

As  mentioned  i 


e involved  in  this  system. 


is  vertical  variations  in  Kp.  The  Kp  values  Cor  PCE 
over  nearly  an  order-of-magnicude  among  10  cm  depth 
Interval  of  the  Borden  aquiCer  (Maekay  et  al . , 1986b) . The  existence  of 

variability  in  both  K and  Kp  may  be  responsible  for  the  nonideal 
behavior  exhibited  by  the  data  of  Golta  and  Roberts  (1986) . 

The  value  for  u (0.3),  representing  PNE,  is  an  order-of-magnitu. 

This  suggests  that  PNE  has  a somewhat 


Schwartz  (1981) , discussed 
ocher  factors  (e.g.,  hydrodynamic  dispersion, 
with  the  discussion  presented  in  section  3. 


is  in  agreement 


Synoptic  sampling  of  tl 
also  performed  during  the  Bi 


(Roberts  et 
affinity  of  the  aquifer 


calculated  re 

L.,  1986).  The  temporal  Increase  in  sorption  that  is 
■s  behavior  nay  be  explained  by  one  or  more  of  the 

of  increasing  sorption 
direction  of 


i directional  t: 


.)  resulting  in  an  increasing  R . 
: a result  of  plume  spreading;  3 


>n  decreases  with  time 

nonsingularity  resulting  in  the  existence  of  a 
component;  4.  differential  action  of  transformation  processes  such  that 
is  leading  odges  of  the  plume  are  diminished 
: the  trailing  edges;  5.  PNE;  6.  CNE;  7.  ISD' 


1 variations  ii 


the  average  sorption 
.n  the  horizontal  direction 
irtis  et  al.  (1986)  have  shown  that  the 
compounds  involved  are  approximately  linear. 

)rm  (30  day)  desorption 


and  exhibit  no  significant  hysteresis, 
experiment  performed  by  Goltz  (1986)  exhibited  only  slight  hysteresis 
and  no  evidence  of  the  existence  of  an  irreversibly-sorbed  component. 
There  was  no  evidence  of  greater  solute  disappearance  at  the  leading 
edges  of  the  plume  (Roberts  et  al.  1986).  These  results  eliminate 
alternatives  1 to  4 from  consideration  as  potential  major  contributors 
to  the  observed  non-ideality.  Alternative.  6 and  7 may  be  considered  . 
relatively  insignificant  at  this  scale,  as  discussed  in  the  previous 
section.  Thus,  PNE,  resulting  from  the  spatial  variability  in  K,  and 


vertical  variations  in  Kj,  are  the  only  remaining  probable  causative 
mechanisms  for  the  observed  sorption  nonideality. 

The  existence  of  small-scale  heterogeneity  in  K at  the  Borden  site 
is  veil  documented  (Sudicky.  19B6).  The  influence  of  this  heterogeneity 
on  the  transport  of  nonsorbing  tracers  has  been  investigated  (Sudicky, 
1986;  Freyberg,  1986).  The  enhanced  dispersion  resulting  from  this 
heterogeneity  would  also  affect  the  transport  of  the  sorbing  solutes. 

The  question  is,  is  the  effect  of  this  process  (FNE)  great  enough  to 
explain  all  of  the  nonideal  behavior  exhibited  by  the  sorbing  solutes? 
Goltz  and  Roberts  (1988)  simulated  the  spatial  first  moment  data 
obtained  from  the  synoptic  sampling  with  a model  representing  physical 
mass  transfer  in  a layered  system  (i.e.,  PNE).  If  PNE  were  the  primary 
mechanism  responsible  for  the  nonideal  behavior  (temporal  Increase  in  R) 
exhibited  in  the  experiment,  the  simulations  produced  with  this  model 
should  well-predict  the  data.  In  fact,  the  model  simulations 
underpredicted  the  degree  of  plume  deceleration  (i.e.,  the  magnitude  of 

nonideality  mechanism  is  involved.  Inspection  of  Figure  8 in  Roberts  et 
al.  (1986)  shows  that  the  magnitude  of  the  temporal  increase  in  R 
Increases  with  the  sorptivity  of  the  solute.  This  vould  suggest  that 
the  additional  nonideslity  mechanism  is  sorption-related.  Given  the 
discussions  presented  above  and  in  the  previous  section,  it  appears  that 


m prasenced  above  demonstraces  Chat  soluce  transport 
at  the  field  scale  is  influenced  by  multiple  nonideality  processes  and 
that  the  Impact  of  these  processes  is  scale-dependent.  A model  must 

under  these  conditions.  To  address  this  issue,  a model  that  includes 
multisource,  scale-dependent,  nonideality  was  developed.  The  model 
explicitly  accounts  for  the  following  processes:  hydrodynamic 
dispersion,  physical  nonequilibrium,  sorption  nonequilibrium,  spatially- 
variable  hydraulic  conductivity,  and  spatlally-varlable  equilibrium- 


The  model  is  based  on  che  model  presented  by  Sudicky  <1989).  This 
model  utilizes  solution  in  Laplace  time  domain  to  eliminate  che  need  for 
time  marching.  This  approach  is  much  more  efficient  and  is  much  easier 
to  modify  to  Incorporate  additional  components.  It  Is.  however,  limited 
to  linear  problems. 


s presented  by  Sudicky  (1989)  a. 


for  spatlally-varlable  hydraulic  conductivity  (K). 
with  one  of  two  approaches:  deterministic,  where  a K value  ii 
for  each  element,  or  stochastic,  where  statistical  parameters 
representing  the  K field  are  employed. 

The  model  of  Sudicky  (1989)  was  modified  to  incorporate  spatlally- 
varlable  equilibrium-sorption  constant  and  mulciprocess  nonequilibrium; 
che  formulation  presented  in  Figure  5-1  was  utilized  for  the  latter. 

The  modified  model,  then,  can  account  for  all  of  the  dispersion 
mechanisms  discussed  in  the  previous  section.  This  model  should  be  very 
useful  for  investigation  of  the  relative  importance  of  che  different 


t solute  transport  al 


The  scale  effect  for  nonideality  can  be  illustrated  with  the 
following  two  examples . For  the  first  example  data  for  the  transport  of 
strontium  ( Sr)  through  a model  microaquifer  reported  by  Starr  et  al. 
(1985)  was  used.  The  microaquifer  consisted  of  three  horizontal  layers, 
where  che  middle  layer  was  a sand  and  the  top  and  bottom  loyers  were  the 
same  silt.  A model  accounting  for  advecclve  transport  in  the  sand 
layer,  diffusive  mass  transfer  in  the  silt  layers,  and  local  equilibrium 
governed  retardation  in  all  layers  was  used  to  simulate  the  experimental 
MSt  BTCs  Starr  et  al.,  1985).  The  authors  could  obtain  a fair 

le  data  only  by  optimizing  the  values  for  R;  the  R value 


for  the  silt  was  much  smaller  chan  that 
From  this,  it  may  be  suggested  that  a l 
sorption  nonequilibrium,  is  involved  ii 
variability  of  hydraulic  conductivity  i 

Figure  6-2).  Values  for  P (50)  and  w ( 


aird  nonideality  process,  namely 
addition  Co  che  vertical  spatial 
id  of  tho  sorption  coefficient. 

>.07)  were  obtained  by  fitting 


a chloride  BTC;  a value  for  R for  the  entire 
porous  medium  (l.e.,  weighted  average  R)  was  obtained  using  the  Kp 
values  obtained  from  the  batch  isocherms;  the  values  for  0,  + 0 , (0.2' 
and  03  + 0A  (0.B)  were  based  on  the  mass  fractions  of  che  two  types  ol 
materials  (sand,  silt).  This,  in  essence,  is  employing  the  assumptlo 
that  the  sand  layer  Is  serving  as  a mobile  domain  and  that  che  silt 
layers  are  acting  as  Immobile  domains.  Hence,  the  mobile-immobile  ma 


r represented  by  u 


e spatially-variable 
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hydraulic  conductivity.  The  influei 
sorption  coefficient  (Kp  - 1.65  for 
>r  by  modifying  tl 


« spatially-variable 
and  12.5  ml/g  for  the  silt) 
of  w vith  the  factor 


t for  kj  calculated  from  k* 
a of  0.001  hr-1  reported 
oxyhydroxides  (Melnyk  et 
suggest  that  chemical 


id  example,  simulations  Vi 


b produced  with  the  modi 
‘here  the  realization  ui 


presented  in  this  section.  The  case  used 
heterogeneous  hydraulic  conductivity  fiel 
produced  with  the  turning-bands  technique  (Sudlcky,  1989).  The 
statistical  properties  employed  to  develop  the  realization  vere  based  on 
values  reported  by  Sudlcky  (1986)  for  Che  Borden  aquifer.  Sorption  was 
constrained  to  be  rate-limited;  the  values  for  K,  and  k:  were  selected 
to  represent  those  obtained  for  the  sorption  dynamics  of 
tecrachloroethene  on  the  Borden  aquifer  material  (see  the  discussion 
associated  with  Figure  6-1).  The  relative  importance  of  nonequilibrium 
sorption  was  investigated  by  performing  two  simulations,  one  with  local 
equilibrium  assumed  and  one  with  race-limited  sorption  included  (see 
Figure  6-3). 

Comparing  the  two  plumes,  it  appears  chat  Che  inclusion  of 
nonequilibrium  sorption  did  not  significantly  impact  plume  development. 
This  is  not  surprising  given  chat  the  characteristic  time  of  reaction 
used  for  this  example  is  smaller  chan  the  effective  residence  time  of 


(ui)  aouDjsig  _ 2 


Che  solute.  Thus,  sorption  nonequilibrium  may  not  be  an  important 
process  affecting  the  transport  of  solute  under  relatively  low 
velocities  (i.e. . natural  gradient  conditions)  at  the  field  scale. 
Is  not  to  say,  however,  that  nonequilibrium  may  not  be  important  at 
field  scale  under  high-velocity  conditions,  such  as  those  occurring 
during  the  flushing  of  contaminated  aquifers. 


SUMMARY 

In  Chapter  3,  experimental  methods  for  investigating  the  sorption 
kinetics  of  organic  chemicals  were  developed  and  evaluated.  It  appears 
that  the  gas-purge  technique  is  ideal  for  investigating  the  sorption 
dynamics  of  volatile  organic  chemicals  in  scdlment/water  systems  and  in 
high  organic  carbon  soil  systems.  The  design  used  herein  has  limited 
usefulness,  however,  for  low  organic  carbon  systems  (e.g.,  aquifer 
materials).  The  miscible  displacement  technique  seems  well  suited  for 
investigating  sorption  and  transport  in  low  organic  carbon  systems.  The 
gas-purge  and  miscible  displacement  techniques  were  shown  to  provide 
comparable  results  for  the  determination  of  kinetic  parameter  values. 
These  two  techniques  may  therefore  be  used  to  complement  each  other;  the 
choice  of  which  to  use  being  determined  by  the  nature  of  the  system  of 


A modified  form  of  the  miscible  displacement  technique,  called  flo< 
interruption,  was  developed.  This  technique  should  be  useful  in 
enhancing  che  robustness  of  parameter  value  determination.  The  flow- 
interruption  method  may  also  prove  useful  in  such  large-scale 
applications  as  remediation  of  contaminated  aquifers.  The  methods 
presented  in  Chapter  3 were  employed  to  investigate  the  mechanisms 

Three  objectives  were  pursued  in  Chapter  4.  The  first  objective 
was  an  investigation  of  che  influence  of  sorbate  molecular  structure  on 


sorption  kinocics.  A log-log  Inverse  relationship  was  discovered 
between  the  kinetic  sorption  constant  and  the  equilibrium  sorption 
constant;  this  relationship  was  quantified.  This  relationship  was  shown 
to  be  useful  in  providing  estimates  of  values  for  the  sorption  rate 


The  Influence  of  sorbate  molecular  structure  on  sorption  kinetics 
was  recognized  to  be  important  under  certain  conditions.  The  sorption 
of  so-called  nonhydrophobic  chemicals  (e.g. . chose  containing  reactive 
functional  groups)  was  found  to  be  much  more  constrained  than  that  of 
the  hydrophobic  chemicals.  This  was  suggested  to  result  from  specific 


can  significantly  affect  sorption  kinetics.  Conversely,  the  Influence 
of  molecular  structure  on  kinetics  appears  to  be  of  secondary  importance 
for  hydrophobic  chemicals.  For  example,  the  kj-K,,  relationship  for  a 
homologous  series  of  benzene,  naphthalene,  and  anthracene  was  found  to 
be  essentially  the  same  as  chat  exhibited  by  two  other  series  (di,  tri, 
and  tetra  chloroechenes  and  chloro,  dlchlora,  trichloro,  and  cecrachloro 


An  investigation  of  the  Influence  of  organic  cosolvent  on  the 
nonequilibrium  sorption  of  organic  solutes  was  the  second  objective  of 
Chapter  A.  A model  describing  a log-linear  relationship  between  the 
sorption  race  constant  and  volume  fraction  cosolvent  was  developed. 
Experiments  were  performed  to  evaluate  this  model;  che  results  appeared 
to  support  the  validity  of  the  proposed  model.  It  was  also  shown  that 
results  obtained  from  experiments  performed  with  mixed  solvents  could  be 


systems . The  use  of  mixed  solvents  should  prove  very  useful  in  reducing 

Identification  of  the  mechanism  responsible  for  sorption 
nonequilibrium  was  the  third  objective  of  Chapter  A.  Diffusive  mass 

was  suggested  to  be  the  most  probable  cause  of  sorption  nonequilibrlua 


organic  carbon  content  of  a soil  waj 
y supportive  of  i 


•is  hypothesis.  For  example,  the 
.sm  responsible  for  rate-limited 

involving  diffusive  mass  transfer 

nonequilibrium.  The  results  of  the  calcium  and  quinoline  experiments 
are  also  supportive  of  the  involvement  of  this  type  of  mechanism. 

A model  for  simulating  solute  transport  under 

Chapter  5.  This  model  will  be  useful  for  investigating  and  modelling 
solute  transport  in  complex  systems  where  nonequilibrium  may  be  caused 
by  more  than  one  process.  It  will  also  be  useful  for  elucidating  the 


nonideality  nay  result  fron  several  processes  (hydrodynamic  dispersion, 
physical  nonequilibrlua,  sorption  nonequllibriua,  spatially-variable 
hydraulic  conductivity,  and  spatially-variable  sorption  constant)  was 
presented  in  Chapter  6.  This  nodel,  along  with  the  nodel  presented  in 
Figure  5,  was  used  to  simulate  several  cases  of  nonideal  solute 
transport  that  have  been  reported  in  the  literature.  The  results  of 
these  simulations  suggest  chat  the  process(es)  responsible  for 

laboratory),  nonequilibrium  sorption  may  be  the  predominant  nonideality 

heterogoneous  soil/aquifer  properties  (e.g. , hydraulic  conductivity, 
equilibrium  sorption  constent)  may  predominate. 
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